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CHAPTER I 


INTRODUCTION 


Antenna radiation and target scattering parameters, characteristics and pat- 
terns are measured on an instrumented facility called a “range” [1]. In the case 
where the scattering from targets is to be measured, it is required that the scat- 
tering body be illuminated by a plane wave of constant amplitude. Plane wave 
illumination is similarly required to measure the far field radiation parameters of 
an antenna. The measurement range thus has to be able to produce a uniform 
plane wave with no cross-polarization errors (or a good approximation thereof) 
over a specified volume where the target or antenna will be located, referred to as 
the “target zone”. A local plane wave in the target zone can be created in two 
ways; viz., 

1. Locating the target zone so that it is in the actual far zone of the source 
antenna[2,3]. In order to meet the requirements for a plane wave approxi- 
mation one finds that the actual distance separating the source antenna and 
the target zone becomes very large, so that the measurement range can only 
be implemented outdoors. 

2. Using a collimator to transform the spherical wave radiated by a point source 
into a plane wave[4]. This conversion can be accomplished in a relatively 
short space, so that these types of ranges can be accommodated indoors, 
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and are referred to as “compact ranges”. Indoor compact ranges are nor- 
mally housed in so called anechoic chambers, in order to simulate a free 
space environment [5]. In true free space all energy radiated from a source 
would propagate to infinity since there is per definition nothing to scatter 
the energy back. The walls, floor and ceiling of an anechoic chamber are con- 
sequently lined with electromagnetic absorbing material, so that (ideally) all 
electromagnetic energy incident thereon is absorbed rather than scattered. 
Two types of collimators have been proposed; viz., parabolic reflectors[6] 
and lenses[4]. The lens concept has not proven to be very succesful, due to 
the fact that metallic lenses are very frequency dependent and the manufac- 
turing difficulties that are associated with large dielectric lenses. Parabolic 
reflectors, on the other hand, have been very succesful and are widely used 
in modern compact ranges. Compact ranges offer many advantages over 
other types of ranges, and as a result a lot of effort is being directed towards 
the improvement of their performance. For a general discussion on compact 
ranges one is referred to references [4], [T]— [10] . A typical focus-fed compact 
range configuration is shown in Figure 1. 

This study proposes several design improvements that can be made to the 
antenna system of a compact range to enhance the system performance; viz., 

1. Shaping the edge contour of the parabolic main reflector. 

2. Adding blended rolled surfaces to the edges of the parabolic main reflector. 

3. A dual chamber configuration in which the main reflector and target zone are 
contained in one chamber and the feed assembly together with a Gregorian 
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Figure 1: Typical focus-fed compact range configuration. 

subreflector in the other. This study is mainly concerned with the absorber 
fence separating the two chambers. 

The first two suggested improvements are aimed at reducing the diffracted fields 
from the main reflector. These diffracted fields interfere with the plane wave and 
degrade system performance in that they create a ripple in the reflected plane wave. 
A dual chamber configuration has several advantages. The use of a tilted Gregorian 
subreflector reduces the amplitude taper in the reflected field considerably, and 
virtually eliminates the cross-polarization errors. Since the main reflector and 
target zone are isolated from the feed assembly and subreflector by the absorber 
fence, minimal scattering from the subreflector or spillover from the feed illuminate 
the main reflector or target zone. The suggested design improvements will thus 
enable one to design a compact, range with very small amplitude taper, ripple and 
cross-polarization errors in the illumination of the target zone. 
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Chapter II discusses the techniques that are used to analyze the antenna and 
dielectric structures in the rest of this study. A method, whereby the end point 
contributions resulting from the application of the physical optics method to three- 
dimensional bodies can be eliminated, is derived. Chapter III discusses the criteria 
for and creation of plane waves and it is shown how rolled edges in general improve 
the field quality in the target zone. Elliptical and blended rolled edges are analyzed 
and compared for two-dimensional cases. The relationship between the order of 
a blending function and the resulting continuity in the derivatives of the surface 
radius of curvature across the junction is derived. An algorithm based upon this 
relationship is proposed for designing optimized blended rolled edges. Design ex- 
amples are given which clearly show the superiority of blended rolled edges over 
elliptical rolled edges. Three-dimensional parabolic main reflectors are discussed in 
Chapter IV. It is shown that a substantial reduction in the edge diffracted fields can 
be achieved if the edge contour of the parabolic section has a concave shape rather 
than the conventional convex shape. The design equations for three-dimensional 
main reflectors with concave edge contours and blended rolled edges are given. Sev- 
eral design examples are discussed to show that these design improvements result 
in increased performance in the case of three-dimensional reflectors. Chapter V 
discusses the benefits of a Gregorian subreflector feed system and proposes a dual 
chamber compact range configuration in which the use of an absorber fence pre- 
vents unwanted fields that are diffracted from the subreflector edge and spillover 
from the feed from reaching the main reflector and target zone. In Chapter VI a 
technique is described to compute the forward scattered field from a homogeneous 
lossy dielectric slab with a wedge termination. This type of structure is used in the 
implementation of the absorber fence in the dual chamber configuration. Chap- 
ter VII contains the summary and conclusions that can be made from this study 
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and suggests further topics that can be studied. Two appendices are included, 
dealing with applications of differential geometry and the derivation of closed form 
expressions for the aperture integration over a half plane in free space. 
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CHAPTER II 


ANALYTICAL TECHNIQUES 


This chapter briefly describes the existing techniques that are used to ana- 
lyze the reflector and absorber structures in this study, with the exception of the 
moment method. Since the derivation of these solutions are not the purpose of 
this section, the equations relating to Geometrical Optics (GO) and the Uniform 
Theory of Diffraction (UTD) are simply stated without a detailed discussion, and 
the reader is referred to the references for more details. A new technique to calcu- 
late the forward scattering from a homogeneous lossy dielectric slab with a wedge 
termination is discussed in Chapter VI. 

The moment method is based on a field integral equation solution, and is 
assumed to be exact within the constraints of numerical accuracy and in so far 
as the scattering body can be accurately modelled. A detailed understanding of 
the moment method is not required for the analysis of structures discussed in 
this study, nor for the physical interpretation of results. The reader is referred to 
references [11,12] for details on the theory and application of the moment method. 

Throughout this study an e JU>i time convention is assumed and suppressed. 
All reflector surfaces are assumed to be perfectly conducting and smooth. In 
two-dimensional problems, transverse electric (TE) polarization indicates the case 
where the electric field vector is perpendicular to the plane of incidence, and trans- 
verse magnetic (TM) polarization indicates the case where the magnetic field vector 
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Figure 2: Astigmatic ray tube. 

is perpendicular to the plane of incidence. Certain symbols are used throughout, 
including the free space impedance ( Zq ), free space admittance (Vo), free space 
wavelength (A), and free space propagation constant (k = Vector quanti- 

ties are indicated by boldface symbols (eg. E), unit vectors by boldface symbols 
with hats (eg. n) and dyadic quantities by boldface symbols wih an arrow on top 
(eg. R). 

2.1 Geometrical Optics 

Geometrical Optics is a ray tracing procedure in which the wavelength is 
allowed to become infinitesimally small, ie., a high frequency limit is assumed. A 
ray optical field[13] can be described in terms of an astigmatic ray tube as shown 
in Figure 2. The electric field at l associated with this ray tube is given by 

E(/) = E(0 )A[ Pi , P2 ,l)e-> u (2.1) 
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where E(0) is the field at the reference (0), A is a spread factor and pi,P2 are the 
radii of curvature of the wavefront. The spread factor is given by 




I PI P 2 
( Pi + 1){P2 + 0 


( 2 . 2 ) 


and accounts for the conservation of energy in a ray tube. Note that there is a 
congruence of the rays at the so called “caustics” (1,2 and 3,4 in Figure 2), so 
that the caustics can be viewed as virtual sources of the wavefront. Examination 
of Equation (2.2) shows that the GO field is unbounded at the caustics (l = 
—p\,—p2)- Although Equations (2.1) and (2.2) describe a three-dimensional field, 
the two-dimensional version can easily be derived by setting the appropriate radius 
of curvature equal to infinity. 

The ray paths follow Fermat’s principle [14] which requires that the path tran- 
sit time be an extremum. It is found that the electric and magnetic fields are 
perpendicular to each other and the ray path. 

In applying Equations (2.1) and (2.2) to the field reflected from a surface, one 
finds that the reflected field can be expressed as[15] 


E r (r^ p ) = E*(r rp ) • R 


Pi P 2 


,-jks r 


(P\ + 


'M 


+ s r ) 


(2.3) 


where 


E r (r /P ) = 

E (r rp ) - 

r _ 

Pi, 2 ~ 


reflected electric field at the field point (r f p ) 

incident electric field at the reflection point (r rp ) 

the principal radii of curvature of the reflected 

wavefront emanating from r rp 

distance from r rp to r f p 

unit vector normal to the surface at r rp 
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s' = direction of incident propagation at r rp 
s = direction of reflected propagation at r rp 


R a ,h = acoustic soft/hard reflection coefficients 


(TE/TM polarization) 

R = dyadic reflection coefficient = i2jej_e^ + e r 
s x n 
e± ~ I s' x n| 

Cy = s 1 x ej_ | and 

efj = s x ej_ . 


(2.4) 

(2.5) 

( 2 . 6 ) 
(2.7) 


Note that R s h. = ^ 1 in the case of a perfectly conducting surface whose radius of 
curvature is large in terms of a wavelength. In this case, one finds that 


s = s ( - 2(n • s^h , and 
E* • R = -E* + 2(n • E*)n . 


( 2 . 8 ) 

(2.9) 


The incident ray, reflected ray and unit vector normal to the surface at the point of 
reflection all lie in the same plane (plane of incidence)[15], and the law of reflection 
requires that the incident angle be equal to the reflection angle. A reflected field 
at Tf p will only exist if a reflection point (r rp ) on the surface exists which satisfies 
these requirements. 

2.2 Wedge diffraction 

The ray optical GO field discussed in the previous section only exists in those 
regions where the line of sight between the source and the target point is not 
obstructed. Consider a wedge with an angle (WA) illuminated by a line source as 
shown in Figure 3. Three regions can be identified; viz., 
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Region 1: Incident and reflected GO fields present. 

Region 2: Only incident GO field present. 

Region 3: No GO fields present. 

These regions are bounded by the so called “incident shadow boundary” (ISB) and 
“reflection shadow boundary” (RSB). The GO incident and reflected fields exist 
only within these boundaries, and are zero outside. However, electromagnetic fields 
must be smooth and continuous everywhere, so that the discontinuities across the 
shadow boundaries are not allowable. The GO incident and reflected fields are thus 
not capable of representing the true total field in a region where scattering bodies 
are present. Keller developed the Geometrical Theory of Diffraction[16] (GTD) to 
supplement the GO fields and correct these deficiencies. The GTD postulates a 
diffracted field that emanates from the edge of a scattering body, so that the to- 
tal field (ie., GO plus diffracted fields) is continuous. However, Keller’s diffracted 
fields are not uniform since they are unbounded along the shadow boundaries. 
Kouyoumjian and Pathak corrected the GTD when they developed the Uniform 
Theory of Diffraction[17] (UTD). They found that when Keller’s diffraction co- 
efficients are multiplied by a transition function, the total fields are continuous 
everywhere, including across the shadow boundaries. A brief description of the 
two-dimensional UTD will now be given. 

The high frequency UTD first order diffracted field from a point (Qe) on the 
edge of a two-dimensional wedge is given by 

v d M = vim > d, a (2.io) 

v 3 

where U z and U\ represent the diffracted and incident transverse electric (soft) and 
magnetic (hard) fields respectively, D s are the soft/hard diffraction coefficients, 
and s is the distance from Qe to the field point. If a two-dimensional line source 
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generates the incident field, then the two-dimensional diffraction coefficients and 
associated parameters are given by 


D tth {4>, 4>', L\ L rn , U ° , R, th , n) = D 1 + D 2 + R s<h (D 3 + D 4 ) (2.11) 


where 


D i 


D 2 


d 3 


d 4 

a±(/J) 

n 

R»,h 


cot ( z+Sl) 

2ny/2irk \ 2n ) 

~ e - j */ 4 

2ny® l 2n / 

-e-M* t (* + P + \ 

2 nVSri l 2n ) 

( *-(3 + \ 

2nV2irk \ 2n / 


F[fcl i a+ (/?")] 
F[fc£*a _ (/?“)] 
F[ibi7 n a + (/?+)] 
F[fcL ro a- (/?+)] 


2 (2mrN ± - (3\ 

2 cos ( ^ ) 

„ WS4 

2 , and 

7T 

soft/hard reflection coefficient of the surface at Qe • 


(2.12) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 
(2.1T) 

(2.18) 


The parameters and angles are indicated in Figure 4. Note that are integers 
which most nearly satisfy the equations 


2imN + — (3 = 7r , and (2.19) 

2imN~ — fl = — 7r . (2.20) 

The transition function is given by 

_ l r°° ■ 2 

F(x) = 2j^e )kx e~ JT dr (2.21) 

Jy/x 

and is plotted in Figure 5. Note that | F(x) | < 1 in the transition region (x < 2 tt) 
and F(x) & 1 outside the transition region (x > 27t). If x < 0, then one finds that 
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1 wedge diffraction geometry. 
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Figure 5: Wedge transition function. 


F(x) — F*(|.r|) where * indicates the complex conjugate. The transition function 
is quite suitable for numerical evaluation[18]. The rest of the relevant parameters 
are as follows: 


AV # 

3 +3 
J TO _ Po 3 

Po + s 

J_ _ 1 | 2 

Po s' a Q cos 0* 

L rn — EjL s 

Pn + 5 

i_ _ I + 2 

Pn s> a n COS 0^ 

a Q — radius of curvature of the o-face at Qg 
a n — radius of curvature of the 7i-face at Qg 
Po n = principal radii of curvature of the reflected 
wavefront at Qg in the plane of incidence 
from the o and n faces, respectively 


( 2 . 22 ) 

( 2 . 23 ) 

( 2 . 24 ) 

( 2 . 25 ) 

( 2 . 26 ) 
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<t> = the incident angle with respect to the o-face 

4> — the diffraction angle with respect to the o-face 
s' — the distance between the source and Q g 
s — the distance between the field point and Qg 
0 l o = the angle defined by n 0 • s / — — cos 0 l o 

with cos 0 l o >O (2.27) 

0 l n = the angle defined by n n • s' — — cos 0 l n 

with cos 0' n > 0 , and (2.28) 

ho,n = the unit vectors normal to the 

o and n-faces at Qei respectively . 


The transition regions around the incident and reflected shadow boundaries are 
shown in Figure 6. Note that these shadow boundaries occur at. <j> ^ $ — 7 r, 
respectively. The diffracted fields have maximum amplitudes along the shadow 
boundaries, and are strong inside the transition regions. Outside the transition 
regions the diffracted fields decay rapidly. 

As shown in Equation (2.10), the first order diffraction solution is directly 
dependent on the magnitude of the field incident on the edge of the wedge. This 
expression may lead one to believe that if a null field is incident on the edge, then 
the edge diffracted field is zero. This is not the case, however, in that the slope of 
the incident field will create currents on the wedge surface which result in a slope 
diffracted field[19,20,21]. This diffraction is known as “slope diffraction”, and is 


given by 


U?(s) 


1 dD sJ) dUl , e~J k 
jks' d <!>' ' d<f>' (QE) ^ 


(2.29) 
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Figure 6: Transition regions around the shadow boundaries. 
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where dD a ^/ d<j)' are the slope diffraction coefficients for the soft and hard cases, 
respectively, and dUl/dfi is the derivative of the incident field in the direction of 
the edge at Q E as shown in Figure 7. From Equation (2.29) it is clear that the 
magnitude of the slope diffracted fields are proportional to the slope of the incident 
fields. The general form of the slope diffraction coefficients are given by 


_ csc 2 ( Sr ') F,[ * :iio ' (r)1 

-R,,h esc 2 +(/?+)] 

- CSc2 (^) F >l* r °“" ( ^ + )] } 


where /?* and a* are given in Equations (2.16) and (2.17), and 


F a {x) = 2jx[l - F(z)] . (2.31) 

2.3 Physical Optics 

2.3.1 Theory of physical optics 

An electromagnetic wave that impinges upon an object induces electric and 
magnetic currents on that object. These induced currents flow in accordance with 
Maxwell’s equations and the prevailing boundary conditions and in turn act as the 
sources of a secondary field known as the “scattered field”. The Physical Optics 
(PO) method obtains the scattered field by summing the contributions of all the 
current elements. For example, the field radiated by a reflector antenna can be 
calculated by the PO method if the currents induced on the reflector surface by 
the incident illumination from the feed can be found. 
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Consider now the magnetic field resulting from the vector potential (A), which 
is given by [22] 

H = V x A (2.32) 

where 

•Mr') 


A(r ^ - // 

J s = induced surface current , and 


(2.33) 


r f p = field point 


In the PO format, the integration is carried out only over the illuminated portion 
of the surface. The far zone scattered electric field is then given by [23] 

E *( r /p) ~ ~j kZ 0 A* (2.34) 

where the far zone transverse vector potential A* is the component of A which is 
transverse to the direction of propagation; ie., 

A< = A - (A • r*)r s , and (2.35) 

f® = direction of scattered propagation . 


Note that this definition for the electric field requires that the field point be in the 
far field of the source current. 

In general, the current induced on the surface is given by 


J, = nxH ( 


(2.36) 


where 


n — unit vector normal to surface , and 
H* = total magnetic field on the surface . 
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The total magnetic field on the surface can be approximated by making a tangent 
plane approximation; ie., the surface fields are assigned the values that they would 
have if the body had been perfectly smooth and flat at the surface patch of inte- 
gration. This tangential approximation restricts the application of the PO method 
to cases where the surface radius of curvature is large compared to the wavelength; 
ie., there is a lower frequency limit. Since it is assumed that the reflector surface 
is perfectly conducting, it follows from the application of the boundary conditions 
at the surface and image theory that the following conditions hold on the surface: 


n x E < = 0 (2.37) 

n x H r = nxH 1 (2.38) 


where 


E* = total electric field on the surface 

H l = incident magnetic field on the surface , and 

H r = reflected magnetic field on the surface 


so that 


J,= 


2 n x H* in the illuminated region 


(2.39) 


(0 in the shadow region . 

This is known as the “physical optics approximation”, and is valid when the ra- 
dius of curvature of the surface is large compared to a wavelength. The incident 
magnetic field on the surface is taken to be the geometrical optics field. 

Consider now an illuminated body as shown in Figure 8. Let the magnetic 
field incident on the body be a plane wave which is given by 


H‘(r) = H 0 


(2.40) 
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Figure 8: Illuminated portion of scattering body used for PO calculation. 
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where r is a position vector and k the direction of propagation of the incident field. 
The illuminated and shadow regions on the surface are determined by 


< 0 

in the illuminated region 


= 0 

on the shadow boundary , and 

(2.41) 

>0 

in the shadow region . 



Note that the PO approximation assumes that J s is zero in the shadow region; 
ie., there is an abrupt discontinuity in the assumed value of J a across the shadow 
boundary, although the true current is smooth and continuous across this bound- 
ary. This discrepancy introduces an error in the result obtained from the PO 
method which creates false scattering centers at the incident shadow boundaries. 
A way to compensate for this error will be discussed in Section 2.3.3. A similar 
problem occurs in bodies with edges, in that the PO assumption for J s near the 
edges is not correct. The Physical Theory of Diffraction (PTD) [24] has been 
developed to compensate for this error. 

2.3.2 Scattered field from a surface patch 

In order to calculate the scattered fields of a reflector antenna using a PO 
approach, one normally models the surface of the antenna by small planar patches. 
Triangular patches are useful, since the three corners of such a patch lie in a flat 
plane, and a unique vector normal to the patch can thus be constructed. This is 
not always possible with patches that have more than three corners. The patches 
are considered to be very small in terms of a wavelength, so that it can be assumed 
that each patch is illuminated by a local plane wave, and that the field point is 
in the far zone of the patch. The total field at a field point is then obtained by 
summing the contributions of all the illuminated patches on the reflector surface. 
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Let a patch be located in its own local coordinate system as shown in Figure 9, 
with r ; the position vector of a point on the surface of the patch in local coordinates. 
Consider now a triangular patch in the far zone of a point source. The incident 
field on the patch is considered to be a local plane wave, so that it can be expressed 
as 


where 


v ’ R l 


(2.42) 


R { 

r‘ 

e‘(M) 




distance from the feed to the local origin on the patch 
unit vector in the direction from the feed to the local patch origin 
polarization and magnitude of the incident electric field 
independent of R l , and 

angular dependence of the incident electric field . 


Since the incident electric field is assumed to be a local plane wave, the incident 
magnetic field on the patch is given by 

H*(r' ) = y 0 r* x eV) = (2.43) 

where 

h‘ = r 1 x e‘ . (2.44) 

Making the PO approximation, the induced surface current on the patch is 
given by 

•Mr') = 2 n(r') x H l (r') = x h * . ( 2 . 45 ) 

K l 

The vector potential at the field point is then given by 

A(r ^> = // dS ' < 2 - 4 «) 
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SOURCE AT 
FOCUS OF 
MAIN REFLECTOR 


Figure 9: Surface patch coordinate system. 
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where the integration is performed over the entire surface of the patch, and 

R a = distance from the local origin to the field point , and 
r 4 = unit vector in the direction from the 
local origin to the field point . 



am = vector describing the length and orientation of 
the m-th edge of the patch, arranged tip to tail 
in a counter-clockwise direction 
r m = position vector of the midpoint of the m-th edge 


w 

— i .- 9 i . 1 

= r — r 


(2.50) 


n x w 



p 

— 1 ~ 


(2.51) 


|n x w| 


n 

= unit vector normal to the plate = 

aj x a 2 
|ai x a 2 

(2.52) 


T = length of the projection of w onto the 
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sinc(x) = 


(2.53) 

(2.54) 


plane of the patch = |w - (w • n)n| , and 
siu® 
x 

Note that the closed form expression above can be used for a flat plate with an 
arbitrary number of edges, as long it is illuminated by a plane wave and the field 
point is in the far zone of the patch. 

The far zone transverse vector potential is thus given by 

A ‘ (r ^ ) = ^- >HS ' +R ^ x h ‘ ~ x h ‘ < 2 - 55 > 

so that the resulting scattered electric field at the field point can be expressed as 

ikiR<+R V (2 ' 56) 

where 

e 5 = n x h* — (n x h* • r s )r 4 . (2.57) 

The total scattered electric field at the field point is then given by the contributions 
from all the surface patches making up the illuminated part of the surface; ie., 

E total( r fp ) = X] ® patch( T fp ) ' (2.58) 

reflector surface 

2.3.3 Subtraction of the end point contributions 

The abrupt termination of the induced surface currents at the incident shadow 
boundary causes the creation of false scattering centers which will corrupt the true 
pattern of the scattered field. The end point contributions have been known to 
exist for some time. However, most authors have assumed that the end point 
contributions from curved surfaces are insignificant in that the PO integrand goes 
to zero for grazing incidence and scatter. Gupta and Burnside[26] have shown 
that this is not the case, however. In order to demonstrate how this error can be 
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corrected, consider a 3-dimensional circular cylinder with radius (a) and length 
(/) as shown in Figure 10. Let the cylinder be illuminated by a plane wave. The 
electric field incident on the cylinder (E*) is then given by 


E*(n) = e -i*(* l +?‘-n) e ;( n ) 


(2.59) 


where 

e *( n ) 

IV 

n (^n? z n) 
n (^m z n) 


polarization and amplitude of E* at n 

distance from the far zone source phase reference to the origin 
position on the cylinder, or 

a cos <f> n x + a sin <f> n y + z n z , and (2.60) 

incident field propagation direction, or 

—(sin 0i cos <t>ik + sin sin <^y + cos 0jz) . (2.61) 


The incident electric and magnetic fields are related by 


H* = r 0 f* x E* 


such that the magnetic field incident on the cylinder can be expressed as 


(2.62) 


H*(n) = Yq e ~ jkRt e jk l a sin 0 i cos (<A»-<M+*n cos x e »j _ ( 2 . 63 ) 

The current induced on the surface of the cylinder (J 4 ) is given by the PO approx- 
imation as 


J s (n) = 2n(^ n ) x H*(n) 


(2.64) 


or 


J s (n) = 2Foe n(<^ n ) x (r* x e*)e^f as * n0 i cos i < i > i~ < t > n)+z n cos9-\ (2.65) 
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Figure 10: Finite circular cylinder geometry. 
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where the unit vector normal to the surface of the cylinder (n) is given by 

n ((j>n) = cos <j>rik + sin <f > n y . (2.66) 

Notice that n(<^ n ) x (r* x e*) is independent of z n . 

The vector potential at the far zone field point is given by Equation (2.33), 
so that 

A(r s , il s ) = JJ (2.67) 

where the integral is taken over the illuminated part of the cylinder, and 

R a — distance from the origin to the far zone field point 
r a (O s ,(f ) 3 ) = scattered field propagation direction, or 


r a {9 3 ,<f> 9 ) = 

sin 0 S cos <f) 3 x + sin 0 9 sin <j> s y + cos 0 s i , and 

(2.68) 

dS n = 

cl dzfi • 

(2.69) 

Thus, the vector potential is given by 



A= 

2itR s * 

(2.70) 

where 



Iz 

= f'l/2 ejkZn(C ° S 9i+C ° S 6,) dz " = 1 sine (k^) 

l 

= —(cos Oi cos 0 3 ) 

(2.71) 


(2.72) 

T a 

_ ^2 x ^igjfco/3 ) an( j 

(2.73) 

P 

= [sin 0 S cos (4> s ~ 4>n) + sin^j cos(<fo - <f> n )\ 

(2.74) 

The angles (<fr n li ( t ) n2) define the illuminated part of the cylinder. 

Note that in 


some cases, eg., axial backscatter from a body of revolution, it is advantageous to 
treat the finite cylinder as an infinitesimal cylinder and then evaluate /. in the 
appropriate coordinate system. 
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The scattered electric field in the far zone is given by Equation (2.34) such 

that 

E s « -j kZ 0 A * . (2.75) 

Let e a (<j> n ) be the component of n(</> n ) x (r l x e l ) that is transverse to r a ; ie., 

e s (<£ n ) = n(<^ n ) x (r* x e*) - [n(^„) x (r* x e‘) • r a ]r a . (2.76) 


The scattered electric field is then given by 

-jkae-M*^) I+h 


E a = 


2irR s 


(2.77) 


where 


f^ n2 ^^^jkalsinesCos^a-^+sineiCos^i-ipn)]^^ _ (2.78) 

J<t>n 1 


The integral (1^) can be evaluated by the method of stationary phase [27], 
which can briefly be described as follows. Consider the integral given by 


/#) = I*" 

J< Pn\ 


(2.79) 


where 


k 

Fo(<f>n) 


real variable 
real, positive and large 
slowly varying amplitude function that is 
well behaved in the range of integration , and 
$(</>«) = phase term that is real and continuous with 

continuous derivatives in the range of integration 


For large k , I<f,{k) can be asymptotically approximated by 

I<f> ~ Isp + Ie(<f>n 2 ) + h{4>n\) 


(2.80) 
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where 


lap = stationary point contribution , and 
I e ( 4>n\, n2) = end point contributions . 


The stationary point contribution is given by 

t _ Fn( <}>«.)<[-- 2ir J{ k *(4>n*)+%sgn*"(<t>n 3 )} 
s P 0{M )I k\$"(<j> n3 )\ 

where the stationary point (at <f> na ) satisfies 

= *'(M = o • 

d(p n 

The end point contributions are given by Gupta and Burnside[26] as 

oo 

/«(*) = E 

m—0 

where 

Fm(<f>n) = TT~ 
d(pn 

Applying the expressions above to I one finds that 


Fm—lifin) 


[ jkVitn) 


for m > 1 




1 <Pn 2 




(2.81) 


(2.82) 


(2.83) 


(2.84) 


F 0 {<t>n) = e'(^„) (2.85) 

$(^ n ) = a[sin 9 a cos (<j> a - <f> n ) + sin cos(<fo - <f> n )\ (2.86) 

$'(<t>n) = a [sin 6 a sin(</> s - <f> n ) -f sin0jsin(<fo - <f > n )] , and (2.87) 

$"(<M = ~*(<M • (2.88) 


Note that the vector notation merely serves to combine three scalar equations into 
one vector equation. 

The numerical computation of the PO integral as described in Section 2.3.2 
yields 1^ as expressed in Equation (2.80). The “true” scattered field is represented 
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by the stationary point contribution and can be obtained by subtracting the end 
point contributions from 1^, since the end point contributions represent the false 
scattering centers at the incident shadow boundaries. For the purpose of this 
derivation it is assumed that the shadow boundary of interest occurs at </> n 2 - The 
first two terms of the end point contributions (I® and Ig) are determined next. 
a) First term (m=0) 

The first term of the end point contribution is given by 


Wn 2) = 


;**'(*» 2 ) 


(2.89) 


or 


o WM sin cos ( < /’«- < An2)+ sin 0; «*(<£,— <£ n2 )] 

T 0/ J \ _ e yyn2; e (2.90) 

n2 jka[sh\0 a sin (<f> 3 - <f> n 2 ) + sin#; sin(<fo - <j> n z)] 

The first term of the end point contribution at (f > n 2 to the far zone scattered 
electric field (Eq) is thus given by 

o-jks 


E5(»2) = 


(2.91) 


where 


and 


do(<£n2) = 


DS - Izdl 


■e s (<f>n 2 ) 


n 2 = 
s' = 
s' = 


s = 


27r[sin0, sin($ 4 - 0 n 2) + sin0j sin (<f>{ - (f> n 2 )] 
n(<£n 2 > 0 ) = 0 cos<£ n2 x + asin^ n2 y 


(2.92) 


(2.93) 

(2.94) 


distance from the far zone source phase reference to 112 , or 
iZ* — asin0jcos(<^ — <}> n 2 ) , and (2.95) 

distance from 112 to the far zone field point, or 
R s - asin0 s cos{4> s - <j> n 2 ) • (2.96) 
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The term d e will be referred to as the “modified scattering coefficient”. In 
the case of a finite scattering cylinder, one thus finds that 

n e r , v -Z sinc[y (cos 6 a + cos ^)]e a (<^ n2 ) 

0 n2 27r[sin^ s sin(<^ s - <^ n 2 ) + sin^j sin(<^j - 0 n 2 )] 

Note that this notation casts the end point contribution of the scattered electric 
field in terms of an end point scattering coefficient (Dq), evaluated on the surface of 
the cylinder, at the shadow boundary in the a3 n 2/n _ pla«e • The scalar components 
of the end point contribution of the scattered electric field are given by 

®oV n 2 ) = E5(n 2 ) ■ p = e-> k ‘‘ D' 0 (M • p (2.98) 

where p represents each of the coordinate unit vectors; ie., x, y, and z. It is 
evident from Equations (2.91) and (2.97) that Dq (and therefore Eg) is independent 
of the cylinder radius (a). 

In general, let the axis of the cylinder be orientated along a unit vector (e). 
The previous equations can then be expressed in vector notation by using the 
following relationships: 



cos = — r l • e , and 

(2.99) 


cos 9 a = r a • e 

(2.100) 

so that 


cos 0 a + cos Si = e • (r* — r l ) . 

(2.101) 

Furthermore, since 

0 < Si a < 180°, one finds that 



sin#j = | r* x e | , and 

(2.102) 


sin# s = | r s x e | . 

(2.103) 
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Let the component of r l that is perpendicular to e and directed away from e be 
given by 


^ = -(** - (r* • e)e) 

(2.104) 

or 


.j r* - (r* • e)e 

| r* — (r* • e)e | 

(2.105) 

Similarly, let the component of r s that is perpendicular to 

e and directed away 

from e be given by 


rj_ = r s - (r s • e)e 

(2.106) 

or 


r s - (r 3 • e)e 
r - L | r 3 — (r s ■ e)e | 

(2.107) 

where the unit vectors and rj_ are shown in Figure 11. Furthermore, one finds 

that 


sin (<f) s - <j>n 2 ) = (n (<j>n2) x r j_) • e 

(2.108) 

sin(^j - <f> n 2 ) = 2 ) x r j_) * ® 

(2.109) 

cos{<f> s - <j> n2 ) = f j_ • n(^> n2 ) , and 

(2.110) 

cos(</»j - <f> n 2 ) = r± • “(<^n 2 ) • 

(2.111) 

The end point scattering coefficient can thus, in general, be 

expressed as 

D e 0 = I z d e 0 

(2.112) 


where 

cWn2) = 


-e*(^n2) 


2tt(| r‘xe| (n(</>„ 2 ) x r *_) • e+ | r s x e | (h(<^ n2 ) x r*_) • e) 

In the case of a finite scattering cylinder, one thus finds that 

-l sinc[^e • (F - r ? )]e*(<ft w2 ) 

0 n2 ' ~ 2tt(| r‘ x e | (n(</>„ 2 ) x r l ± ) • e+ | r 3 x e | (n {<f> n2 ) x r^_) • e) 


. (2.113) 


. (2.114) 
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Figure 11: The .T n y n -projection of the scattering cylinder. 
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The vector representation of Equation (2.114) is very convenient from a computa- 
tional viewpoint, since the vectors need not necessarily be known in the cylinder’s 
local coordinates. 

Note that the incident field on the cylinder is considered to be a plane wave. 
If, however, the far zone source is a point source (as in Section 2.2) and the incident 
field on the cylinder is considered to be a local plane wave, then e* can be replaced 

by 

e i = f i (2.115) 

where e l a is the angular dependent part of e*. In this case, the far zone source 
phase reference will be the location of the point source. 
b) Second term (m—1) 

The second term of the end point contribution is given by 

-Fi(^ n2 )e>‘*W”2> 






where 


such that 


Fi(<M 


Wn) 


l^n.2 


d(t> n [jk$'(<j) n )\ 


I'Mnl) 


fc 2 ($') 3 


l<f>n 2 


(2.116) 


(2.117) 


(2.118) 


The second term of the end point contribution at <f > n 2 to the far zone scattered 
electric field (Ej), is thus given by 

e -jks 


E5(n 2 ) = e-> k ‘ DUM'- 


(2.119) 


where 


Df = I z df 


( 2 . 120 ) 
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and 


> -i\ gWj+giMl . 

1 ” 2 27rfca[sin 9 a sin(<£ s - <f> n2 ) + sin sin (<fo - <j> n 2 )] 3 

In the case of a finite scattering cylinder, one finds that 

_ -jl sinc[y(cos 9 S + cos ^)][B'(^> n2 ) + B (<f> n2 )} 

1 n2 27rfca[sin0 s sin(<£ s - <f> n 2 ) + sin0j sin(<fo - <^ n 2 )] 3 

where 

B'(<£„2) = I sin ^ sin (^« “ 4 > n2) + sin^j sin (fa - <f> n2 )] — ~ 

d<p n 

and 


( 2 . 121 ) 


( 2 . 122 ) 


(2.123) 


B(<£ n2 ) = [sin^ s cos(^» s - (f> n2 ) + sin^cos(0i - <j> n 2)W i^) • (2.124) 

In terms of a general cylinder with its axis orientated along a unit vector e, the 
scattering coefficient can be expressed as 


Df = I z d\ 


(2.125) 


where 

dUM = 


-WiM + BiM) 


2-rrka(\ r 1 x e | (n(^„ 2 ) x r* ) • e+ | r» x e | (n(<£„ 2 ) x r* ) • e) 3 


(2.126) 


In the case of a finite scattering cylinder one thus finds that 

-jl sinc[ye ■ (r s - r*)]^^) + B(<ft n2 )] 


D!(^. n2 ) = 


2?rA:a(| r* x e | (h(<f> n2 ) x r^) • e+ | r J x e | (n(<^ n2 ) x r^) • e) 3 

(2.127; 


where 


B'(<£n2) = (I r‘xe | (h(<f> n2 )xr l ± )-e+ \ r s x e | {h{<t> n2 ) x r]_)-e) de (2.128) 

d(pn 
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and 


B (<f> n2 ) = (| ? x e | (ri{<j> n2 ) • r*j> | r* x e | (h(<t> n2 ) • ri))e a (^ n2 ) . (2.129) 


If P(<£) and Q(<£) are vectors, then 


^< p - Q > - ^■ Q + p '^’ and 
>xQ) = £ x Q + p x^ . 


(2.130) 

(2.131) 


It thus follows from Equation (2.76) that 


de* 


dn 


x (r* x e 


.• / dn 

>-U x(r 


x e*) • r‘ 


d(j>n d<t> n 

since r* x e* is independent of <f> n . Differentiating Equation (2.66) yields 

dn 


(2.132) 


d<j>n 


— — sin + cos <f>ny • 


(2.133) 


In terms of the general cylinder with its axis orientated along a unit vector e, this 
term can be expressed as 


dn 

d<j>n 


= e x h((j> n ) . 


(2.134) 


The vector dh/d<f> n is a unit vector, since e and n are perpendicular to each other. 
Note that Df, and thus E|, are functions of the radius of curvature of the cylinder 
(a). The larger the radius of curvature becomes, the smaller the contribution from 
the second term will be. In fact, one finds that the higher order terms in the series 
given in Equation (2.83) are functions of (1 /ko) n . Since A: is a large number, higher 
order terms in the series thus become less significant. 

In the case of axial backscattering from bodies of revolution, it was found that 
the first term of the end point contribution (Dg) goes to zero, but that the second 
term (D^) gives the correct expression for the false scattering from the incident 
shadow boundary [28]. 
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The total end point contribution at the incident shadow boundary contour on 
a reflector illuminated by a point source at the focus (E end P ° int ) ca n be found by 
evaluating 

E endpoint = / je e jkce.(i* -P) dc ( 2 . 135 ) 

J C 3 

in a counter-clockwise direction, where 

C = incident shadow boundary contour on the surface 
d e = modified end point scattering coefficient. 
r l = direction of incident propagation on C at the element (dc) 
r s = direction of scattered propagation from C at the element (dc) 
s = distance from the element (dc) on C to far zone field point , and 
s' = distance from the far zone source phase reference 
to the element (dc) on C . 

In the case of a numerical computation where the three-dimensional surface 
is modelled by triangular patches, the finite scattering cylinders are fitted onto 
the sides of the patches making up the incident shadow boundary contour. Fig- 
ure 12 shows how the finite scattering cylinders are fitted onto the incident shadow 
boundary contour of the rolled edge of compact range main reflector. E en dpoint j s 
then found by evaluating 

... -jfc(s4-s') 

E endpomt = £ ( jje + D e } (2>136) 

C 3 

where 

D e = end point scattering coefficient 
s = distance from the shadow boundary on the cylinder 

to the far zone field point (see Equation (2.96)) , and 
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s' — distance from the far zone source phase reference 

to the shadow boundary on the cylinder (see Equation (2.95)) . 

The PO solution in which the end point contributions have been eliminated will 
forthwith be referred to as the “corrected PO method”. 

This analysis of the PO method has been presented in a three-dimensional 
coordinate system. A similar analysis can be done in two dimensions, in which 
case the two-dimensional Green’s function is used in the PO integration, and the 
fields have a 1 /\/R amplitude dependence[26]. 
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FINITE SCATTERING 



FINITE SCATTERING 



(b) Finite scattering cylinders along the entire 

shadow boundary contour of a blended rolled edge. 


ure 12: Location of the finite scattering cylinders at the incident shadow 
boundary of a rolled edge of a parabolic reflector. 



CHAPTER III 


EDGE TREATMENT FOR THE MAIN REFLECTOR 
3.1 Plane wave criteria 

Due to the nature of their operation, one is usually interested in the far zone 
characteristics of radar and communications antennas, although in some special 
cases the near zone characteristics are of interest. Since there are only transverse 
field components present in the far zone, it is mathematically convenient to ap- 
proximate the far zone fields by a plane wave; ie., the spherical phase front radiated 
by an antenna is approximated locally by a planar phase front. This can only be a 
local approximation, since true plane waves require infinite energy and can there- 
fore not exist in reality. If all the field components in the planar phase front have 
the same amplitude, then the plane wave is known as a “uniform plane wave”. 

Consider an Antenna- Under- Test (AUT) being operated in the receive mode. 
Measurements relating to the far zone characteristics of the AUT are based on 
the assumption that it is illuminated by a uniform plane wave of given polariza- 
tion. The definition of radar cross-section (RCS) also assumes that the target is 
illuminated by a uniform plane wave. Although true plane waves do not exist, it 
is nevertheless convenient to use them as a measurement standard. In practice, 
a portion of a spherical or cylindrical wavefront is considered to be a local plane 
wave if the phase and amplitude deviations over that portion of the wavefront do 
not exceed specified limits. 
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Figure 13: Plane wave approximation over an aperture. 

Consider now a point source radiating a spherical wave as shown in Figure 13. 
A plane wave is required over an aperture of width ( D ), that, is removed a distance 
(R) from the point source. Let this plane wave be approximated by the spherical 
wave, such that the maximum phase deviation over the aperture does not exceed 
nax- It is required to find the minimum separation ( R m i n ) where this will 
occur. The phase of the spherical wave ($) at a distance (r) from the point source 
is given by 

${r) = kr (3.1) 

so that the phase at the center of the aperture is given by 

<P{R) = kR (3.2) 

and the phase at the extreme end of the aperture is given by 

$(R + d) = k(R + d) . (3.3) 
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Note that d is the radial distance discrepancy between the spherical and planar 
phase fronts at the extreme end of the aperture as shown in Figure 13 and is given 


by 


R + d = sjR 2 + (D/2) 2 

(3.4) 

so that 


d = R(\j\ + {D/2R) 2 - 1 ) 

(3.5) 

or 


D 2 

d & — — for D <C R . 
8R 

(3.6) 

The maximum 

phase deviation over the aperture is thus given by 



x — kd 

(3.7) 


so that 


R 


min 


i rD 2 

4A A $moa: 


(3.8) 


Due to the quadratic nature of Equation (3.6), the resulting phase error is some- 
times referred to as a “quadratic error”. Clearly, the larger R is, the smaller 
A 4* m ax will be, yielding a better plane wave approximation. Table 1 shows vari- 
ous values of R m i n versus A4> maiS obtained from Equation (3.8). Note that in the 
case of RCS measurements, the phase error is included in both the incident and 
scattered paths, so that 


R = 7rI)2 

min 2AA$ 


(3.9) 


max 


Traditionally an amplitude deviation of 1 dB and a phase deviation of 


7T 


A$max — g (—22.5°) 


(3.10) 


in a phase front have been considered as an acceptable approximation for a plane 
wave[29], resulting in 


Rr 


2D 2 

~T~ ' 


(3.11) 
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Table 1: R m { n versus A $max- 


^4* m ax 

■^mtn 

45° 

d 2 /x 

22.5° 

2D 2 /\ 

11.25° 

4D 2 /X 

10° 

4.5D 2 /X 

9° 

5D 2 /X 

7.5° 

6D 2 / A 

5° 

9D 2 /X 

4.5° 

iod 2 /x 


Recent studies [30] have shown, however, that although these criteria may be ac- 
ceptable for measuring the patterns of antennas with -25 dB sidelobes, the mea- 
surement of low RCS targets and antennas with low sidelobes (-30 dB and lower) 
require that measurements be made at distances much larger than that given 
in Equation (3.11), indicating a phase deviation much smaller than 22.5°. The 
quadratic errors result in the following discrepancies[31]: 

• Reduced main lobe amplitude. 

• Increased sidelobe levels. 

• Filling in of nulls. 

In order to accurately measure low RCS targets and the patterns of antennas 
with low sidelobes, it is thus imperative that the amplitude and phase deviations of 
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the plane wave in the target zone of a compact range he as small as possible. The 
design improvements proposed in this study are aimed at reducing these deviations. 

3.2 Creation of a plane wave by a compact range 

A plane wave can be created in the target zone of a compact range by illu- 
minating a parabolic reflector with a point source located at the focus. This is a 
result of the following geometrical characteristics of a paraboloid: 

1. The two caustic distances for the reflected fields are both equal to infinity [32], 
so that all the reflected rays are parallel to the axis of symmetry. This implies 
that there is no amplitude decay after reflection. 

2. The path lengths of all rays that originate at the focus and are reflected by 
the parabolic surface to a plane that is perpendicular to the axis of symmetry 
are equal. This plane can thus be viewed as the phase front of a plane wave. 
Since all the reflected rays are parallel, this planar phase front will propagate 
into the target zone. It will be shown that for a plane located at H > fc 
(where f c is the focal length of the parabolic reflector), the path length is 
equal to the constant zt + f c - 

Notice, however, that due to the spatial attenuation associated with the point 
source at the focus, the rays illuminating the reflector have a 1/R amplitude de- 
pendence before reflection. Since all points on the surface are not an equal distance 
away from the focus, this spatial decay results in an amplitude taper in the re- 
flected field. Although the reflected field is a plane wave, it is thus not a uniform 
plane wave. This effect is shown in Figure 14(a). Compact ranges usually employ 
offset main reflectors in order to reduce aperture blockage by the feed. The pri- 
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mary feed antenna at the focus is then tilted to decrease the amplitude taper in 
the reflected field[33,34]. 

The reflected field exists as a plane wave only over the surface of the parabola, 
and stops abruptly at the surface termination. Since electromagnetic waves must 
be smooth and continuous everywhere, there will be a diffracted field emanating 
from the terminating edge to compensate for the discontinuity in the reflected field 
and make the total field (ie., reflected plus diffracted fields) smooth and continuous. 
This diffracted field interferes constructively and destructively with the plane wave, 
causing unwanted amplitude and phase ripples in the total field in the target zone. 
The edge diffracted fields are shown in Figure 14(b). In fact, the antenna and/or 
target sizes that can be measured in presently available compact ranges are directly 
related to the edge treatment of the main reflector. Burnside and Peters [35] have 
recently shown that in order to measure low backscatter targets, a ripple of less 
than 0.1 dB is required; whereas, a taper of less than 1 dB is satisfactory. Phase 
deviations of 5° or less have also been suggested[36]. 

In order to relate the magnitude of the ripple in the total field to the level of 
the diffracted field, consider Figure 15. Let e<, e go and e d denote the total, GO 
and diffracted fields, respectively. The maximum deviations of e< from e go will 
occur when e go and e d are either completely in phase or completely out of phase. 

Case 1: e go and e d are in phase 
The ripple amplitude is given by 

R\ — (201og e t - 201og e go ) dB (3.12) 

or 

Rl = 20 log(l + 10 (E d -Ego)/20 ) dB ( 3 . 13 ) 
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TARGET 




(a) Taper due to spatial decay of feed radiation. 


EDGE 

DIFFRACTION 


I — ■ — 



(b) Stray signal due to edge diffraction. 


Figure 14: Top view of a compact range, showing the reflected plane wave and 
diffractions from the edge terminations. 
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TOTAL FIELD 



Figure 15: Ripple amplitude related to the diffracted field level, 
where 

E<i = amplitude of the diffracted field (dB) , and 
Ego = amplitude of the GO field (dB) . 

Case 2: e go and e j are out of phase 
The ripple amplitude is given by 

R 2 = — 201og(l - i 0 (^-^o)/2°) dfi ( 3>14 ) 

Since R 2 > R\ for a given value of E j - E go , the ripple amplitude (R) will be 
defined as 

R = -201og(l - lO^rf - ^? 0 )/ 20 ) dB . (3.15) 

Notice that this definition implies that 

E g o — R < Ei < E g o + R (3.16) 
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Table 2: Ripple amplitude ( R ) versus (E^ — E go ). 


Ed. ~ E go (dB) 

R (dB) 

-10 

3.30 

-20 

0.83 

-30 

0.27 

-40 

0.09 

-50 

0.03 

-60 

0.01 


where 


Et = amplitude of the total field (dB) . 

Several values for R versus (E^ — E go ) are shown in Table 2. 

In order to eliminate the edge diffracted fields completely, one would have to 
use an infinite paraboloid. However, since this is not practical, some form of termi- 
nation of the parabolic section has to be found such that edge diffractions into the 
target zone are minimized. Various options to reduce the energy that is scattered 
from the edge termination of the main reflector into the target zone have been 
investigated in the past, including the use of absorber material as well as serrated 
and rolled edges[33,37]. Serrated edges serve to disrupt the caustic effect of a con- 
vex edge, and are being used on various compact range main reflectors [38,39,40]. 
It was found, however, that the addition of a large curved surface to the termi- 
nation, as shown in Figure 16, can increase the performance of the main reflector 
considerably [41]. The addition of such a curved surface to the parabolic reflec- 
tor reduces the energy diffracted from the termination of the parabolic section, 
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REFLECTED FIELD 
FROM ROLLED SURFACE 


Figure 17: Continuity of reflected field using rolled edges. 


since the rolled edge creates a reflected field which goes more smoothly from the 
parabolic surface to the rolled surface (as shown in Figure 17) than in the case of 
an abrupt termination with no edge treatment. Since the surfaces on either side 
of the edge junction are not the same, there will still be a discontinuity in the re- 
flected field across the junction. However, the discontinuity in the reflected field is 
smaller, and the diffracted field weaker, than in the case with no edge treatment. 
The rolled surface is added in such a way that the surface slope is continuous 
across the junction. Furthermore, the rolled edge continuously curves away from 
the focus and continues into the shadow region. This ensures that the GO rays 
that are reflected from the rolled edge do not enter the target zone. 

An elliptic rolled edge was added to the reflector at the Ohio State Univer- 
sity ElectroScience Laboratory, resulting in greatly improved performance[41,42]. 
In gaining experience with rolled surfaces, the following niles of thumb were 
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developed [32]: 


1. The radius of curvature of the rolled surface at the junction between the 
parabolic reflector and the rolled surface should be greater than a focal 
length. 

2. The minimum radius of curvature of the rolled surface should never be less 
than a quarter of a wavelength at the lowest frequency of operation (A maa: /4). 

The lowest frequency of operation of a compact range main reflector is thus related 
to the minimum radius of curvature of the rolled edge. The upper frequency limit 
is generally dictated by the surface tolerance of the reflector[43]. 

In addition to the adding of rolled surface terminations, it can be shown that 
shaping the edge contour separating the parabolic section from the rolled surface 
terminations will also enhance the performance of the main reflector. Such a 
shaping has also been proposed by Kelleher[44] in conjunction with the control of 
sidelobes from parabolic antennas. 

3.3 Elliptic rolled surface terminations 

3.3.1 Reduction of diffracted fields using elliptic rolled edges 

In order to examine the effect of elliptic and blended rolled surface termina- 
tions on the total field in the target zone, the two-dimensional case will be discussed 
in the rest of this chapter. 

Consider an elliptic rolled edge attached to a parabolic reflector at Pj(yj,Zj) 
as shown in Figure 18. Let the focal length of the reflector be f c . The ellipse has 
semi-major and semi-minor axes of lengths (a e ) and (6 e ), respectively. The ellipse 
is attached in such a way that the surface slope at the junction is continuous. The 
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parabolic surface can be described by 


z = 


4/c 


(3.17) 


where the 2 -axis is the axis of symmetry. In anticipation of the three-dimensional 
case that will follow, the y-axis is chosen as the vertical axis. The coordinates 
of the junctions located at the top and bottom of the parabolic section of the 
reflector, respectively, are then given by 


z 3 = 


4 fe 


(3.18) 


where 


Vj = Vtop indicates the top junction , and 
Vj = ybot indicates the bottom junction . 

In the case of the offset reflectors considered in this study, it is assumed that y^ot 
and yt 0 p are both greater than zero. 

It is convenient to define a local coordinate system for the ellipse such that its 
major axis is aligned with the * e -direction, and its minor axis with the y e - direction. 
Let the origin of this coordinate system be at the junction ( Pj ). In order to make 
the slope of the surface continuous across the junction, the tangential vectors to 
the parabola and ellipse at the junction should be parallel. Let the y e -axis thus be 
defined by 

Ye = — junction (3.19) 

where ivj unc ti on is the unit vector normal to the parabola at the junction, as given 
in Appendix A. The ellipse coordinate axes are then given by 

x e = x p 2 y + x p 3Z , and 
y e = y V 2Y + y P 3z 
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(3.20) 

(3.21) 



where the direction cosines are given by 


Vp 2 

yj 

yjy) + 4/| 

(3.22) 

y P 3 

— 2/c 

s]y) + 

(3.23) 

x p2(ytop,bot ) 

= Ty P 3 , and 

(3.24) 

x p3(ytop,bot) 

= ±y P 2 • 

(3.25) 


This will ensure that x e is parallel to the vector that is tangential to the parabola 
at the junction. The transformation from the * e y e -system to the yz-system is used 
to specify a point on the rolled surface in the yz-system. It is given by 


y(i) ' 

__ 

/ 

x p2 

\ 

y P 2 

*e(7) ^ 

+ 

( \ 
yj 

k *( 7 ) J 


K X P3 

yp3 ) 

, Ve( 7) j 


\ z j > 


where the parametric angle 7 of the rolled surface is valid in the region 


(3.26) 


0 < 7 < 7m 


such that 


(3.27) 


7 = 0 at the junction , and (3.28) 

7 = 7m at the extreme end of the rolled surface . (3.29) 

The ellipse can be described in the x e j/ e -system by the following curve: 

r e(7) = x e{l)*e + y e {l)y e = («e sin 7)3^ + b e (l - cos7)y e . (3.30) 

In terms of the yz-system, the elliptic rolled surface can thus be expressed as 

r ellipse( 7) = Vellipsei 7 )y + z ellipse(l (3.31) 
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where 


V ellipse (7 ) = (o e sin7)x p2 + M 1 ~ cos 7)Vp2 + Vj , and (3.32) 

z ellipse (t ) = («e sin 7 )x p3 + 6 e (l - cos 7 )y p3 + zj . (3.33) 

There is no need for a z e -component, since the ellipse lies in a flat plane. 

It is shown in Appendix A that the radius of curvature of the parabola at the 
junction is given by 

r / \ 21 § 

Rr ato %j) = 2/c 1 +{§-) (3-34) 

and that the radius of curvature of the ellipse at the junction is given by 

■ (3.35) 

It is thus apparent that, although the slope of the surface is continuous across the 
junction, the radius of curvature there is discontinuous. In fact, the parabola is a 
concave surface (when viewed from the focus); whereas, the ellipse is a convex sur- 
face. The radii of curvature on either side of the junction can thus be interpreted 
as having different signs in addition to being discontinuous in magnitude. Conse- 
quently, the reflected fields from the two surfaces are not equal, and a diffracted 
field is created at the junction which compensates for this discontinuity. The ex- 
pression for this diffracted field is given by Chu[45]. As mentioned before, this 
diffracted field is much smaller than the one resulting from an untreated edge and 
accounts for the increase in performance obtained by using rolled edges. 

3.3.2 Design example (focus feed) 

Consider now the parabolic reflector with f c = 12 ' and elliptic rolled surface 
terminations shown in Figure 19. The values of a e and b e were chosen such that 

°e 

— = f c , and (3.36) 
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Table 3: Surface characteristics for S2.1. 


No Blending 
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mm 


RADIUS OF CURVATURE (ROLLED EDGE) 

&e 

H 


Junction 

Minimum 

Shadow Boundary 

be 

H 


bottom 

top 

bottom 

top 

bottom 

top 

Vbot 

5' 

7 

0° 

0° 

90° 

90° 

89° 

81° 

ytop 

11' 

R c 

12' 

12' 

0.188' 

0.188' 

0.188' 

0.308' 

z pat 

24' 









Table 4: Radii of curvature at the junctions of a parabola with f c = 12'. 


RADIUS OF CURVATURE (PARABOLA) 


bottom(y 6o< = 5') 

top (ytop = 11') 

Rc 

25.579' 

31.947' 


a 


max 

4 


(3.37) 


where A max corresponds to the wavelength at the lowest frequency of operation, 
in this case 1 GHz. The surface is designated as S2.1 with characteristics given 
in Table 3. The radii of curvature at y< op and y^ 0< of a parabola with f c = 12' 
are given in Table 4. Comparison of Tables 3 and 4 shows the discontinuity in the 
radius of curvature of the reflector across the junctions. The radii of curvature 
of the elliptic rolled surfaces can be calculated using the expressions derived in 
Appendix A and are shown in Figures 20 and 21 as a function of y. 

The total field in the target zone was calculated using the corrected PO 


58 





method, with the end point contributions at the incident shadow boundaries of 
the rolled surfaces subtracted. The parabolic section stretches from y bot to y top , as 
does the target zone. The reflector is illuminated by a transverse magnetic (TM) 
line source with a uniform pattern located at the focus. The target zone is consid- 
ered to be located at a distance (zpat) of 24' away from the main origin along the 
2 -axis. Although Section 2.3.3 describes the three-dimensional case of the PO end 
point subtraction, the principle can also be applied to the two-dimensional case. 
Investigation of reflectors with rolled edges has shown that results obtained with 
the two-dimensional corrected PO method compares very favourably with those 
obtained by the moment method[26|. In the case of two-dimensional surfaces, the 
PO integral is evaluated numerically by applying Simpson’s rule. Samples are 
taken at intervals of 0.03A along the surface. For the purpose of this calculation 
it is assumed that there is no direct illumination from the source into the target 
zone. 

Once the total field is known, the GO reflected field can be calculated and 
then subtracted from the total field. The residual field will by default be the fields 
diffracted from the junctions at y bot and yt op , since no other significant scattering 
mechanisms contribute to the total field in the target zone. 

The two-dimensional GO reflected field ( E r ) at a point (y t ,z t ) in the target 
zone, is given by 

e ->fc(a'+«) 

E r {y t ,z t ) = R sh — — (3.38) 

where 

R s,h = Tl (3.39) 

, yf 

3 = + /c , and (3.40) 
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Figure 22: Two-dimensional GO reflected field from a parabolic reflector. 


so that 



+ 9 = fc + Zf . 


(3.41) 

(3.42) 


The fact that s' + s i§ a constant indicates that the reflected field is indeed a plane 
wave, since all the reflected field components in the target plane at zt will have 
the same phase. Note that the sign of the reflection coefficient ( R s / t ) is associated 
with the incident polarization. The geometry associated with Equation (3.38) is 
shown in Figure 22. 

The normalized magnitude and phase of the total field in the target zone, cal- 
culated at 3 GHz, are shown in Figures 23 and 24. Figure 25 shows the magnitude 
of the reflected field and the diffracted fields from the junctions. The corresponding 
magnitudes and phases calculated at 6 GHz and 10 GHz are shown in Figures 26- 
31. 
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Figure 23: Normalized magnitude of the total field from S2.1 (3 GHz). 
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Figure 24: Phase of the total field from S2.1 (3 GHz). 
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Figure 25: Normalized magnitude of the reflected and diffracted fields from 

S2.1 (3 GHz). 


65 



o 



Figure 26: Normalized magnitude of the total field from S2.1 (6 GHz). 
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Figure 27: Phase of the total field from S2.1 (6 GHz). 
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Figure 28: Normalized magnitude of the reflected and diffracted fields from 

S2.1 (6 GHz). 
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Figure 29: Normalized magnitude of the total field from S2.1 (10 GHz). 
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Figure 30: Phase of the total field from S2.1 (10 GHz). 
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Figure 31: Normalized magnit ude of the reflected and diffracted fields from 

S2.1 (10 GHz). 
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In the case of pure elliptic rolled surface terminations, the UTD can also 
be used to calculate the scattered fields in the target zone. However, the UTD 
diffraction coefficients for the higher order blended junctions discussed in the next 
section are not yet available, so that an alternative method, like PO, has to be 
used. 

It is evident from the results that the performance of this reflector is not. up 
to standard. There is a peak-to-peak phase deviation of 14° at 3 GHz (11° at 
10 GHz). The diffracted field is barely 20 dB below the GO level over a 2 ' area of 
the target zone at 3 GHz, and 20 dB below the GO level over a 4 ' area of the target 
zone at 10 GHz. The ripple amplitude level is thus in the order of 0.8 dB. From 
these results it is also evident that the performance of the reflector improves with 
increasing frequency. By examining the low frequency case one can thus perform 
a worst case analysis. 

3.4 Blended rolled surface terminations 

3.4.1 Reduction of diffracted fields using blended rolled edges 

The discontinuity of the radius of curvature across the junction of the elliptic 
rolled surface causes a discontinuity in the reflected fields across the junction which 
gives rise to an unwanted diffracted field. The creation of a more continuous 
reflected field across the junction, and hence a smaller edge diffracted field, can be 
accomplished by smoothly blending the surface of the rolled edge from a parabola 
into an ellipse. The rolled surface is thus transformed from a pure ellipse into a 
blended surface according to the expression [32,46]: 

fblend(l) = fparabolai 7)[1 ~ Hi)] + fellipse(l)Hl) (3.43) 

where 7 is a parametric angle (see Equation (3.27)) and 6 ( 7 ) is a blending function 
which in general varies between zero (at the junction) and one (at the extreme end 
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of the rolled surface); ie., 


6 ( 0 ) = 0 , and 

f>(7m) 1 • 


(3.44) 

(3.45) 


The term f parabola implies the extension of the parabola beyond the junction. The 
extended parabola is shown in Figure 32. Note that the blended surface is attached 
to the parabola at Pj in the same way as the ellipse, so that the slope of the surface 
is smooth across the junction. Equations (3.19)-(3.25) define the local coordinate 
system of the rolled edge. 

It is shown in Appendix A that the use of a blending function as described, 
will ensure that the surface radius of curvature is continuous across the junction; 
ie., the radius of curvature of the blended surface at the junction is equal to that 
of the parabola at the junction in both magnitude and sign. It will be shown later 
in this section that the following conditions will also hold at the junction: 

R[ n ^ is continuous if 6 ^ m ) ( 0 ) = 0 , m = 1 , 2 , . . . ,n - 1 (3.46) 

where Rc ^ is the (n — l)-th derivative of the radius of curvature with respect 
to y, and 6 ^ m ^ is the m-th derivative of the blending function with respect to 7 . 
The selection of a blending function that has 6( m )(0) = 0 , will create junctions 
with increasing degrees of smoothness in terms of the radius of curvature as m 
increases. A smoother junction will, in turn, tend to minimize the diffracted fields 
in the target zone. There is, however, a point of diminishing returns where this is 
no longer true. This point will be discussed again in the next section. A blending 
function with 

6 ( n )( 0 ) / 0 , and 
6 < m )( 0 ) = 0 , m = 1,2 , . . . ,n — 1 
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(3.47) 

(3.48) 



Figure 32: The .r e 2 /e-coordinate system for a blended rolled surface termination. 
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will be referred to as an n-th order blending function. Reference to Equation (3.46) 
shows that the n-th derivative of the radius of curvature of a blended surface 
incorporating an n-th order blending function is discontinuous across the junction, 
while the m-th order derivatives (m = 1, 2 , . . . , n - 1) will all be continuous across 
the junction. 

The following blending functions have been considered: 

1 ) Linear blending: 


Hi) = 

&'( 7 ) = 

i (n, (7) = 

2) Square blending: 

6(7) = 

b '(l) = 
b"( 7 ) = 
(,«( 7 ) = 

3) Cosine blending: 


1_ 

1m 

(3.49) 

— , and 
1m 

(3.50) 

0 for n > 2 . 

(3.51) 


(i) 

27 

(3.52) 

ill 

(3.53) 

2 

—5- , and 
1m 

(3.54) 

0 for n > 3 . 

(3.55) 


Hi) = 

b\j) = 



(3.56) 

(3.57) 
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(3.58) 



(3.59) 

(3.60) 

(3.61) 

(3.62) 


(3.63) 

(3.64) 

(3.65) 

(3.66) 

(3.67) 

(3.68) 

(3.69) 


The case where b( 7) = 1 implies that no blending occurs, so that the blended 
surface remains a pure ellipse. The elliptic curve has already been described as a 
function of 7 in Equations (3.32) and (3.33). However, in order to describe the 
general blended surface, it is necessary to define the extended parabola resulting 
from 6(7) = 0 as a function of 7 as well. Since the blending algorithm has to be 
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applied ultimately to the design of a three-dimensional concave edged reflector as 
well as to the two-dimensional case considered here, the extended parabola will be 
defined as follows: 


Let x m be a maximum distance measured along the positive x e -axis, 
starting at the junction. Let the parametric angle (7) be transformed 
into a linear displacement along this axis, such that 


Ml) = 7— for 0 < 7 < 7m . 
' Ym 


(3.70) 


Using Equations (3.18), (3.26) and (3.70), the extended parabola 
can then be defined by 

Vparabolail ) = 7 ^ x p2 T Uj , and (3.71) 


z parabola{l ) — 


Vparabolah) 


^fc 


(3.72) 


The general two-dimensional blended rolled surface can thus now be expressed in 
the yz-system as 

r blend(lf ) = Dblend^)^ T z 6/end(7)® (3.73) 

where 


yblendilf ) — 


7 1 , 7 ^) Xp2 + Vj 

+ [(a e sin 7 )xp 2 + 6 e (l - cos 7 )y p2 + yj)b( 7) (3.74) 


and 


ZbltndM = 

+ [(fl e sin7)*p 3 +b e {l - cos7)y p3 + z j ] 6(7) . (3.75) 


The direction cosines (x p 2 , x p $, y p %, y p s) are given in Equations (3.22)-(3.25). 
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The radii of curvature of the blended surface and parabola as well as their 
derivatives with respect to y are derived in Appendix A. In order to investigate 
the discontinuities in the derivatives of the radius of curvature across the junction, 
the expressions for the derivatives of the radius of curvature of the blended surface 
and the parabola will now be compared at the top junction (except for the first 
order case where both top and bottom junctions are considered). Consider a first 
order blending function; ie., b'(Q) ^ 0. At 7 = 0, one finds that 

'dR c \ blend (dR c \ parabola 


(dRc\ 

\*y 


dy ) 


where 


±12/c [l + (t/j/2/ c ) 2 ] : 


a e 




+ 


+ b\ 0) C1 


/A^i + faj/2/,) 2 


(3.76) 


(3.77) 


{ x m/ym) { x mMm)^ 2 , 

and yj is the y-coordinate of the junction. The (±) sign corresponds to the top and 
bottom junctions, respectively. In the special case of linear blending, &'(0) = l/7 m , 
so that 

1 2 / 


fr'(0) Cl = 


±12f c [l + ( yj /2f c ) 2 \ 

x m 


dp 


i x m/lfm) 


+ 


fchefi + iVjWc)* 

(®m/7m) 2 


(3.78) 

Consider now a second order blending function; ie., £/(0) = 0 and 6"(0) ^ 0. At 
7 = 0, one finds that 

blend 

/ d“ \ * 

(3.79) 


(d*R c \ bl ' nd (iiR c \r‘ rM <‘ , 

{ iy* ) ~ [WJ + ’’ (0)£2 


where 


€ 2 = 


24 fc [l + {ytop/2fc)“\ 2 / a e 

{ x m/jm)^ \( x m./lim) 


+ 


fcb e \J^ + (y top/2 fcY 

{ x m/lm)2 


(3.80) 
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In the special case of square blending, b"( 0) = 2/7^, so that 

&"(0)e 2 = 48 -k + (yt°p/ 2 fc) ] ^ I a e /Ae yA + (i/fop/2/ c ) 2 A 

X m y ( x m/jm) ( x mhm ) 2 2 J 

( 3 . 81 ) 

In the special case of cosine blending, 6"(0) = tt 2 / 27^, so that 

b"(0)e 2 = U8A4fc t 1 + (yt°p/ 2 fc) 2 ] 2 f a e fcbetfl + (ytop/2f c )2 ~ A 

\(®m/7»n) { x m/ 7m) 2 

( 3 . 82 ) 

Although square blending and cosine blending are both second order blending 
functions, square blending yields a smoother junction in terms of the continuity 
of the second derivative of the radius of curvature at the junction, since the value 
of the second derivative of the square blending function at the junction is smaller 
than that of the cosine blending function for a given set of blending parameters 
(a e ,be,x m , j m ). However, these parameters will not necessarily yield a blended 
surface that adheres to all the design criteria, so that one cannot say that a square 
blend will be better than a cosine blend in general. Consider now a fourth order 
blending function; ie., 6'(0) = 0,6"(0) = 0,6"'(0) = 0 and b iv (0 ) / 0. At 7 = 0, 
one finds that 

(d 4 R c \ blend (d 4 R c \ parahola 

=(-*£?-) +^4 ( 3 . 83 ) 

where 

.. - 6 °/<[ 1 + (ntop/2/.) 2 ] * ( a e ScW + (t/«>p/2/ c ) 2 , \ 

(»m/7m) 4 \(W7m) + 2 j (3 ' 84) 

In the special case of cosine square blending, 6"’(0) = In 4 /2-y 4 n , so that 

b tv ( 0)e 4 = 8767 fc t 1 + ^qp/ 2 /c) 2 ] 2 / ae fcK \Ji + ( ytop/2f c ) 2 A 

x m y (*m/7m) (®m/7m) 2 2 J 

( 3 . 85 ) 
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Based upon the expressions above, one now can make the following general- 


izations: 


1 . 


In the case of an n-th order blending function, the following will hold at the 


junction (7 = 0,y = ytop)' 


U n R c ) 

blend 

(dTRcV 

\ dy n J 

1 = 1 

{ dy 11 ) 

(d m R c \ 

blend 

\ ~ 1 

( d m R c \ 

V dy m ) 

I - ! 

V dy m ) 


+ 6^ n ^(0)e n > and 


parabola 


, m = 1,2 , . . . ,n — 1 


(3.86) 

(3.87) 


where 

an/c[l -Kyj72/c ) 2 ]( n+3 ) / 2 / a e ScW ft + ( yj/2f c ) 2 _ A 

( Xmhm) n \(*m/7m) (®m/7m) 2 2) 

(3.88) 


and a n is a constant. 


2. The error term (e n ) will be zero if 

^ = a e -I- \jai + 2f c b e (yj/2f c )^ . (3.89) 

where 7 m is evaluated in radians. This will lead to continuity in the n-th 
derivative of the radius of curvature across the junction. For the design 
cases considered, however, it was found that blending parameters that sat- 
isfy this equation do not yield blended surfaces that adhere to the design 
requirements. 

Table 5 indicates the orders of various blending functions with their derivatives at 
the junction. 

3.4.2 Design of blended rolled edges 

This section describes a design algorithm for blended rolled edges. As of yet 
there are no closed form design equations available that will yield optimum blended 
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Table 5: Orders of various blending functions and their derivatives at the 

junction. 



First order 

Second order 

Fourth order 


Linear 

Square 

Cosine 

Cosine squared 

fe'(O) 

l/7m 

0 

0 

0 

b"( 0) 

0 

2/7m 

7T 2 / 2 7m 

0 

b"'( 0) 

0 

0 

0 

0 

6* v (0) 

0 

0 

-* 4 / 2 7 m 

3?r 4 /27 m 


edges, so that the design is an iterative process which involves the calculation of 
the scattered fields for every given set of blending parameters (a e , b e , x m ,-y m ). 
The blending parameters are changed from one iteration to the next according to 
whether the diffracted fields become smaller or not. The final design will then be 
based on that set of parameters that yields the smallest edge diffracted field in the 
target zone. This cut-and-try approach is very unsatisfactory since one can never 
be sure that the design is actually an optimum. Furthermore, the calculation of 
the scattered field is a very cumbersome process, especially as a step in an iterative 
process. 

A new design algorithm was developed in the course of this study which is 
based upon an analysis of the continuity of the higher order derivatives of the 
radius of curvature of the surface across the junction. It allows the user to see how 
the surface characteristics that influence the diffracted field change as a function of 
the blending parameters. This enables one to select a set of blending parameters 
that will yield an optimum design. Once the optimum blending parameters have 
been selected, the scattered fields are calculated to verify the design. 
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Figure 33: Axially symmetric reflector with a blended rolled edge. 

This discussion centers around the design of a blended rolled edge for an 
axially symmetric reflector as shown in Figure 33(47,48). The yz-plane can thus 
be considered without loss of generality. It is assumed that the focal length (f c ) 
and the height of the junction between the parabola and the blended edge (ytop) 
are given. The size of the parabolic section is usually equal to or larger than ten 
wavelengths at the lowest frequency of operation ( lOA^x )[49|. The blended edge 
is to be designed such that the minimum radius of curvature in the lit region of the 
rolled edge ( R™ tn ) is equal to a quarter of a wavelength at the lowest frequency of 
operation (Ama.r/4) and the maximum height of the reflector, measured from the 
axis of symmetry of the paraboloid to the top of the rolled edge, is equal to y m ax • 
Note that the skirt shown in Figure 33 is not included in the definition of y m ax- 
The design problem is to find that set of blending parameters (a e ,b e ,x m ,y m ) that 
will satisfy these requirements and results in the smallest diffracted field from the 
junction into the target zone. 

The first step in the design is to select a blending function. It is desirable to 
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select a blending function with as high an order as possible, since this will ensure 
that the higher order derivatives of the radius of curvature are continuous across 
the junction. However, an arbitrarily high order of blending function may make 
it impossible to satisfy the requirements placed on R™ tn and y m ax since a higher 
order blending function requires more space ( y max ~ ytop) in order to obtain a 
smooth junction. Furthermore, there is a point of diminishing returns where the 
continuity of higher derivatives of the radius of curvature across the surface no 
longer yield smaller diffracted fields over the entire target zone. This is due to the 
fact that the transition region spreads out as the junction becomes smoother, as 
shown in Figure 34. Since the diffracted field is strongest in the transition region, 
a higher order blending function consequently results in a weaker diffracted field 
which is prevalent over a larger area of the target zone. In the course of this study 
it was found that a cosine blending function is a good choice when 



ymax Vtop < 5 A max 

(3.90) 

and that 

a cosine squared blending function should be used when 



ymax ~ ytop > 5 \max • 

(3.91) 

If 

ymax ~ ytop ~ 5 A max 

(3.92) 


it is worthwhile to consider both cosine and cosine squared blending functions and 
then to select the one that results in the smallest edge diffracted fields in the target 
zone. Equations (3.90)-(3.92) can be used as a rule-of-thumb when selecting 
a blending function. Since higher order blending functions inherently result in 
smoother junctions, the manufacturing difficulties will also be a consideration. 

In order to facilitate the iterative design process a computer program was 
developed that determines the following values in table format: 
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BLENDED 

ROLLED 



Higher order blending: Weaker diffracted fields 

Lower order blending: Stronger diffracted fields 

Figure 34: The influence of the blending function on the transition region around 

the reflection shadow boundary. 
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1. Radius of curvature and the first four derivatives with respect to y on ei- 
ther side of the junction. Analytic expressions are used to evaluate these 
quantities (see Appendix A). 

2. Radius of curvature and corresponding 7 - value at the incident shadow bound- 
ary. The location of the shadow boundary is found by making use of Equa- 
tion (2.41). 

3. Minimum radius of curvature and corresponding 7 -value along the entire 
blended edge (lit and shadow regions). 

4. Maximum height of the reflector (including parabolic section and blended 
edge, ie. y ma x ) and maximum ^-dimension. 

The values above are determined for a range of x m and parameters when a e and 
fee are fixed. A set of tables are then generated for various combinations of a e and 
fe e . From these tables the sets of blending parameters that satify the constraints 
(R™ tn , y m ax an< l Vtop) can be determined. A typical output using a cosine squared 
blending function is shown in Table 6 . 

There are usually several sets of blending parameters that will satisfy the 
design constraints. However for an n-th order blending function, the first n — 1 
derivatives of the surface radius of curvature are continuous across the junction 
but the n-th derivative is not. The optimum set of blending parameters will be 
the set that satisfies the design constraints and minimizes the discontinuity in the 
n-th derivative of the radius of curvature across the junction. Examination of 
Equation (3.88) indicates that large values of x m , small values of 7 m (ie. a large 
*m/ 7 m ratio), and small values of a e and b e will generally yield better results, since 
they will tend to minimize the error term (e n ). Note that the absolute minimum 
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Table 6: Typical output of the optimization program for the design of blended 

rolled edges. 


COSINE SQUARED BLENDING 

FOCUS=24' 

YTOP=15' 

DESIGN REQUIREMENTS: 

RCMIN = 0.246' (1 GHz) 

YMAX = 20' 


RC RC1 



* = optimum design 

FOCUS = focal length of parabolic reflector 

YTOP = height of junction between parabolic section and blended rolled edge 
RCMIN = minimum radius of curvature along the lit part of the blended rolled edge 
YMAX = maximum height of surface 
XM = maximum linear blending distance (x m ) 

GM = maximum parametric blending angle (7 m ) 

RC = radius of curvature at the junction 

RCn = n-th derivative of the radius of curvature at the junction 
RSH = radius of curvature at the shadow boundary 
GSH = angle at which shadow boundary occurs 
RMIN = minimum radius of curvature along the surface 
GMIN = angle at which minimum radius of curvature occurs 
Lengths in feet, angles in degrees 
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radius of curvature along the entire blended surface usually occurs in the shadow 
region, so that the radius of curvature at the shadow boundary is the minimum 
radius of curvature in the lit region (i2™ m ). However, the absolute minimum radius 
of curvature along the entire blended surface can also occur in the lit region, in 
which case that value will then be i2’ nm . The actual minimum radius of curvature 
in the lit region of the blended surface that is obtained for a given set of blending 
parameters should not be smaller than R™ tn , since that implies a higher, lowest 
frequency of operation. On the other hand, an actual minimum radius of curvature 
in the lit region that is larger than R™ m indicates an over-design and consequently 
not an optimum. An actual maximum height that is smaller than y m ax similarly 
is not an optimum design, since it was found that, in general, smoother junctions 
can be obtained if more height is available to fit the blended rolled edge in. 

The iterative design process thus involves the following steps: 

1. Determine the focal length of the parabola and the height of the junction 
between the parabolic section and the blended edge (yt op ). 

2. Select a maximum total height (y m ax)- 

3. Select the lowest frequency of operation and let R™ tn = A max /4. 

4. Select a blending function of order n. The rules-of-thumb given in Equa- 
tions (3.90)-(3.92) are useful guidelines in this choice. 

5. Generate data tables for a wide range of blending parameters (a e , b e , x m , 7 m ) 
to gain a general idea of the orders of magnitude of the blending parameters 
that will satisfy the design constraints (R ” un , yt op and y m ax)- 

6. Generate data tables using the blending parameter values in all the regions 
around those found in Step 5 until blending parameters are found that satisfy 
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the design constraints. Select that set of blending parameters that has the 
smallest discontinuity in the ra-th derivative of the radius of curvature of the 
surface across the junction This step is repeated until an optimum is found. 

7. Calculate the scattered fields. 

Although it has not been done at this time, the algorithm is suitable for automa- 
tion, so that the complete iteration can be done by a computer. 

As an example, consider the data in Table 6. The design requirements were 
as follows: 


fc 

= 24' 

(3.93) 

ytop 

= 15' , and 

(3.94) 

l /max 

- 20' . 

(3.95) 


The lowest frequency of operation is 1 GHz, which corresponds to a free space 
wavelength of \ ma x = 0.984'. Since the minimum radius of curvature of the 
blended edge in the lit region has to be equal to A max /4, the design requirement is 
jfmin _ 0.246'. Table 6 indicates that an optimum design, using a cosine squared 
blending function, was obtained for the following blending parameters: 


de 

= 2.60' 

(3.96) 

be 

= 2.40' 

(3.97) 

x m 

= 11.62' , and 

(3.98) 

“Yin 

= 116° . 

(3.99) 


The design is an optimum since the blending parameters in Equations (3.96)-(3.99) 
result in a blended edge that satisfies the design constraints in Equations (3.93)- 
(3.95), and a small change in the parameters result in a larger discontinuity in 
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R c across the junction. Note that for these parameters, the minimum radius 
of curvature along the entire blended surface occurs at the shadow boundary. 
Since cosine squared blending was used, the first three derivatives of the radius of 
curvature are continuous across the junction, but the fourth derivative is not. 

The effect of the error term defined in Equation (3.86) is clearly demonstrated 
in Table 6; viz., 

1. For a e ,b e and x m fixed, (R^fend approaches ( R™ ^parabola ag ^ decreases 

2. For a e ,b e and 7m fixed, (R^fend approaches ( R iv ^parabola ag ^ increases . 

3. Increasing and decreasing the values of a e and b e (with x m and 7m fixed) 

results in enlarging and reducing the discontinuity between ( R™) blend an d 
^ R iv yxirabola j unc tion, respectively. 

Note that this algorithm is used to optimize the blended rolled edge in the yz- 
plane, and is applicable to the three-dimensional reflector in Figure 33 because it 
is axially symmetric. In order to expand the method to optimize edges in planes 
other than the yz-plane, the equations for the derivatives of the radius of curvature 
in Appendix A would have to be derived in three dimensions. 

3.4.3 Design example (focus feed) 

Consider now the parabolic reflector described in Section 3.3 (S2.1), but with 
a cosine blended edge and blending parameters x m = 7' and 7m = 120°. Note that 
this blended surface is not optimized. The design illustrates the improvement that 
can be obtained by blending when the dimensions of the ellipse axes (a e and b e ) 
are kept constant. This surface is designated as S2.2 and shown in Figure 35. The 
surface characteristics are given in Table 7. The radii of curvature at the junctions 
Vbot and y top of a parabola with f c = 12' are given in Table 8. The quantities R' c 
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Figure 35: Two-dimensional parabolic reflector with cosine blended rolled edges 

(S2.2). 
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Table 7: Surface characteristics for S2.2. 


Cosine Blending 


fc 

12' 

RADIUS OF CURVATURE (ROLLED EDGE) 

cig 

3' 

■ 

Junction 


Shadow Boundary 

be 

0.75' 

5 

bottom 

top 



bottom 

top 


7' 

7 

0° 

0° 

O 

oo 

O 

OO 

81° 


7m 

120° 

Rc 

25.579' 

31.947' 

0.234' 

0.209' 

0.248' 

0.256' 

Vbot 

5' 

R'c 

0.638 

1.513 





Vtop 

11' 

R'c 

40.655 

68.451 





Zpat 

24' 



and R" denote differentiation with respect to y. Comparison of Tables 7 and 8 
shows that in the case of the cosine blended edge, R c and R' c are continuous across 
the junction, but R " is not. Based on Equation (3.46), this result is to be expected 
since cosine blending is a second order blending function. The radii of curvature 
of the blended rolled surfaces are shown in Figures 36 and 37. Note that the rolled 
surface is concave at the junction but convex at the extreme end, so that there is 
an inflection point somewhere along the surface where the surface changes from 
concave to convex. The radius of curvature of the surface goes to infinity at this 
point as shown in Figures 36 and 37. Consequently, there is no discontinuity in 
the reflected field in the area of this inflection point such that no diffracted field 
results. 

The reflector is illuminated by a TM line source with a uniform pattern located 
at the focus. The target zone is considered to be located at a distance ( z pa t ) of 
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Table 8: Radii of curvature at the junctions of a parabola with f c = 


RADIUS OF CURVATURE (PARABOLA) 


bottom (y bot = 5') 

to P (mop = n') 

Rc 

25.579' 

31.947' 

R'c 

0.638 

1.513 

R'c 

0.133 

_j 

0.161 



Figure 36: Radius of curvature of the bottom blended surface of S2.2 



LO - 
CO _ 



Figure 37: Radius of curvature of the top blended surface of S2.2. 

24' away from the main origin along the z-axis. For the purpose of this analysis 
it is assumed that there is no direct illumination from the source in the target 
zone. The total field in the target zone can be calculated using the corrected PO 
method, with the end point contributions at the shadow boundaries of the rolled 
surfaces subtracted to yield the “true” total field. The GO reflected field can 
then be subtracted from the total field to determine the fields diffracted from the 
junctions. 

The normalized magnitude and phase of the total field in the target zone, 
together with plots comparing the reflected fields with the fields diffracted from 
the junctions are shown in Figures 38-46, calculated at frequencies of 3 GHz, 6 GHz 
and 10 GHz, respectively. Comparison of Figures 38-46 with Figures 23-31 show 
the improvement in performance of the cosine blended rolled surface (S2.2) over 
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the pure elliptic rolled surface (S2.1). Note that the blended rolled edges of S2.2 
are based on an ellipse that has the same dimensions as the rolled edge of S2.1, so 
the only difference between S2.1 and S2.2 is the blending of S2.2. It is thus clear 
from these results that the blended rolled edges have increased the performance of 
S2.2 substantially compared to that of S2.1. At 3 GHz the diffracted fields from 
S2.2 are 20 dB below the GO field over the entire target zone and 30 dB below 
over a 3' area of the target zone. The ripple amplitude of S2.2 is thus 0.27 dB 
over a 3' area of the target zone compared to 0.83 dB over a 2' area for S2.1. The 
peak-to-peak phase variation of S2.2 is 6° over the entire target zone compared to 
14° for S2.1. At 10 GHz the performance improvement of S2.2 over S2.1 is even 
better. The diffracted fields of S2.2 are 26 dB below the GO field over the entire 
target zone, 40 dB below over a 4.5' area of the target zone, 50 dB below over a 3' 
area and 60 dB below over a 2' area. These levels correspond to ripple amplitudes 
of 0.45 dB, 0.09 dB, 0.03 dB and 0.01 dB respectively. 

3.4.4 Design example (subreflector feed) 

Consider now the case of a parabolic reflector with / c = 7.25' and a cosine 
squared blended edge with blending parameters a e = 3.4', b e = 0.75', x m = 6.8' 
and 7 jji = 120°. This surface is designated as S2.3, and is shown in Figure 47. 
Note that the surface has not been completely optimized, since the optimiza- 
tion algorithm was developed after pattern calculations corresponding to a three- 
dimensional version of S2.3 were done. Nevertheless, the example serves to illus- 
trate the improvements that can be obtained with higher order blending functions. 
The surface characteristics are given in Table 9, and the radii of curvature at y top 
anti Vboi of a parabola with f c = 7.25' are given in Table 10. The quantities R' c , i?", 
R'c and R l c v in Tables 9 and 10 denote differentiation with respect to y. Compar- 
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Figure 40: Normalized magnitude of the reflected and diffracted fields from 

S2.2 (3 GHz). 
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;ure 41: Normalized magnitude of the total field from S2.2 (6 GHz). 



Figure 42: Phase of the total field from S2.2 (6 GHz). 
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Figure 43: Normalized magnitude of the reflected and diffracted fields from 

S2.2 (6 GHz). 
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44: Normalized magnitude of the total field from S2.2 (10 GHz). 



Figure 45: Phase of the total field from S2.2 (10 GHz). 
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Figure 46: Normalized magnitude of the reflected and diffracted fields from 

S2.2 (10 GHz). 
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ison of Tables 9 and 10 shows that in the case of the cosine squared blended edge, 
R c , R'c R” and R"' are all continuous across the junctions, but is not. Based 
on Equation (3.46), this result is to be expected, since cosine squared blending is 
a fourth order blending function. 

In this case, the main reflector is illuminated by a Gregorian subreflector 
system, which in turn is illuminated by a magnetic line source with a uniform 
pattern. A discussion of the equations governing the illumination of the main 
reflector by the subreflector is given in Section 4.5 for the three-dimensional case. 
The two-dimensional case discussed here is analogous to the yz-plane of the three- 
dimensional case, with the provision that the fields in the two-dimensional case 
have a 1 /\/R amplitude dependence. 

The target zone is considered to be located at a distance ( Zpa t ) of 20' away from 
the main origin along the z-axis. In calculating the scattered field in the target zone 
it is assumed that there is no direct illumination from the feed or diffraction from 
the subreflector into the target zone and that the reflector is illuminated only by 
the GO field reflected from the subreflector. Fields diffracted from the edges of the 
subreflector are deliberately ignored in order to establish the performance of the 
main reflector when illuminated by an ideal subreflector system. The motivation 
for this assumption will be given in Chapter V. 

The total field in the target zone can be calculated using the corrected PO 
method, with the end point contributions at the shadow boundaries of the rolled 
surfaces subtracted to give the “true” total field. Subtraction of the GO reflected 
field from the total field will give the fields diffracted from the junctions. The 
normalized magnitude and phase of the total field in the target zone together 
with plots comparing the reflected fields with the fields diffracted from the junc- 
tions, calculated at 3 GHz, are shown in Figures 48—56 for the cases where no 
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Table 9: Surface characteristics for S2.3. 


Cosine Squared 


Blending 


fc 

7.25' 


3.4' 

be 

0.75' 

x m 

6.8' 

7 m 


Vbot 

5.5' 

Vtop 

11.5' 

a ar 

5.25' 

bar 


far 

3' 

P 

5.5° 

a 


z pat 

20' 



RADIUS OF CURVATURE (ROLLED EDGE) 


Junction 

Minimum 

Shadow Boundary 


bottom 

top 

bottom 

top 

bottom 

top 

7 

0° 

0° 

mm 

mm 

84° 

75° 

Rc 

R'c 

Rc 

R'c' 

r ? 

17.739' 

1.217 

0.249 

0.0145 

56.069 

30.148' 

3.037 

0.356 

0.0232 

265.649 

0.267' 

0.188' 


0.332' 







Table 10: Radii of curvature at the junctions of a parabola with f c = 7.25'. 


RADIUS OF CURVATURE (PARABOLA) 


bottom (y& ot = 5.5') 

t°p(y<op = H-5') 

Rc 

17.739' 

30.148' 

R'c 

1.217 

3.037 

R'c 

0.249 

0.366 

R'c 

0.0145 

0.0232 

R? 

0.00211 

0.00087 


end point contributions are subtracted, only first order end point contributions 
are subtracted, and finally, first and second order end point contributions are sub- 
tracted. The effect of the erroneous scattering terms that are imbedded in the 
results obtained from the numerical PO integration are clearly visible when the 
magnitude and phase of the total field without the subtraction of the end point 
contributions (see Figures 48 and 49) are compared to the magnitude and phase of 
the total field with the end point contributions subtracted (see Figures 54 and 55). 
The high frequency ripples in Figures 48 and 49 are created by the false scattering 
centers at the incident shadow boundaries on the blended rolled edges. Note that 
the frequency of the ripple is much lower when these terms are subtracted (see 
Figures 54 and 55). This indicates that the resulting diffractions are originating at 
the junctions between the parabola and the blended edges. Since these junctions 
have a smaller spatial separation than the false scattering centers at the incident 
shadow boundaries on the blended edges, the frequency of the ripple is lower. The 
plot of diffracted fields with the end point contributions subtracted in Figure 56 
clearly shows that the diffracted fields from the parabola/rolled edge junctions 
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have maxima on the reflection shadow boundaries. The normalized magnitude 
and phase of the total field in the target zone together with plots comparing the 
reflected fields with the fields diffracted from the junctions, calculated at 1 GHz 
and 10 GHz, are shown in Figures 57-62. Comparison of Figures 54-56 and 60-62 
corresponding to S2.3 with Figures 23-25 and 29-31 corresponding to S2.1, and 
Figures 38-40 and 44-46 corresponding to S2.2 show the improved performance 
achieved by the reflector with a cosine squared blended edge and subreflector feed 
over the other designs. The subreflector feed of S2.3 resulted in a dramatic de- 
crease in the taper of the reflected field. This concept will be explored further in 
Chapter V. Even though the focal lengths of S2.1 and S2.2 are longer than that of 
S2.3, it is clear that the performance of S2.3 is significantly better than that of the 
other two. This is mainly the result of the higher order blending function that has 
been employed in the rolled edges of S2.3. At 3 GHz the diffracted fields from S2.3 
are 30 dB below the level of the GO reflected field over almost the entire target 
zone and 40 dB below over a 3.5' area of the target zone. These levels correspond 
to ripple amplitudes of 0.27 dB and 0.09 dB, respectively. The peak-to-peak phase 
variation is less than 4° over the entire target zone compared to 6° for S2.2 and 14° 
for S2.1. At 10 GHz the diffracted fields of S2.3 are 42 dB below the GO level over 
the entire target zone and 60 dB below over a 4.5' area. These levels correspond 
to ripple amplitudes of 0.09 dB and 0.01 dB. 

From Table 9 it is seen that the minimum radius of curvature of the rolled 
edge in the lit region is 0.188'. Using the A mo , T /4 rule-of-thumb, this corresponds 
to a minimum operating frequency of about 1 GHz. Figures 57-59 indicate that at 
1 GHz the diffracted fields are 20 dB below the GO level over the entire target area 
and 28 dB below over a 2' area. These levels correspond to ripple amplitudes of 
0.83 dB and 0.35 dB respectively. The peak-to-peak phase variation is 9° over the 
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entire target, zone. Even though this reflector thus does not satify the 0.1 dB ripple 
amplitude requirement at 1 GHz, it will be shown in Section 4.6 that optimizing the 
design of the blended rolled edges allows one to achieve the desired performance. 

Based on the numerical results and the experience gained while calculating 
them, one is justified in saying that the use of blended rolled edges enables one 
to design a reflector which meets the ripple amplitude requirement at higher fre- 
quencies without too much difficulty, as shown in Figure 62. It is, however, much 
more difficult to obtain an acceptable performance over a useful area of the target 
zone at lower frequencies. This is due in part to the fact that the transition region 
around the reflection shadow boundary spreads out as the frequency decreases, so 
that a larger part of the target zone will be intercepted by the transition region 
as shown in Figure 63 . Since the diffracted fields are strongest in the transition 
region, this creates a larger ripple in the plane wave. It is for this same reason 
that it is advantageous to have the target zone as close as possible to the main 
reflector. 
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Figure 48: Normalized magnitude of the total field from S2.3 (3 GHz). No end 

point contributions subtracted. 



Figure 49: Phase of the total field from S2.3 (3 GHz). No end point contributions 

subtracted. 
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Figure 50: Normalized magnitude of the reflected and diffracted fields from 
S2.3 (3 GHz). No end point contributions subtracted. 
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Figure 51: Normalized magnitude of the total field from S2.3 (3 GHz). First 
order end point contributions subtracted. 



Figure 52: Phase of the total field from S2.3 (3 GHz). First order end point 

contributions subtracted. 
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Figure 53: Normalized magnitude of the reflected and diffracted fields from 
S2.3 (3 GHz). First order end point contributions subtracted. 
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Figure 54: Normalized magnitude of the total field from S2.3 (3 GHz). First and 
second order end point contributions subtracted. 



Figure 55: Phase of the total field from S2.3 (3 GHz). First and second order end 

point contributions subtracted. 
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Figure 56: Normalized magnitude of the reflected and diffracted fields from 
S2.3 (3 GHz). First and second order end point contributions subtracted. 
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;ure 57: Normalized magnitude of the total field from S2.3 (1 GHz). 



Figure 58: Phase of the total field from S2.3 (1 GHz). 
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Figure 59: Normalized magnitude of the reflected and diffracted fields from 

S2.3 (1 GHz). 
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* 60: Normalized magnitude of the total field from S2.3 (10 GHz). 



Figure 61: Phase of the total field from S2.3 (10 GHz). 
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Figure 62: Normalized magnitude of the reflected and diffracted fields from 

S2.3 (10 GHz). 
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Figure 63: The frequency dependence of the transition region around the 

reflection shadow boundary. 
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CHAPTER IV 


A CONCAVE EDGED REFLECTOR WITH BLENDED ROLLED 

SURFACE TERMINATIONS 


4.1 Introduction 

The blended rolled surface terminations can now be applied towards the design 
of a three-dimensional main reflector. The main reflector surface is a paraboloid 
with focal length (/ c ), and it can be described by 


x 2 + y 2 
4/e 


(4.1) 


where the z-axis is considered to be the axis of symmetry. The reflector is assumed 
to be symmetric around the t/z- plane, and offset in the positive y-direction. The 
blended rolled surface terminations are based on a blend between the paraboloid 
and ellipses with semi-major and semi-minor axes of lengths of a e and b e , re- 
spectively. In the three-dimensional case, one finds that additional performance 
enhancement can be obtained by shaping the the edge contour separating the 
paraboloid and blended rolled surface terminations. 


4.2 Shaping the edge contour 

The creation of the edge contour at the termination of the parabolic sec- 
tion can be visualized by considering the edge contour to be the intersection of a 
paraboloid by a flat plane as shown in Figure 64(a). It is found that the orientation 
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of this plane has a great influence on the shape of the edge contour and the re- 
sulting edge diffracted fields. A plane parallel to the zy-plane will create a convex 
edge contour as shown in Figure 64(a), with the corresponding zy-projection of the 
reflector shape shown in Figure 64(b). The case of an intersecting plane parallel to 
the zz-plane is illustrated in Figure 65(a). Planes parallel to the zz-plane and the 
yz-plane create “straight” edges; ie., the edges are perceived to be straight when 
projected onto the zy-plane, as shown in Figure 65(b). The edge contours are ac- 
tually curved (see Figure 65(a)), but for focal lengths and reflector sizes typically 
associated with compact range main reflectors, the curvature will be slight. This 
edge contour will forthwith be referred to as “rectangular”. Concave edge contours 
are created when the planes parallel to the zz-plane and the yz-plane are tilted 
as shown in Figure 66(a). The concave edge contours of this reflector shape are 
clearly visible in the zy-projection of the edge contour as shown in Figure 66(b). 
The unit vector normal to the edge contour (n e ) at a point Q E (in the yz-plane) 
directed away from the center of curvature but in the plane of the curvature is also 
shown in Figures 64(a), 65(a) and 66(a). It will be shown that the direction of n e 
and the radius of curvature of the edge contour at Q E have a profound effect on 
the magnitude of the edge diffracted fields. 

The zy-projection of the edge contours of conventional compact range reflec- 
tors have been circular (convex) or rectangular in shape. Edge diffraction analysis 
indicates, however, that convex or rectangular edge contours may not be the best 
choice. A concave contour in fact, results in smaller edge diffracted fields in the 
target zone. 

To substantiate the above claim, consider the UTD expression for the electric 
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(a) Creation of a rectangular edge contour. 



(b) XY-projection of a rectangular edge contour. 


Figure 65: Rectangular edge contour. 
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<c 



(4.2) 



n e = unit vector normal to the edge contour 
at Qe directed away from the center of 
curvature, but in the plane of the curvature 
a = distance from Qe to the field point 
E^(s) = diffracted electric field at the field point 

E 1 {Qe) = incident electric field at Qe 
D = dyadic diffraction coefficient 

p e = radius of curvature of the incident wavefront in the edge-fixed 
plane of incidence 

p g = radius of curvature of the edge at Qe 
s' = unit vector in the direction of the incident propagation 
s = unit vector in the direction of the diffracted propagation , and 
00 — angle between s' and the tangent to the edge at Qe • 

The parameters in Equations (4.2)— (4.4) are shown in Figure 67. Let Q E 
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be fixed at a point in space through which the edge will pass. It is then found 
that all parameters in the spread factor, except p g and n e , are also fixed if the 
source and field points remain fixed. These two parameters determine the shape of 
the edge contour in that they control the radius of curvature of the edge contour 
and the orientation of the plane in which the edge contour lies. It is clear from 
Figures 64(a), 65(a) and 66(a) that the different edge contours result in different 
directions of n e and magnitudes of p g . By choosing an appropriate edge contour 
through Qe, one can thus manipulate p g and n e to minimize the spread factor and 
thereby reduce the energy diffracted from Q e into the target zone. The reason for 
this can be found by examining the physical significance of the spread factor as 
given in Equation (4.3). This factor accounts for the conservation of energy along 
a ray tube by enforcing the condition that all the energy entering the ray tube at 
one end be transmitted to the other end. By changing the size of the edge caustic 
(p c ), rays diffracted from adjacent positions on an edge can be made to be more 
divergent. This will alter the shape of the ray tube and decrease the value of the 
spread factor. 

Consider now a convex, straight and concave edge in the zy-plane, all of which 
are illuminated by the same point source located along the z-axis. The diffracted 
rays from the edges are shown in Figures 68-70. It is seen that the horizontally 
diffracted rays from two adjacent points (Qi and Q 2 ) are parallel in the case of the 
convex edge. By reciprocity, when the edge is illuminated by a plane wave from 
the target zone, rays diffracted from the entire edge contour will intersect at every 
point on the z-axis as shown in Figure 68, so that the edge is a caustic. In contrast 
to the convex edge, the straight edge will have a divergent ray pattern as shown in 
Figure 69. The rays diffracted from the concave edge will diverge even more than 
those from the rectangular edge (see Figure 70). The change from the divergent 
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(b) TOP VIEW 

Figure 68: Diffraction from a convex edge. 
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HORIZONTALLY 
DIFFRACTED RAYS 


(b) TOP VIEW 


Figure 69: Diffraction from a straight edge. 
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rays in the case of the concave edge to the highly convergent rays in the case of 
the convex edge, is caused by the change in direction of the unit vector tangent to 
the edge as the edge contour is traced from Qi to Q 2 - It is clear from these figures 
that a more divergent ray pattern can be obtained by using a concave edge contour. 
It is seen from Equation (4.2) that this will result in a reduction of the diffracted 
field magnitude. 

In order to compare the spread factor and diffracted fields for the three types 
of edges discussed, consider a paraboloid with focal length equal to 12', and an 
offset target zone of width 8' in the ^-dimension and height 6' in the y-dimension. 
Let the center of the target zone be 8 1 high, and located at a distance 24* away 
from the origin along the 2 -axis, as shown in Figure 71. Let the shapes of the 
three offset reflectors be that of the rectangular target zone projected onto the 
paraboloid, and enclosed by the convex, rectangular and concave edge contours, 
respectively. In the case of the convex edge, the target zone is modified to be 6' 
wide in the x-dimension due to the symmetry required by the circular shape of 
the edge contour. The xy-projection of the three reflector shapes are shown in 
Figure 72. Note that the xy-projection of the rectangular edge contour coincides 
with the xy-projection of the target zone, and that these reflectors have no rolled 
edge terminations. 

The concave edge contour is formed by considering a rectangle in the xy-plane, 
hereafter referred to as the “defining rectangle”, of dimensions (x m ; n , x max ) and 
( ymin,ymax ) around the xy-projection of the target zone. The xy-projection of the 
target zone is considered to have dimensions (x; e y<,x r i„^) and (yboti Vtop)- The 
center of the defining rectangle coincides with the center of the xy-projection of 
the target zone. Radial lines are drawn from the center of the defining rectangle 
to its border. The concave edge is then formed by the locii of points on the radial 


129 










lines that are a fixed distance (a e ) from the border of the rectangle as shown in 
Figure 73. The relation between the dimensions of the defining rectangle and the 
xy-projection of the target zone are thus given by 


X max 

~ x right a e 

(4.5) 

x min 

= x left ~ a e 

(4.6) 

ymax 

— ytop + Ge > and 

(4.7) 

Vmin 

Vboi a e • 

(4.8) 


Note that the parameters (y m i n ) and ( Umax ) that are used in this context are not 
the same as those in Section 3.4.2. The center of the xy-projection of the target 
zone, and thus the center of the defining rectangle, is located at (0,y atlff ), where 
Vavg is given by 

Vmax + ymin ytop + Vbot , . „ , 

Vavg - = • (4.9) 

Let (x ax , y a x) be a point on the perimeter of the defining rectangle as shown in 
Figure 73. The corresponding point on the concave edge contour of the paraboloid 
[Pj(xj,Vj,Zj)] f is given by 



The entire concave edge contour of the paraboloid can be generated by tracing 
(xax,y<i; t) around the perimeter of the defining rectangle. Note that the actual 


maximum and minimum dimensions of the reflector will not necessarily be equal 
to the values of x m i n , x max , y m in ailf l Vmax ■ The parameters defining the concave 
edge used in this comparison are given in Table 11. 
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Figure 73: Defining rectangle for a concave edge contour. 
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Table 11: Parameters for the concave edge used in the comparison of the edge 

contours. 


x min 

-7' 

X max 

1' 

ymin 

2' 

Umax 

14' 

fc 

12' 

&e 

3' 


The NEC Reflector Antenna Code[50] was used to analyze the three reflectors, 
using the UTD. Numerical values of a modified spread factor, which is given by 


A 9 = 


Pc 


Pc + s 


(4.13) 


were calculated for rays diffracted into the target zone from points common to all 
three edges. Since s is a constant when is fixed, A 9 is a good measure of p c . 
The common diffraction points are shown in Figure 72, and are defined as follows: 


Q A — to P °f the target zone projection, in the yz-plane 
Qb ~ bottom of the target zone projection, in the yz-plane 
Qc = left of the target zone projection, at the height of the centerline 
(concave and rectangular edge contours) , and 
Qd — left of the target zone projection, at the height of centerline 
(convex edge contour) . 


The spread factor magnitudes are given in Tables 12 and 13 for vertical cuts 
through the center of the target zone and in Table 14 for a horizontal pattern cut 
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through the center of the target zone, for the various edge contours. The magnitude 
of the y-polarized fields diffracted from the three reflectors into the target zone are 
shown in Figure 74. The patterns were taken along a vertical cut through the 
center of the target zone at 10 GHz and also calculated with the NEC Code. The 
following observations and conclusions can be made from these results: 

• The caustic effect of the convex edge is clearly demonstrated by the sharply 
increasing values of A s towards the center of the target zone. This is espe- 
cially true in the case of the horizontal pattern cut. Figure 74 shows the 
strong diffracted field from the convex edge in the center of the target zone. 
This is a very undesirable characteristic. 

• The numerical values of the modified spread factor for the “rectangular edge” 
tend to stay close to 1. Since the edge is almost straight one finds that 

p g -> oo (4.14) 

so that it follows from Equation (4.13) that 

A a « 1 . (4.15) 

• In the case of the concave edge, A s is seen to decrease significantly across the 
target zone, with the highest values along the reflection shadow boundaries. 

• Figure 74 clearly shows that the fields diffracted from the concave edged 
reflector into the target zone are smaller than those diffracted from both 
the rectangular and convex edged reflectors. This is due to the fact that 
the concave edge has a smaller spread factor than the other two edges as 
indicated in Tables 12-14. 
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Table 12: Modified spread factor for diffraction from Q&. 


VERTICAL PATTERN CUT 

Diffraction point: $^=(0. 00, 11. 00, 2. 52) 

Pattern position 

Modified spread factor (A s ) 

xpat 

ypat 


concave 

rectangular 

convex 

0.00 

5.00 

24.00 

0.563 

0.983 

0.960 

0.00 

5.50 

24.00 

0.581 

0.986 

1.052 

0.00 

6.00 

24.00 

0.601 

0.988 

1.175 

0.00 

6.50 

24.00 

0.623 

0.990 

1.353 

0.00 

7.00 

24.00 

0.645 

0.992 

1.646 

0.00 

7.50 

24.00 

0.673 

0.994 

2.275 

0.00 

8.00 

24.00 

0.701 

0.996 

7.353* 

0.00 

8.50 

24.00 

0.734 

0.997 

2.540 

0.00 

9.00 

24.00 

0.772 

0.998 

1.748 

0.00 

9.50 

24.00 

0.815 

0.999 

1.417 

0.00 

10.00 

24.00 

0.865 

0.999 

1.224 

0.00 

10.50 

24.00 

0.925 

1.000 

1.094 

0.00 

11.00 

24.00 

0.999+ 

l.ooot 

0.998t 



Dimensions in 

feet 



*Edge diffraction caustic 
t Reflection shadow boundary 


I 

I 
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Table 13: Modified spread factor for diffraction from Qg. 


VERTICAL PATTERN CUT 


Diffraction point: Qg = 

(0.00,5.00,0.52) 

Pattern position 

Modified spread factor (A 5 ) 

xpat 

ypat 

zpat 

concave 



0.00 

5.00 

24.00 

l.ooot 

l.ooot 

0.998 t 

0.00 

5.50 

24.00 

0.926 

1.000 

1.092 

0.00 

6.00 

24.00 

0.867 

0.999 

1.218 

0.00 

6.50 

24.00 

0.818 

0.999 

1.400 

0.00 

7.00 

24.00 

0.776 

0.998 

1.697 

0.00 

7.50 

24.00 

0.740 

0.997 

2.318 

0.00 

8.00 

24.00 

0.709 

0.996 

6.211* 

0.00 

8.50 

24.00 

0.682 

0.995 

2.741 

0.00 

9.00 

24.00 

0.657 

0.993 

1.854 

0.00 

9.50 

24.00 

0.635 

0.991 

1.494 

0.00 

10.00 

24.00 

0.615 

0.989 

1.287 

0.00 

10.50 

24.00 

0.597 

0.987 

1.148 

0.00 

11.00 

24.00 

0.581 

0.985 

1.047 



Dimensions in 

feet 



*Edge diffraction caustic 


t Reflection shadow boundary 


137 







Table 14: Modified spread factor for diffraction from Qq and Qd- 


HORIZONTAL PATTERN CUT 


Diffraction point: Qc- 

=(-4.00, 8.00, 1.67) 


Diffraction point: Qjj= 

=(-3.00,8.00,1.52) 

Pattern position 

Modified spread factor (A a ) 

xpat 

ypat 

zpat 

concave 

rectangular 





( Qc ) 

(Qc) 

( Qd ) 

-4.00 

8.00 

24.00 


l.ooot 


-3.00 

8.00 

24.00 


0.999 

0.998+ 

-2.00 

8.00 

24.00 


0.998 

1.206 

-1.00 

8.00 

24.00 

0.807 

0.996 

1.639 

-0.50 

8.00 

24.00 

0.783 

0.995 

2.165 

0.00 

8.00 

24.00 

0.762 

0.993 


0.50 

8.00 

24.00 

0.736 

0.991 

2.954 

1.00 

8.00 

24.00 

0.718 

0.989 

1.948 

2.00 

8.00 

24.00 

0.688 

0.985 

1.309 

3.00 

8.00 

24.00 

0.659 

0.979 

1.052 

4.00 

8.00 

24.00 

0.631 

0.974 


Dimensions in feet 


*Edge diffraction caustic 
t Reflection shadow boundary 
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• The performance of the concave edged reflector is thus clearly superior to 
the other two, and one can ascertain that the use of such an edge contour 
will, in general, lead to smaller edge diffractions in the target zone. 


4.3 Three-dimensional rolled surface terminations 

The ellipses forming the basis of the blended rolled surface terminations lie in 
flat planes as described in Sections 3.3 and 3.4.1. In particular, the ellipse attached 
at Pj is defined in a local * e j/ e -coordinate system, with Pj as its origin. In order 
for the ellipse attached at Pj{xj,yj, zj) to be aligned with the parabolic surface it 
is required that 

= — junction (4.16) 

where injunction * s the unit vector normal to the paraboloid at the junction {Pj)- 
The expression for injunction derived in Appendix A. Thus, one finds that 


y e(Pj) = t/pi* -f y p2 y + y p3 z 


(4.17) 


where the direction cosines are given by 


Vpl — 
Vp2 = 

Vp3 = 


\/ x 2 j + y] + 4 fc 

Vi 

n/ x j + y ) + 

-2 f c 


, and 


(4.18) 

(4.19) 

(4.20) 


\J X ) + v) + 4 /c 

Let p be a unit vector in the zy-plane at the projection of Pj onto the xy- 
plane such that p is perpendicular to the radial line connecting the projection of 
Pj and the center of the defining rectangle as shown in Figure 73; ie., 

-{Vax ~ Vavg ) A , X ax » . , . . „ , 

y= Pl x + p 2 y. (4.21) 


P = 


, ^=x H — i - ■ — 

\J x< ax + ( Vax ~ Vavg ) 2 \J x%. x + ( Vax ~ Vavg ) 2 


139 




Figure 74: Magnitude of the diffracted fields from the reflectors used in the edge 

comparison (10 GHz). 
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The unit vector x e can then be defined as 



Ye X p 

Xe ~ |y e x p| ~ XplX + Xp2y + XpZZ 

(4.22) 

where the direction 

cosines are given by 


x pl 

~Up3P2 

VP9 

(4.23) 

x p2 

ypwi 

ypg 

(4.24) 

x p3 

y P iP 2 - y V 2 P\ 

— , and 

yp9 

(4.25) 

yp9 

- \J{yp3P2) 2 + (y P 3Pi) 2 + (y p iP 2 - y P 2 Pi) 2 • 

(4.26) 


Consider now a projection in the zy-plane. The surface normal (h Junctlon ) points 
toward the axis of symmetry of the paraboloid and in particular to the point 
(0,0, 2/ c ). In the case where the target zone is symmetric around the axis of 
symmetry of the paraboloid, the center of the defining rectangle lies on the axis 
of symmetry of the paraboloid, so that all the * e y e -planes will be perpendicular 
to the *y-plane as shown in Figure 75(a). In the case of an offset target zone, 
however, only the two x e y e - planes attached to Pj{0, yj ,zj) will be perpendicular 
to the xy-plane as shown in Figure 75(b). This is due to the fact that the radial 
lines of the defining rectangle are parallel to n junction i* 1 the case of the symmetric 
target zone, and in the case of the offset target zone, they are not. The two cases 
are illustrated in Figure 76. The dark line in the interior of the reflectors shown 
in Figure 75 indicate the junction between the parabolic section and the rolled 
surface terminations. Note that the grid lines on the surfaces of the reflectors in 
these figures are not contour lines in that contour lines generally indicate some 
parameter being constant. 

In order to specify the rolled surface attached at. Pj(xj,yj,zj) in the , rye- 
system, it is necessary to use the transformation from the x e y e -system to the 
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(a) Symmetric target zone. 
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(b) Offset target zone. 


Figure 75: The x e t/ e -planes for symmetric an offset target, zones. 
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(a) SYMMETRIC TARGET ZONE 

iY 



(b) OFFSET TARGET ZONE 

Figure 76: Radial lines for symmetric and offset target zones. 
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zyz-system which is given by 

/ . . \ / 


*( 7 ) 

y( 7) 
*(7) 


\ 


* P i y P 1 
x p2 y P 2 
\ ^3 y P 3 ) 


1 *e(7) ^ 


y e (7) 


+ 


/ \ 

x j 

yj 


\ z i } 


(4.27) 


Following the discussion in Section 3.4, it is seen that in the case where 6 ( 7 ) = 1 , 
the pure elliptic rolled surface attached at Pj is given by 

*( 7 ) = (a e sin 7 )a:pi + 6 e (l - cos 7 )y pl + (4.28) 

y(7) = ( a e sin 7 )iEp 2 + M 1 ~ cos l)y P 2 + Vj > and (4.29) 

z( 7 ) = (a e sinqr)x p3 + b e (l - cos 7 ) 2^3 + Zj ■ (4.30) 

Similarly, in the case where b( 7 ) = 0, the resulting extended parabola attached at 

Pj is given by 

x m \ 


Tl ,^r pi+ 


/ \ I 

*(7) 

y(7) = 7 j x p 2 + Vj , and 
-(7) 


x 2 + y 2 


4 fc 


(4.31) 

(4.32) 

(4.33) 


In the case of the general blending function, the blended rolled surface attached 
at Pj(xj,yj,Zj) can thus be described by 


*(7) 


( x m \ 

= rUJ 


*pi + 


[1 - 6 ( 7 )] 


+[(oe sin 7 )x p i + b e ( 1 - cos 7 )y p i + x j]Hl) 


(4.34) 


y( 7) 




[1 - 6 ( 7 )] 

+[(a e sin 7 )xp 2 + b e { 1 - cos 7 )y p2 + yj]b( 7 ) , and 


(4.35) 
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z( 7) = 


[7(gm/7w)gpl + Sj] 2 + h{xmhm)x p 2 + Vj] 


2l 


4/e 

+[(ac sin 7 ) 2^3 + b e ( 1 - cos 7)^3 + Zj] 6 ( 7 ) 


[1 - 6(7)] 


(4.36) 


The method described here by which the blended rolled surfaces are attached 
to the parabolic section is not restricted to concave edge contours, but can be 
applied to other edge contours as well. Figure 77 shows a parabolic reflector with 
a rectangular edge contour and blended rolled edges. A parabolic reflector with a 
convex edge contour and blended rolled edges will be described in Section 4.6. 

4.4 Main reflector design with focus feed 

Two three-dimensional reflectors with f c = 12' and concave edge contours will 
now be considered. They are designated as S3.1 and S3. 2 , and correspond to the 
two-dimensional surfaces S2.1 and S2.2 in that S2.1 and S2.2 are the yz-plane cuts 
of S3.1 and S3. 2, respectively. The surface characteristics of S3.1 and S3. 2 are 
given in Table 15. 

The target zone is considered to be 8 ' wide and 6 ' high, with the center of 
the target zone 8 ' above the axis of symmetry of the paraboloid and situated at 
a distance ( Zpat ) of 24' away from the main origin along the z-axis, as shown in 
Figure 72. In calculating the scattered fields from the reflectors it is assumed that 
there is no direct illumination from the feed into the target zone. 

The reflectors are illuminated by a Huygens source located at the focal point. 
The Huygens source approximation is based on the assumption that it best rep- 
resents the radiation and polarization properties of a horn feed antenna. A more 
complete motivation for the use of a Huygens source will be given in Chapter V. 

Consider the feed antenna to be located in its own xjyfZj - coordinate system, 
with origin located at the focal point of the main reflector so that Xf = — x as 
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Table 15: Surface characteristics for S3.1 and S3. 2. 



S3.1 

S3. 2 

Blending 

None 

Cosine 

fc 

12' 

12' 

a e 

3' 

3' 

be 

0.75' 

0.75' 

x m 


7' 

7 m 


120° 

Vbot 

5' 

5' 

Vavg 

8' 

8' 

Vtop 

11' 

11' 

x left 

-4' 

-4' 

x right 

4' 

4' 

a 

36.87° 

36.87° 

z pat 

24' 

24' 
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Figure 78: Focus feed coordinate system. 


shown in Figure 78. Let the z^-axis be the axis of symmetry of the feed. The 
z^-axis is tilted at an angle (a) relative to the — z-axis. In this case the angle (a) 
is chosen so that the zj - axis points towards the center of the projection of the 
target zone on the main reflector; ie., 


tan a 


/c-y«V 4 /c ' 

The transformation from the ®yz-system to the xjyjzj - system is given by 


(4.37) 


■1 0 


0 cos a sin a 


0 sin a — cos a 


y + —fc sin a 
z / f c cos a 


(4.38) 


It is assumed that the Huygens source is y^-polarized, so that the radiated 


far zone electric field is given (in feed coordinates) by 


E f{R,6f,4>f) = 


~-jkR 
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where 


e ‘(^/»^/) = 

+ fy(0f t 4f)4f(0f,4f) 

(4.40) 

F e{0fAf) = 

^sin (f>f cos 2 F${O f ,<t>f) , and 

(4.41) 

F^fAf) = 

^cos <f> f cos 2 yj F£(0f,<f>f) . 

(4.42) 


In this case, it is assumed that Ffi = 1 and F£ = 1. The spherical unit vectors are 
given by 


/ x \ 


\ *f / 


cos Of cos (f>f cosOfsinfif —sin Of 
— sin 4>f cos (j)f 0 



( » \ 

\ 

x / 


y/ 

) 



l z / / 


(4.43) 


so that 




f cos Of cos <f)f 


cos Of sin <f>f cos a 
— sin Of sin a 


sin 


+f 


cos a cos 




cos Of sin <f>f sin a 
+ sinOf cos a 

sin a cos <f> f 


( - \ 
x 

y 

l * > 

/ 

(4.44) 


Given a point P p (a; p , y p , z p ) on the surface of the main reflector, Equation (4.38) 
can be used to transform it to P p(xf,yf,Zf) in feed coordinates. The unit vector 
from the feed origin to the point on the surface of the main reflector is thus given 
in feed coordinates by 


One further finds that 



(*/n»y/ n? z fn ) • 


(4.45) 


R = \jx 2 f + yj + z'f 
Of = arccos(z^ n ) , and 


(4.46) 

(4.47) 
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+f = 


(4.48) 


a,rcta.n(y f n / x f n ) if x fn ± 0 
90° if x^ = 0,yf n > 0 , and 

( -90° if xf n = 0,y fn < 0 . 

Once 9f and <f>j are determined, the field incident on P p can be calculated in terms 
of x, y and z, by substituting Equation (4.44) into Equation (4.39). 

4.4.1 Design example (elliptic rolled edges) 

Consider now S3.1 with pure elliptic rolled surface terminations, as shown in 
Figures 79 and 80 for two different viewing angles. 

The patterns of the three-dimensional reflectors can be computed using the 
corrected PO method. As mentioned before, this method is used since the UTD 
diffraction coefficients for higher order blended junctions are not yet available. 
The numerical PO integration is performed by dividing the entire reflector surface 
into small triangular patches as described in Section 2.3.2, so that the maximum 
dimension of a patch does not exceed 0.1 A. Each of the trapezoidal patches shown 
in Figures 79 and 80 are divided into two triangular surface patches. Since the 
patches are very small, it can be assumed that each patch is in the far zone of 
the source, and that each pattern point in the target zone is in the far zone of 
a patch. This implies that the incident field on a patch can be considered to be 
a local plane wave, and that the field scattered from a patch can be calculated 
using the far zone approximation of the PO integral. The total field at a pattern 
point is obtained by summing the contributions of all the illuminated patches 
on the reflector surface and then subtracting the end point contributions around 
the incident shadow boundary contour on the blended rolled surface terminations 
as described in Section 2.3.3. Note however that the second term of the end 
point scattering coefficient depends on the radius of curvature of the surface at 


150 


FEET 



-8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 

FEET 


Figure 79: Three-dimensional parabolic reflector with elliptic rolled edges 

( S3. 1)( Front view). 
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Figure 80: Three-dimensional parabolic reflector with elliptic rolled edgi 

(View angle: 9 — 30 °,<j) = 150°). 


(S3.1 ) 
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the incident shadow boundary. Due to the twisted planes in which the rolled 
surface ellipses are located (see Figure 75(b)), this is a very tedious calculation; 
consequently, the radii of curvature at the top and bottom shadow boundaries are 
approximated by their values in the yz-plane. These values are obtained as part 
of the optimization algorithm described in Section 3.4.2. The second order PO 
correction terms have not been applied for the left and right rolled surface shadow 
boundaries. Since the actual numerical computation of a PO pattern for a three- 
dimensional antenna was found to be almost prohibitively expensive in terms of 
computer time, patterns were calculated at 3 GHz only. 

Once the total field is known, the GO reflected field can be calculated and 
subtracted from the total field. The residual field will by default be the field 
diffracted from the concave edge junction between the parabolic and rolled surfaces, 
since no other significant scattering mechanisms contribute to the total field in the 
target zone. 

Keeping in mind that the caustic distances for the reflected field from a 
paraboloid due to illumination from a point source at the focus are both equal 
to infinity, one finds that the GO field (E r ) reflected from a point Q R on the 
parabolic surface to the target zone is given by 

E r (s) = [-E* + 2(n • E*)n]e->** (4.49) 

where 

E* = incident electric field at Q R 

n = unit vector normal to the surface at Q R (See Appendix A) , and 

s = distance from Q R to the field point . 

Patterns were calculated along a cut through the j/z-plane, located at a dis- 
tance ( Zpat ) of 24' away from the main origin along the 2 -axis. The normalized 
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magnitude and phase of the total field together with plots comparing the reflected 
fields with fields diffracted from the concave junctions are shown in Figures 81-89 
for the cases where no end point contributions are subtracted, only first order end 
point contributions are subtracted, and, finally, first and second order end point 
contributions are subtracted. The tickmarks on the plots indicate the sampling 
points in the target zone used in the numerical calculations. By comparing Fig- 
ures 83 and 89 it is seen that, in this case of pure elliptic rolled surfaces, the 
difference in diffracted fields for the cases where no end point contributions are 
subtracted and first as well as second end point contributions are subtracted, is 
small. This implies that the end point contributions are overshadowed by the rel- 
atively large diffractions coming from the edge junction between the paraboloid 
and the rolled edge. Comparison of Figures 87-89 with Figures 23-25 indicate 
the good agreement between the patterns for the three-dimensional (S3.1) and 
two-dimensional (S2.1) reflectors. 

4.4.2 Design example (cosine blended edge) 

Consider now S3. 2 with cosine blended rolled surface terminations as shown 
in Figures 90 and 91 for two different viewing angles. The scattered fields from 
this reflector were calculated along a cut through the yz-plane, located at a dis- 
tance ( Zp a t ) of 24* away from the main origin along the z-axis. The normalized 
magnitude and phase of the total field in the yz-plane of the target zone together 
with plots comparing the reflected fields with fields diffracted from the concave 
junctions are shown in Figures 92-100 for the cases whex-e no end point contribu- 
tions are subtracted, only first order end point contributions are subtracted and, 
finally, first and second order end point contributions are subtracted. Comparison 
of Figures 98-100 with Figures 38-40 iixdicate the good agreement between the 
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Figure 81: Normalized magnitude of the total field from S3. 1(3 GHz). No end 

point contributions subtracted. 
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Figure 82: Phase of the total field from S3.1(3 GHz). No end point contributions 

subtracted. 
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Figure 83: Normalized magnitude of the reflected and diffracted fields from 
S3. 1(3 GHz). No end point contributions subtracted. 
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Figure 84: Normalized magnitude of the total field from S3. 1(3 GHz). First order 

end point contributions subtracted. 
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Figuie 85: Phase of the total field from S3.1 (3 GHz). First, order end point 

contributions subtracted. 
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Figure 86: Normalized magnitude of the reflected and diffracted fields from S3.1 
(3 GHz). First, order end point contributions subtracted. 
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Figure 87: Normalized magnitude of the total field from S3.1 (3 GHz). First and 
second order end point contributions subtracted. 
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Figure 88: Phase of the total field from S3.1 (3 GHz). First and second order end 

point contributions subtracted. 
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Figure 89: Normalized magnitude of the reflected and diffracted fields from S3.1 
(3 GHz). First and second order end point contributions subtracted. 
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Figure 91: Three-dimensional parabolic reflector with cosine blended rolled edges 

(S3. 2) (View angle: 0 = 30 °,<j) = 150°). 
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Figure 92: Normalized magnitude of the total field from S3. 2 (3 GHz). No end 

point contributions subtracted. 

patterns for the three-dimensional (S3. 2) and two-dimensional (S2.2) reflectors. 
Comparison of Figures 94 and 100 indicate that the false contributions from the 
end point terms do have an effect on the pattern in this case. This is also evident 
from the high frequency ripple visible in Figures 92 and 93. Since the diffractions 
from the junction between the paraboloid and blended rolled edge are smaller than 
in the case with the pure elliptic rolled edge, the false scattering from the incident 
shadow boundaries are more significant. The subtraction of the contributions from 
the end point terms reduce the high frequency ripple errors as shown in Figures 98 
and 99. 

4.5 Main reflector design with Gregorian subreflector 

Consider now S3. 3, a parabolic reflector with f c — 7.25', cosine squared 
blended rolled surface terminations and a concave edge contour. S3. 3 corresponds 
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to S2.3 in that it is the t/z-plane cut. of S3. 3. The surface characteristics of S3. 3 are 
given in Table 16, and the surface is shown in Figures 101-103 for three viewing 
angles. The target zone is considered to be 8 1 wide and 6' high, with the center of 
the target zone 8.5' above the axis of symmetry of the paraboloid, and situated 20' 
away from the main origin along the z-axis. The main reflector is illuminated by 
a Gregorian subreflector system. For the purpose of this discussion it is assumed 
that the main reflector is illuminated by the GO field reflected from the subre- 
flector only. Fields diffracted from the subreflector are thus deliberately ignored. 
The purpose of this deletion is to investigate the performance of the main reflector 
when illuminated by an ideal subreflector system. It is further assumed that there 
is no direct illumination from the feed or diffraction from the subreflector into 
the target zone. This is done in order to simulate the dual chamber configuration 
described in Chapter V. 
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Figure 95: Normalized magnitude of the total field from S3. 2 (3 GHz). First 
order end point contributions subtracted. 



Figure 96: Phase the of total field from S3. 2 (3 GHz). First order end point 

contributions subtracted. 
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Figure 97: Normalized magnitude of the reflected and diffracted fields from S3. 2 
(3 GHz). First order end point contributions subtracted. 
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Figure 98: Normalized magnitude of the total field from S3.2 (3 GHz). First and 
second order end point contributions subtracted. 



Figure 99: Phase of the total field from S3. 2 (3 GHz). First and second order end 

point contributions subtracted. 
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Figure 100: Normalized magnitude of the reflected and diffracted fields from S3. 2 
(3 GHz). First and second order end point contributions subtracted. 





Table 16: Surface characteristics for S3. 3. 
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Figure 101: Three-dimensional parabolic reflector with cosine squared blended 

rolled edges (S3.3)(Front view). 
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Figure 102: Three-dimensional parabolic reflector with cosine squared blended 
rolled edges (S3.3)(View angle: 0 = 30°, <f> = 150°). 
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The Gregorian subreflector is a spheroid which is defined in terms of its own 
aj^i/^z^-coordinate system. The origin of this coordinate system is located at the 
center of the spheroid. The subreflector has a semi-major axis of length (a ar ) in 
the direction of its axis of symmetry (z a -axis), and semi-minor axes of length ( bar ) 
in the directions of the £ a -axis and y a -axis, such that a ar > b sr . The focal points 
of the subreflector are located at a distance (f sr ) from the origin along the positive 
z a - axis (f+) and negative z a -axis (f“) respectively, as shown in Figure 104. The 
focal length is given by 

fsr = V“? r - Wr ■ (4.50) 

The main axis of the subreflector is tilted at an angle (/?) with respect to the 
axis of symmetry of the paraboloid, and the subreflector is placed such that 
coincides with the focal point of the parabolic main reflector (/ c ). The primary 
feed is placed at the other focal point (f“). All rays originating from the point 
source at f~ and reflected from the subreflector will pass through f+, so that f+ 
is thus a caustic for all rays reflected from the subreflector. 

The primary feed antenna is assumed to be a Huygens source and is orientated 
such that its axis of symmetry (zf- axis) is tilted at an angle (a) relative to the main 
axis of the subreflector as shown in Figure 104. Chapter V is devoted towards the 
design of a dual chamber compact range configuration, where the assumptions that 
the diffracted fields from the subreflector and spillover from the feed can be ignored 
are justified, and the purpose of the tilt angles is explained. The transformation 
from main coordinates to subreflector coordinates is given by 
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/ 1 n \ 
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1 0 0 
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where 


Ve\ 

Zel 

Pe 

Pec 


r 0 


y\ 

*1 


r a sin P e 
—r 0 cos f3 e 
Pec + P 

arctan(— if z\ ^ 0 

< 

90° if zi = 0 

\/y i + z i 

~ far sin P , and 
fc ~ far cos/? . 


(4.52) 

(4.53) 

(4.54) 

(4.55) 

(4.56) 

(4.57) 

(4.58) 


The geometry associated with this transformation is illustrated in Figure 105. 

In order to calculate the scattered fields from the main reflector, it is necessary 
to determine the incident field at any point P p (x p ,y p ,z p ) on the main reflector. 
The first step is to convert P p to subreflector coordinates [P p (x p3 ,y pa , z ps )} using 
the transformation given in Equation (4.51). The next step is to find the point 
(P r ) on the subreflector that is the reflection point for the ray intercepting P p . 

Let the subreflector be described in its own coordinate system by the following 
parametric equation (see Appendix A): 

r apheroid( u i v ) = (&sr sinu cosv)x 3 + ( b sr sinusinv)y, + ( a<r cosu)z, (4.59) 


so that P r can be expressed as P 7 .(u r ,v r )- Since the spheroid is a body of revo- 
lution around the 2 ,-axis, the parameter (i>) is the spherical angle (</>), while («) 
corresponds to the spherical angle (0). 

Let r be the vector from P p to f^, such that 

r = fjJ — P p = r ms r (4.60) 
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and 


r = rjx, + r 2 y a + r 3 z a 


where 

r l 

r 2 

r 3 

r ms 

It then follows that 
where 


x pa 

r ms 

~Vpa 

r ma 

far ~ Zp S 

= — , and 

r ma 

~ \f x pa "h Vp S + ( far ~ z pa)^ • 

P r = + It , and 

» = |Pr-Ct| 


(4.61) 

(4.62) 

(4.63) 

(4.64) 

(4.65) 


(4.66) 

(4.67) 


so that 



b ar cos v r sin u r 

= lr\ 

(4.68) 


b sr sinv r sin u r 

= Zr 2 , and 

(4.69) 


&ST COS Uy 

II 

> 

+ 

«■*» 

(4.70) 


' 

arctan(r 2 /ri) 

if ri 7 ^ 0 


v r = < 

90° 

if rj = 0 ,r 2 > 0 

(4.71) 


-90° 

if ri = 0 ,r 2 < 0 , and 


l = < 

(bar 1^2) sin v r sin u r if r 2 7 ^ 0 

(bar / r l ) cos v r sin u r if r 2 = 0 . 

(4.72) 


The unit vector (r) and distance (1) are shown in Figure 106. The reflection 
caustic distance (pc) has the same magnitude as /, but has a negative sign. Note 
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Figure 106: Reflection point on the subreflector. 

that tt r = 180° in the case where r\ = r 2 — 0 since u r = 0 is not permitted by 
the dual chamber configuration described in Chapter V. From the above it is thus 
seen that 

fljr COS ttf ~ f st “i” d$r sin Ur 


where 


dsr — \ 


^sr(r 3 /r 2 )sinv r if r 2 / 0 

bsr(r3/ri)cosv r if r 2 = 0 


and 


sinu. r = 


a sr •Jfy + r — fsrds 

a sr + <( 5r 


(4.73) 


(4.74) 


(4.75) 


The positive sign of the square root, is chosen since sin u r > 0. The reflection point, 
on the subreflector is thus given by 


P r = (6 sr sinti,. cos v r )x 3 + (b S r sin u r sint? r )y s + ( a sr cost/. r )z a 


(4.76) 
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or 


P r — Brs^-a "h VraYa + Z rs Z t 


(4.77) 


where v r is given in Equation (4.71) and u r in Equation (4.75). 

Once P r is known, the field incident on P r from the Huygens source at 
can be calculated. To do this, it is necessary to convert P r to feed coordinates by 
using the following transformation: 
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Let the distance (Rf) from the Huygens source to P r be given by 


(4.78) 



(4.79) 


so that 


Of = arccos (zf/Rf) , and 

(4.80) 


arctan(y^/a:^) 

Xf ^0 


4>f = < 

90° 

Xf — 0 and yf > 0 , or Of = 0° 

(4.81) 


-90° 

Xf — 0 and yf< 0 . 



Note that the choice of <f>f = 90° when Of — 0° indicates that the Huygens source 
is y^-polarized. The field incident on P r can now be expressed in feed coordinates 
as 

—jkRf 

E i?(-R/, #/,<£/) = e l (0f,(t>f )— — — (4.82) 

where e* is given in Equation (4.40). Using the transformation given in Equa- 
tion (4.78), e l can be expressed in subreflector coordinates as 
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(4.83) 


e '(0f,4>f ) = [Fg(0 f ,</> f) cos Of cos <f)f — F^Of, (f>f) sin <j>f]x s 
F[Fq{6 f,<^> f){cos 6 j sin 4>f cos a + sinfy sin a) + 
F^Of^f) c°s <f)f cos a]y 5 

F[Fg(0f, <f>f)(cos Of sin <f>f sin a — sin Of cos a) + 

Ff,{0 f,4> f) cos <f>fsma]z s . 

The incident field at. P p on the main reflector is then given by 

E r (P p ) = [-Efl + 2(n • E l R )h\ A e~> ks (4.84) 

where n is the unit vector normal to the subreflector surface at P r (see Ap- 
pendix A), and 

* = I Pp - Pr | (4.85) 


or 


(4.86) 


3 = \Azps - Xrs ) 2 + ( y p3 - yrs ) 2 + (zps - Zrs ) 2 . 

Since s > /, one finds that the spread factor (A) is given by 

1/(1 -s) in three dimensions , and 

e W 2 \///( a _ /) j n two dimensions . 

As the electric field incident at Pp is now known in subreflector coordinates; ie., 


A = 


(4.87) 


® (Pp) - F x x s -f- E y y s -f- E z z 9 


(4.88) 


it can be converted to main coordinates by using the following transformation: 


/ - ) 
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X, 


y a 

= 

\ ia ) 

V 


10 0 x 

0 cos (3 -sin j3 y • (4.89) 

0 sin ft cos (3 J y z ^ 

This transformation will give the incident electric field on the surface of the main 
reflector in main coordinates, so that it is in a suitable form for use in the PO 
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integration. Note that the total path length of a ray from the phase center of 
the Huygens source to the target plane, which is located at a distance (zpat) from 
the main origin along the z-axis, is constant and equal to 2 a sr + f c + Zpat • The 
reflected field from the main reflector into the target zone is thus still a plane wave. 
The scattered fields from S3.3 were calculated along a cut in the yz-plane, located 
at a distance ( z^t ) of 20' away from the main origin along the 2 -axis, at 3 GHz 
using the corrected PO method. The normalized magnitude and phase of the total 
field together with plots comparing the reflected fields with fields diffracted from 
the concave junctions are shown in Figures 107-115 for the cases where no end 
point contributions are subtracted, only first order end point contributions are sub- 
tracted and, finally, first and second order end point contributions are subtracted. 
Comparison of Figures 107-115 with Figures 48—56 indicate the agreement between 
the patterns for the three-dimensional (S3.3) and two-dimensional (S2.3) reflectors. 
The high frequency ripple in Figures 107 and 108 and the discrepancy in diffracted 
fields between Figures 109 and 115 indicate that the contributions from the false 
scattering centers at the incident shadow boundaries on the blended rolled edges 
have a significant effect on the pattern. Since the diffraction from the junction 
between the paraboloid and the rolled edges is small, the contributions from the 
end points terms are much more pronounced. Comparison of Figure 115 with Fig- 
ure 56 indicates that the two-dimensional prediction of the diffracted field is much 
smaller than that obtained by the three-dimensional results. This discrepancy 
most likely results from an inaccurate subtraction of the end point contributions 
in the three-dimensional case. 

Based on the results obtained from the two-dimensional analysis it was claimed 
in Section 3.4 that a reflector with blended rolled edges will have a significantly 
better performance than one with elliptic rolled edges. In this section it has been 
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Figure 107: Normalized magnitude of the total field from S3. 3 (3 GHz). No end 

point contributions subtracted. 

shown that this claim holds true for three-dimensional reflectors as well. The 
performance improvement that can be obtained by choosing higher order blending 
functions as discussed in Sections 3.4.1 and 3.4.2 have now been verified for the 
three-dimensional case. 

4.0 Comparative design example 

A final design example, intended as a 1/3 scale model, will now be considered. 
A reflector with a concave edge contour is compared to one with a convex edge 
contour. Both reflectors have the same optimized cosine squared blended rolled 
edges. The reflectors are designated as S3.4 and S3. 5 and shown in Figures 116- 
117, respectively. Note that the t/z-profile shown in Figure 118 is the same for both 
reflectors, so that the dimensions of their target zones in the yz - plane are equal. 
The 2 -dimension of the target zone of S3.5 is smaller than that of S3.4 in order 
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to accomodate the circular edge contour. The surface characteristics are given in 
Table 17. 

The reflectors are illuminated by a Huygens source at the focus, tilted at an an- 
gle (a) as given in Equation (4.37). Field patterns were calculated at 3 GHz along 
a cut in the yz-plane, located at a distance ( Zpat ) of 8' away from the main origin 
along the 2 -axis. Fifty sample points were taken in the target zone. Figure 119 
shows the normalized magnitude of the total field in the target zone; whereas, 
Figure 120 shows the magnitudes of the normalized reflected and diffracted fields. 
The phases of the total field in the target zone are shown in Figure 121. 

These results clearly indicate that a concave edge contour has superior perfor- 
mance compared to a convex edge contour. The diffraction pattern of S3. 5 shows 
a large lobe in the center of the target zone, which can be attributed to the caustic 
effect of the convex edge contour. The concave edge contour does not exhibit such 
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Figure 109: Normalized magnitude of the reflected and diffracted fields from S3. 3 
(3 GHz). No end point contributions subtracted. 
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Figure 110: Normalized magnitude of the total field from S3.3 (3 GHz). First 
order end point contributions subtracted. 



Figure 111: Phase of the total field from S3. 3 (3 GHz). First order end point 

contributions subtracted. 
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Figure 112: Normalized magnitude of the reflected and diffracted fields from S3. 3 
(3 GHz). First order end point contributions subtracted. 
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Figure 113: Normalized magnitude of the total field from S3.3 (3 GHz). First 
and second order end point contributions subtracted. 



Figure 114: Phase of the total field from S3. 3 (3 GHz). First and second order 

end point contributions subtracted. 


191 


CO . 
O O 
— Of) 
i 


UJ O 

• — ■ o 
li_ ? 


- 






- 

REFL 

DIFF 

ECTED FI 
RflCTED F 

ELD 

' IELD — 

— 

— 

\ 

\ 





/ 

\ 

\ 

/ \ 




/ 

/ 

r 

\ 

/ \ 

1 

i 


A 

\ ■ 
V / 

A 

1 

! 



i 

'/ 


i 

i -i r i — 

-i — r ■■■< t— i 

1' 

— r "i — i — i — 

— 1 — r— 1—1 — 

— n" i — i — 

— i — i — i — i — i 


O ' 

i-5.5 

i 


7.5 8.5 9.5 

HEIGHT T (FEET) 


10.5 


11.5 


Figure 115: Normalized magnitude of the reflected and diffracted fields from S3. 3 
(3 GHz). First and second order end point contributions subtracted. 
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Table 17: Surface characteristics for S3.4 and S3. 5. 



S3.4 

S3. 5 

Blending 

Cosine Squared 

Edge contour 

Concave 

Convex 

fc 

4' 

4' 

&e 

1.25' 

1.25' 

be 

0.53' 

0.53' 

x m 

3.61' 

3.61' 

7m 

110° 

110° 

Vbot 

1.83' 

1.83' 

Vavg 

2.83' 

2.83' 

ytop 

3.83' 

3.83' 

x left 

-1.33' 

-1' 

x right 

1.33' 

T 

a 

38.96° 

38.96° 

z pat 

8' 

8' 


a caustic effect, and an improvement of about 14 dB is obtained in the center of 
the target zone, exactly where it is most needed. 
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Figure 116: Three-dimensional parabolic reflector with a concave edge contour 
and cosine squared blended rolled edges (S3.4)(Front view). 
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Figure 117: Three-dimensional parabolic reflector with a convex edge contour 
and cosine squared blended rolled edges (S3.5)(Front view) 6 
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Figure 119: Normalized magnitude of the total field from S3.4 and S3. 5 (3 GHz). 
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Figure 120: Normalized magnitude of the reflected and diffracted fields from S3. 4 

and S3. 5 (3 GHz). 
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Figure 121: Phase of the total field from S3.4 and S3. 5 (3 GHz). 
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CHAPTER V 


A DUAL CHAMBER COMPACT RANGE CONFIGURATION 


5.1 Introduction 


This chapter describes a dual chamber compact range configuration, where the 
main reflector and target zone are located in the main chamber and a Gregorian 
subreflector together with the associated feed assemblies in the other. The two 
chambers are separated by an absorber fence, with a small coupling aperture to 
transmit signals from one chamber to the other. It will be shown that such a dual 
chamber system incorporating a Gregorian subreflector has several advantages 
over other configurations, especially in terms of the reduction of the amplitude 
taper and cross-polarization errors in the reflected plane wave. Stray scattered 
fields from the subreflector and spillover by the feed are virtually eliminated by 
the absorber fence. System performance is analyzed with and without the absorber 
fence to show that this is actually the case. It will also be shown how the coupling 
aperture should be shaped to minimize diffractions from the fence itself. 

5.2 Compact range configurations 

A compact range usually employs an offset parabolic reflector to convert the 
spherical wave radiated by a point source located at the focus of the reflector to a 
reflected plane wave for illumination of the target zone, as shown in Figure 122(a). 
An offset feed system is used to reduce aperture blockage, but requires that the 


200 


TARGET ZONE 



EOGES 


(a) FOCUS FED SYSTEM 



SUBREFLECTOR 


TARGET ZONE 



(b) GREGORIAN SUBREFLECTOR FED SYSTEM 



TARGET ZONE 



LOW CROSS-SECTION 
PEDESTAL 


(c) CASSEGRAIN SUBREFLECTOR FED SYSTEM 

Figure 122: Various compact range configurations. 




axis of the feed be tilted at an angle (a) relative to the axis of the reflector[34] for 
proper illumination. Even though an offset system has many advantages, it also has 
several disadvantages. First, the spherical wave radiated by the point source has 
a 1 /R spatial attenuation, and since all the points on the reflector surface are not 
an equal distance away from the focus this gives rise to an amplitude taper in the 
reflected plane wave. The cross-polarization errors inherent with offset focus-fed 
configurations is another shortcoming. Koffman has shown that when an axially 
symmetric parabolic reflector is illuminated by a Huygens source (ie., a — 0), the 
cross-polarized fields in the aperture plane are eliminated[51]. A Huygens source 
results from crossed electric and magnetic dipoles of equal strength, so that the 
field radiated by a y-polarized Huygens source can be expressed as 

e -jk R 

E(R, 0, <f>) — (0 sin <f> + <f> cos <j>) cos 2 (0/2) — — — . (5.1) 

it 

The cross-polarization associated with offset parabolic reflector systems has been 
studied by various researchers, among them Chu and Turrin[52], Dijk et a/[53] 
and Gans and Semplak[54]. In the case of an offset reflector illuminated by a 
balanced feed like a Huygens source, it is found that the cross-polarization orig- 
inates from the tilt angle (a) between the feed and the main reflector axes. The 
cross-polarization vanishes when a — 0 and deteriorates as a increases. 

Many of the problems encountered using a focus-fed system can be relieved 
by using an offset subreflector fed system instead, even though a subreflector feed 
system has its own disadvantages. Two types of subreflector configurations are 
potential candidates for compact range applications; viz., a Gregorian system using 
a concave elliptic subreflector as shown in Figure 122(b), and a Cassegrain system 
using a convex hyperbolic subreflector as shown in Figure 122(c). It will be shown 
that a Gregorian subreflector offers substantial advantages over a Cassegrain. 
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Hannan has shown that the equivalent focal length ( f e ) of a subreflector /main 


reflector configuration is given by[55] 


fe = m f c 


(5.2) 


where 


m = 


e + 1 

|e - 1| 


(5.3) 


In these equations, e is the eccentricity of the subreflector, and f c is the real focal 
length of the main reflector. Note that Equation (5.3) applies to the case where 
the axis of the subreflector is aligned with the axis of the main reflector. Since 
m > 1 for both Gregorian and Cassegrain subreflectors, f e > f c for both systems. 
Consequently, one can use a main reflector with a shorter focal length, which 
ultimately reduces the chamber size. In addition, the edge diffracted fields from 
the main reflector into the target zone are also reduced because the target zone 
can be located closer to the main reflector as shown in Section 3.4.4. A larger 
equivalent focal length leads naturally to a reduced amplitude taper and smaller 
cross-polarization errors in the reflected field. The taper is improved because the 
distance from the virtual focus to the reflector surface is more nearly uniform[56j. 
The cross-polarization is related to the f/D ratio (where D is the diameter of the 
main reflector) and improves if this ratio increases. Since the subreflector/ main 
reflector has a larger equivalent focal length than the focal length of the main 
reflector alone, the angle between the feed and main reflector axes is reduced[52], 
resulting in lower cross-polarization errors. 

The reflected field amplitude taper and cross-polarization errors can both be 
reduced further by tilting the subreflector axis by an angle (/?) relative to the main 
reflector axis and tilting the feed axis by an angle (a) relative to the subreflector 
axis, as shown in Figure 122(b). The tilt (/?) rotates the subreflector pattern 
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relative to the main reflector, resulting in a reduced amplitude taper. Rudge and 
Adatia[57], Dragone[58], Mizugutch et a/[59], Tanaka et a/[60] and Shore[61] have 
shown that the cross-polarization errors in dual offset reflector systems can be 
eliminated if, in addition, the axis of the feed is tilted by an angle (a) relative to 
the subreflector axis such that 


a f3 

tan — = m tan — 

2 2 


(5.4) 


where m is given in Equation (5.3). The primary feed antenna is assumed to be 
a balanced linearly polarized feed. The effect of the two tilt angles is to have the 
depolarization arising from the two offset reflectors cancel each other. Note that 
these equations are valid for both Gregorian and Cassegrain systems. In the case 
where the subreflector axis is tilted by an angle (/?), the equivalent focal length of 
the system is given by[62] 


fe — 


|l-e 2 


fc 


(5.5) 


1 + e 2 — 2e cos f3 

Even though subreflector systems have many desirable qualities, they have 
some serious shortcomings. For example, the scattered fields from the subreflector 
that illuminate the main reflector and target zone as well as the spillover from the 
feed are highly undesirable. Adatia[63] has shown that the diffracted fields from 
the subreflector also have a very detrimental effect on the polarization purity. 
Rader[32] analyzed a Cassegrain subreflector for compact range applications and 
found that the diffractions from the subreflector edges result in an unacceptable 
field quality in the target zone. In order to reduce these diffractions he added 
blended rolled edges to the subreflector, but then found that higher order reflected 
fields from the rolled edges also cause too much ripple in the reflected plane wave. 
These problems make a Cassegrain subreflector unattractive for compact range 
applications. 
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A Gregorian subreflector, on the other hand, allows one to use a dual cham- 
ber approach as shown in Figure 123. Note that this antenna configuration is 
the same as the one described in Section 4.5 with main reflector S3. 3. The yz- 
cut of the antennas is also shown in Figure 47. The main chamber contains only 
the main reflector and target zone. The Gregorian subreflector is oversized and 
placed in the second chamber along with the feed assemblies. The two chambers 
are separated by an absorber fence with a small coupling aperture to transmit the 
signals between the chambers. All the desired geometrical optics (GO) reflected 
rays from the subreflector pass through one of its focal points (see Figure 122(b)) 
which is coincident with the focal point of the main reflector. The coupling aper- 
ture between the two chambers is then positioned at this location so that the 
desired signal can propagate unperturbed from one chamber to the other. The 
diffracted fields from the subreflector and spillover from the feed are attenuated 
by the absorber fence before they can illuminate the main reflector and target 
zone. This approach virtually eliminates the subreflector scattered field and feed 
spillover problems, while retaining the benefits of a subreflector feed (ie., reduced 
amplitude taper and cross-polarization errors). Note that this solution cannot be 
applied to a Cassegrain subreflector, since the reflected rays from a Cassegrain 
diverge instead of passing through a caustic point as shown in Figure 122(c). 

In order to compare the performance levels of the different feed configurations, 
consider the following three-dimensional reflector systems: 

1. A Gregorian subreflector with both the subreflector and feed axes tilted 
(f3 = 5.5°, a = 20°) and illuminated by a Huygens source of the type de- 
scribed in Equation (5.1). The subreflector has semi-major and semi-minor 
axes of lengths 5.25* and 4. SOS', respectively. The dimensions of the sub- 
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Figure 123: Dual chamber compact range configuration. 







reflector and tilt angles (/ 3 and a) are chosen to optimize the taper of the 
reflected field[64], while maintaining the relationship between a and (3 given 
in Equation (5.4). The main reflector has a focal length of 7.25' and cosine 
squared blended rolled edges with a concave edge contour to reduce edge 
diffractions (S3. 3 described in Section 4.5). The yz-cut of this configuration 
is shown in Figure 123. 

2. The same configuration as in the previous case, except that the subreflector 
axis is aligned with that of the main reflector, and the feed axis is tilted 
so that the reflected ray along this direction illuminates the center of the 
projection of the target zone onto the parabolic reflector (/ 3 = 0°,a = 17°). 

3. A focus-fed system with a Huygens source at the focus of a parabolic reflector. 
The feed axis is tilted to point at the center of the projection of the target 
zone onto the parabolic section (a = 39° as shown in Figurel22(a)). The 
main reflector has a focal length of 12* . 

Figure 124 shows the normalized GO reflected field along a vertical cut through 
the center of the target zone for the three cases. The amplitude tapers for cases 1, 
2, and 3 are 0.05 dB, 0.3 dB and 1.3 dB, respectively. From this plot it is clear that 
even though the main reflectors in the first two cases have shorter focal lengths 
than that in the third case, the use of the subreflectors result in a substantial 
reduction in the taper of the reflected field. The taper is reduced even more 
when the subreflector axis is tilted relative to the main reflector axis. The cross- 
polarized GO fields relative to the copolarized GO fields are shown in Figure 125 
along a horizontal cut. through the center of the target zone. Note that the cross- 
polarization of the system where the subreflector axis is aligned with that of the 
main reflector (case 2) is much smaller than that of the focus-fed system (case 3), 
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Figure 124 Normalized copolarized reflected fields in the target zone associated with 
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Figure 125: Cross-polarized fields in the target zone associated with various 
compact, range configurations. Note that the cross-polarized fields are zero in the 
case where both subreflector and feed axes are tilted. 


209 



and that the cross-polarization has been eliminated for the case where both the 
subreflector and feed axes are tilted (case 1). The total field along a vertical 
cut through the center of the target zone for case 1 is shown in Figure 113, and 
includes the reflected field and the fields diffracted from the edges of the main 
reflector. Note that this result is uniform to within a few tenths of a decibel which 
makes for a very good approximation to a plane wave. 

The advantages of using a Gregorian subreflector has been clearly demon- 
strated in this section. In the next section it will be shown how an extension 
of the subreflector surface reduces the edge diffracted fields from the subreflector 
that illuminate the main reflector . The rest of the chapter will be devoted to 
showing how the dual chamber configuration can significantly reduce the effects of 
the diffracted fields from the subreflector and spillover from the feed in the target 
zone. Particular attention is paid to the absorber fence which is used to separate 
the two chambers. 

5.3 Geometry of a Gregorian subreflector 

The Gregorian subreflector is a prolate spheroid defined in terms of its own 
x »y$ ^-coordinate system with the origin located at the center of the spheroid, as 
shown in Figure 104. It has a semi-major axis of length (a 3r ) in the direction of the 
axis of symmetry (z,-axis) and semi-minor axes of length ( b sr ) in the directions of 
the x s — and y s - axes, such that a sr > b sr . The subreflector can thus be described 
by the following parametric equation: 

r(xi,v) = b sr c osv sinu x, + b sr sinv sin u y 3 + a sr cos u z s (5.6) 

for t>i < v < V 2 and u\ < u < U 2 . The focal points of the subreflector are located 
at a distance ( f 3r ) from the origin along the positive 2 ,-axis (f+) and negative z 3 - 
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axis (f ar ), respectively. The eccentricity and focal length are given, respectively, 

by 

« = sr /a ar ) 2 , and (5.7) 

far = \] « ar - Hr ■ (5.8) 

The subreflector is placed such that coincides with the focal point of the main 
reflector (/ c ). The primary feed antenna is placed at the other focal point (f^,). 
All rays originating from the point source at f a ~ pass through f+ r after reflection 
from the surface of the subreflector, so that is a caustic point for all rays 
reflected from the subreflector. The fact that f+ is a caustic implies that the 
subreflector reflected field has a virtual point caustic at the focal point of the main 
reflector. Recall that the main reflector converts the spherical field radiated by a 
point source at its focus to a plane wave propagating in the direction of its axis 
of symmetry; ie., into the target zone. The ray reflected from the bottom edge of 
the subreflector will propagate along a reflection shadow boundary (RSB) at 0 r 
as shown in Figure 126. The ray illuminating the top edge will create a RSB in a 
similar way, so that the reflected field will exist only in the region, 0 r < 0 < 180°, 
and will be zero elsewhere. 

Consider now the bottom edge of the subreflector. A diffracted field will 
emanate from this edge in order to make the total field (reflected plus diffracted 
fields) continuous. A Uniform Theory of Diffraction (UTD)[17) analysis shows 
that the diffracted field has a maximum amplitude along the shadow boundary. 
In order to minimize the effect of these diffracted rays on the illumination of the 
main reflector, the subreflector surface should be extended so that the RSB at 6 r 
is directed away from the main reflector. This can be done by increasing un (see 
Equation (5.6)) which will result in a smaller value of $ r . Note that 0 r should not 
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Figure 126: Scattering from the subreflector edges. 
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be too small, since the RSB will then point to the target zone, which is an equally 
undesirable situation. For the sake of future discussion, the fields diffracted from 
the subreflector edges will be referred to as “primary diffracted fields”. 

5.4 Analysis of the absorber fence 

The purpose of the absorber fence is to isolate the two chambers and thus 
prevent diffraction from the subreflector from illuminating the main reflector as 
well as eliminating scattering from the subreflector and spillover from the feed from 
reaching the target zone. The absorber fence consists of a slab of electromagnetic 
absorbing material with a coupling aperture of width (to) around the focal point 
of the main reflector as shown in Figure 127. Since all GO rays reflected from 
the subreflector (ie., desired main reflector illumination) pass through the focal 
point and hence through the aperture, they will not be perturbed by the absorber 
fence. Note however that the concept of a ray implies a high frequency limit. In 
practice, the frequency is not infinitely high, so that the reflected field forms a 
caustic region which is confined to about a wavelength in diameter. The aperture 
width should therefore be large enough to accomodate the reflected beam at the 
lowest frequency of operation. Furthermore, the edges of the coupling aperture 
create incident shadow boundaries (ISB) for the primary diffracted rays, so that 
diffracted rays will emanate from these edges as shown in Figure 127. For the 
sake of future discussion, the fields diffracted from the periphery of the coupling 
aperture will be referred to as “secondary diffracted fields”. 

In order to develop an idea of the magnitude of this problem, the UTD is used 
to calculate a two-dimensional far field pattern at 10 GHz for the field scattered by 
the subreflector shown in Figure 123, in the direction of the main reflector. Two 
cases are considered; viz., 
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Figure 127: Scattering from the absorber fence. 
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1. No absorber fence present. 


2. Absorber fence with a coupling aperture of width (w) present. For the pur- 
pose of this analysis it is assumed that the absorber material is infinitesimally 
thin and that there is no reflection from the material nor any transmission 
through it. 

An electric line source with a uniform pattern is located at f a ~. The subreflector 
surface characteristics are given in Table 16. 

Consider first the case where no absorber fence is present. The yz-coordinates 
of the top and bottom edges of the subreflector are given by (yt op , z< op ) = (O', 2.253') 
and ( Vbot, z bot ) = (-4.576', -2.574'), respectively, where the origin of the coordi- 
nate system has now been shifted to the focal point (f+). Using the UTD solution, 
the scattered far field pattern from the subreflector in the direction of the main 
reflector is given by 

EiiO) = E r + Ef op + E( ot (5.9) 

where E r is the reflected field, and E? op and Efj ot are the diffracted fields from 
the top and bottom edges of the subreflector, respectively, as shown in Figure 128. 
Using Equation (2.3) one finds that the far field pattern of the reflected field is 
given by 

e —j(2ka ar —3ir/4) 

J(2a sr )/l{9) - 1 ' 

In this case the reflection caustic ( p r ) is given by 


(5.10) 


E r {0) 


p r (9) = -l r (9) (5.11) 

where l r is the distance from the point of reflection on the subreflector to the 
focus (f+). In the case of a half plane (n = 2), the diffraction coefficients in 
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Figure 128: Scattering from subreflector with no absorber fence present. 
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Equation (2.11) can be expressed as 

D,,k = !>*(* -*',£*) + R,, h D r (4 + <6', L r ) 

where 

V : 

£> r (<A + ^,I r ) : 


i7 

R-s,h 


(5.12) 



F[jfcl*a(^ - <t>')) 

(5.13) 

2V2irk 

cos j 

s s' 
s + s' 


(5.14) 

_ e -W 4 

F[kL r a(<f> + <f>')\ 

(5.15) 

2y/2-irk 

cos (&$-) 

3 P \ , and 
s + p r 

(5.16) 


acoustic soft/hard reflection coefficients of the 
surface at the point of diffraction . 


The far field pattern of the diffracted field from the top edge is thus given by 
e ~jkl\ 

EfopW = — - 44) - B'Wl + <fi ,pi)]e’ k ^P^ e (5.17) 


where 


Pt = -Hop , and (5.18) 

h = 2 a ar — \p r t \ . (5.19) 


Similarly the far field pattern of the diffracted field from the bottom edge is given 

by 


eLw 


e~H 

T-mh 


A 


44) - D '(fo + (5.20) 
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where 


(5.21) 

(5.22) 


pi = -'Jylot + z bot > and 

4 = — |p£| . 

This is obviously not the complete solution for the total scattered field, but serves 
to show the first order effects of the indicated field components on the pattern in 
the area of interest. 

Consider now the case where an absorber fence is present. The coupling 
aperture in the fence is symmetrical around the focal point and is 2' wide as 
shown in Figure 129, with the rest of the configuration remaining the same. The 
scattered far field pattern in the direction of the main reflector is then given by 

E *M = E r + Et t + E? ight + Ef eft (5.23) 

where E r is the reflected field, Ef oi is the diffracted field from the bottom edge 
of the subreflector, and E? ight and Ef eft are diffracted fields from the right and 
left edges of the coupling aperture, respectively. The far field pattern of the field 
diffracted from the right edge of the aperture of the absorber fence due to illumi- 
nation from the bottom edge of the subreflector is given by 

where 

E right ~~ 

/* - 
( r — 

Similarly the far field pattern of the field diffracted from the left edge of the 
aperture of the absorber fence due to illumination from the bottom edge of the 


= Kigh, - KX)^ 


(5.24) 


— jkl * 

[D l {4>b - 4>bA) ~ + <i>biPb)\ — 7 — , and (5.25) 


(5.26) 


l b 


+ - z bot ) 2 
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Figure 129: Scattering from subreflector with the absorber fence present. 
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subreflector is given by 


E&flW = 4/< (5.27) 


where 



1 } 



m<j>l-4A)~ D r (<f>l+<f>lp r b)\ 



\jylot + ( 


w 


"b z bot 


(5.28) 

(5.29) 


The diffraction coefficient ( D r ) is not included in the expressions for the diffracted 
fields from the edges of the aperture, since it has been assumed that there will be 
no reflections from the absorber and consequently no discontinuity in the reflected 
field. Note that the diffracted field from the top edge of the subreflector is shadowed 
by the absorber wall so that it does not contribute in the main chamber. It will 
be shown in Section 5.6 that shaping the aperture will cause the diffraction from 
the aperture due to direct illumination from the feed to have a minimal effect on 
the system performance. 

Figure 130 shows the normalized total fields given by Equations (5.9) and (5.23); 
whereas, the normalized reflected and total diffracted fields are shown in Fig- 
ure 131. The normalized component fields are shown in Figure 132, with the phases 
of the total and reflected fields shown in Figure 133. It is clear from Figure 130 
that the addition of the absorber fence has reduced the ripple in the reflected field 
substantiably. Examination of Figure 131 shows that in the case where no ab- 
sorber fence is present, the diffracted fields are about 31 dB below the level of the 
reflected fields in the angular region where the main reflector is illuminated. In 
the case where the absorber fence is present, the diffracted fields are about 40 dB 
below the the level of the reflected field in the same angular region. Examination 
of Figure 132 clearly shows the primary reflection shadow boundaries at 9 = 60.6° 
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and 9 — 180°, as well as the fact that the maxima of the diffracted fields occur 
at the shadow boundary and then decay for directions away from the RSB. The 
incident shadow boundaries of the primary diffracted field from the bottom edge 
of the subreflector at 9 — 54° and 9 = 74° are also visible in Figure 132. These 
results thus show that the use of an absorber fence as described would, in principle, 
improve the quality of the field illuminating the main reflector. 

5.5 Shaping of the coupling aperture in the absorber fence to minimize 
diffraction 

Commercially available absorbing material is far from perfect [65], and rays 
incident upon an absorber slab will be able to penetrate through it. The trans- 
mitted rays will, however, be strongly attenuated. The amount of attenuation will 
depend on the material characteristics (which are frequency dependent), but can 
be controlled by the thickness of the slab. Practical considerations dictated by the 
physical space available in the feed chamber will put a constraint on the maximum 
slab thickness that can be used. In addition to the attenuation by the absorber, 
the primary diffracted rays incident upon the absorber slab will also experience 
a decay due to the nature of the diffraction process as they move away from the 
shadow boundary as shown in the previous section. 

In order to reduce the secondary diffraction, it is necessary to shape the ab- 
sorber forming the periphery of the aperture. Keeping in mind that the amplitude 
of the secondary diffracted fields are proportional to the difference between the 
fields on either side of the ISB [66] , a wedge-like shaping with wedge angle (0) 
is most appropriate. Let 9 e be a pattern angle measured from the subreflector 
as shown in Figure 134, with the ISB at 9 e gb . The wedge ensures that the fields 
just outside the absorber ( E\ at 9 e - 9 e sb - e) and those just inside the absorber 
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Figure 130: Normalized total field from the subreflector (10 GHz). 
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Figure 132: Normalized field components from the subreflector (10 GHz). 


224 






Figure 133: Phases of total and reflected fields from the subreflector (10 GHz). 
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Figure 134: Conceptual wedge-like absorber aperture shape. 
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(i?2 at 6 e = @ e sb + c) are the same as e goes to zero, and that the primary diffracted 
rays are attenuated gradually as they move away from the shadow boundary into 
the absorber; ie., as 9 e increases beyond 6 e ab . Since there is then no difference 
between the fields on either side of the shadow boundary, the secondary diffracted 
rays are not first order diffracted fields, but rather slope diffracted fields. Slope 
diffracted fields correct discontinuities in the slope of a field on either side of a 
shadow boundary. If the difference in the slope is small, the slope diffracted fields 
can be much weaker than first order diffracted fields. Although the UTD slope 
diffraction coefficients for a lossy dielectric wedge do not exist at present, the 
principle still applies, as will be shown below. 

In order for the slope diffracted secondary fields to be as small as possible, 
the difference in slope of the field on either side of the shadow boundary should 
also be as small as possible. Let be the distance that a ray travels through 
the absorber, as shown in Figure 134. To keep the slope diffracted fields small it 
is thus necessary to keep dp ot i >s jdO e small, since the attenuation is exponentially 
dependent on p aba , and the phase shift is linearly related to p abs . This concept can 
be translated into the geometrical consideration that the wedge angle (ip) should 
be small, and that the angle with which the diffracted ray intercepts the face of the 
wedge at the ISB (ipi as shown in Figure 134) be close to 90°. On the other hand, ip 
has to be large enough so that the slab reaches a sufficient thickness to attenuate 
the primary diffracted rays propagating in the direction of the main reflector. 
Desirable and undesirable geometries illustrating these concepts are illustrated in 
Figure 134. 

A two-dimensional slab of absorbing material was analyzed using the method 
of moments[67] to investigate the effect of the size of the wedge angle on the mag- 
nitude of the fields diffracted from the apex of the wedge. (Note that the TE/TM 
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convention used in [67] is different from the convention used in this study). The 
method of moments is considered to be an exact solution within the constraints of 
numerical accuracy and in so far as the scattering body can be accurately mod- 
elled. However, because of the large number of variables involved its practical use 
is limited to bodies whose physical dimensions are small in terms of a wavelength. 
In the interest of computational efficiency, finite symmetrical slabs are considered 
as shown in Figure 135. In order to perform the moment method calculation the 
finite slab is divided into n x divisions in the i-direction and n y divisions in the 
^-direction, resulting in n x n y cells. The incident field is a TE-polarized plane 
wave propagating in the z-direction. The patterns are calculated at a relatively 
low frequency in terms of normal compact range use (3 GHz) because of the lim- 
its imposed by the large number of variables that have to used with the moment 
method solution. In an attempt to isolate the diffraction effects from the apexes 
of the wedge, they are spaced a reasonable distance apart. It is assumed that the 
lossy dielectric material can be characterized by a complex permittivity; ie., 

e — c o( f r ~ J e r) (5.30) 

where eo is the permittivity of free space. The permeability is assumed to be that 
of free space (/zo). As an example, the results for slabs with wedge angles of 90°, 
14° and 26.6° are shown in Figures 136-139. The emphasis of these results is on 
comparing the levels of the diffracted fields. The three slabs are designated as Qq, 
Ql and Q 2 respectively, with the geometrical and material parameters given in 
Table 18. Figure 136 shows the normalized magnitudes of the total fields along a 
pattern cut parallel to the y-axis at a distance ( Zpai ) away from the main origin. 
For the cases considered, the scattered fields can be obtained by subtracting the 
incident fields from the total fields. The scattered fields relative to the incident 
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Figure 135: Finite symmetrical absorber slab. 
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Figure 136: Normalized magnitude of total fields from finite symmetrical slabs 
with different wedge angles (3 GHz). 
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Figure 137: Diffracted fields (relative to incident fields) from finite symmetrical 
slabs with different wedge angles (3 GHz). 
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Figure 138: Detail of scattered fields (relative to the incident field) from finite 
symmetrical absorber slabs with different wedge angles (3 GHz). 






Figure 139: Phases of total fields from finite symmetrical slabs with different 

wedge angles (3 GHz). 
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Table 18: Parameters of the finite absorber slabs for comparison of the effect of 

the wedge angle (rp). 



Qo 

Ql 

q 2 



14° 

26.6° 

*1,2/1 


0,-0.1 

0,-0. 3 

*2,2/2 


o 

o 

bi 

0.1, -0.5 

*3,2/3 


0.1, 0.5 

0.1, 0.5 

*4,2/4 


0.1 

0,0.3 

4 

1.5 

1.5 

1.5 

4 

1 

1 

1 

n x 

6 

8 

8 

n y 

51 

60 

60 

x pat 

1 

1 

1 

Dimensions in meters 
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fields are shown in Figure 137. Note that the scattered fields are equal to the 
diffracted fields for y > 0.5 m in these examples. Figure 138 shows the scattered 
fields in detail in the region 0 < y < lm. Note that the peaks of the scattered 
fields do not occur at the geometrical ISB ( y isb = 0.5m), but are shifted towards 
the interior of the slabs. Figure 139 shows the phases of the total fields. Although 
these results are certainly not conclusive proof, they do indicate that, for the 
cases considered, a smaller wedge angle results in smaller diffracted fields. This 
is consistent with the previous contention that a smaller wedge angle results in 
smaller slope diffracted fields. 

5.0 Design of an absorber fence 

The analysis of a subreflector system with a realistic primary feed antenna 
and complex absorber fence is a very complicated problem, primarily due to the 
complexities involved in determining the scattering from a lossy dielectric wedge. 
However, building on the concepts developed above, one can discuss the general 
considerations that should be taken into account when designing an absorber fence, 
and in particular the coupling aperture. 

The first step is to identify critical ray paths; ie., rays on the boundaries of 
desired and undesired illumination regions. The absorber fence should then be 
designed so as to leave the desired rays undisturbed and attenuate the undesired 
rays as much as possible and in such a way that the secondary diffracted rays 
created by the absorber fence have a minimal detrimental effect on the desired 
pattern. The coupling aperture of the absorber fence in Figure 123 is shown in 
detail in Figure 140. The critical ray paths in Figure 140 are indicated by the 
numbers 1-11; viz., 

(1),(2): Reflected from the subreflector to the bottom/top of the parabolic 
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section of the main reflector. These rays should be well clear of any 
part of the absorber fence, since they constitute the desired main re- 
flector illumination. 

(3), (4): Diffracted from the bottom edge of the subreflector to the 
bottom/top of the parabolic section of the main reflector. These rays 
should be strongly attenuated and therefore pass through a thick piece 
of absorber. 

(5), (6): Diffracted from the top edge of the subreflector to the bot- 
tom/top of the parabolic section of the main reflector. These rays 
should be strongly attenuated and therefore pass through a thick piece 
of absorber. 

(7) : Reflected from the subreflector and passes through the top apex of 
the left side of the absorber fence aperture. As this ray moves deeper 
into the absorber, it approaches the primary RSB from the top edge of 
the subreflector. The ISB created by this ray should be directed away 
from the bottom edge junction. 

(8) >(9) : Diffracted from the bottom edge of the subreflector and passes 
through the left/right side of the absorber fence aperture. Since it is 
assumed that the subreflector has an extended surface, these rays will 
propagate in a direction away from the main reflector. 

(10) : Reflected from the bottom edge of the subreflector along the 
RSB. This ray should not be directed towards either the main reflector 
or target zone. 

(11) : Direct from the primary feed antenna into the target zone; ie., 
spillover from the subreflector. This ray should be strongly attenuated 
and pass through a thick piece of absorber. 
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In order to satisfy the requirements of the critical rays, certain geometrical features 
can be incorporated tito the design of the absorber fence. The fence is terminated 
in wedges in the region of the aperture to facilitate a gradual attenuation of rays as 
they pass into the absorber and to ensure that the secondary diffractions from the 
fence are only small slope diffractions, as discussed in the previous section. The 
proposed absorber geometry is not intended to be an optimum design, but serves to 
illustrate the design approach. Figure 140 shows the yz-cut of a three-dimensional 
antenna configuration. The complete design would of course be three-dimensional, 
so that the left and right sides of the fence as shown in Figure 140 are representative 
of any cut through the coupling aperture. Consider now features (A)-(D): 

(A) : This wedge is needed to gradually attenuate rays (3), (4) and (11). 

Ray (8) propagates along the ISB created by the apex of the corre- 
sponding wedge, so that the resulting field diffracted from A will have 
a maximum amplitude in the direction of ray (8). Point A should thus 
be chosen such that ray (8) propagates in a direction away from the 
main reflector and target zone. Note that the ISB at A due to direct 
illumination from the feed, and hence direction of maximum diffrac- 
tion amplitude, is not pointed at the main reflector. The thickness of 
the absorber should be sufficient to eliminate any effects due to these 
diffractions. 

(B) : This wedge is needed to gradually attenuate ray (7), which prop- 
agates along the ISB created by the apex of the corresponding wedge, 
so that the resulting field diffracted from B will have a maximum am- 
plitude in the direction of ray (7). Point B should thus be chosen such 
that ray (7) does not point in the direction of the parabolic section of 
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the main reflector or the junction. Note that A and B together with 
the slopes of the corresponding wedges should also be chosen such that 
the path length of ray (4) through the absorber is long enough so that 
it is sufficiently attenuated by the time it reaches the parabolic section 
of the main reflector. 

(C) : This wedge is needed to start the gradual attenuation of ray (9). 
This wedge should also be thick enough to attenuate ray (11). 

(D) : The length of absorber from D to the top edge of the subreflector 
is needed to attenuate the primary diffracted rays from the top edge of 
the subreflector. 
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Figure 140: Design of the coupling aperture in the absorber fence. 
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CHAPTER VI 


FORWARD SCATTERING FROM A HOMOGENEOUS LOSSY 
DIELECTRIC SLAB WITH WEDGE TERMINATIONS 

6.1 Introduction 

In this chapter the forward scattered field from a two-dimensional homoge- 
neous slab of electromagnetic absorber with a wedge termination is analyzed. This 
structure is related to the absorber fence described in Chapter V. In order to an- 
alyze the absorber fence a technique is required whereby the forward scattering 
from electrically large structures can be obtained, since the moment method tech- 
nique used in Section 5.5 is only practical for structures that are small in terms of 
a wavelength, and the UTD cannot be used since the diffraction coefficient for a 
lossy wedge is not known. 

The problem of scattering from a dielectric wedge, and particularly a lossy 
dielectric wedge, is a complicated one to solve due in part to the difficult nature 
of the boundary conditions. DeWitt[65] has investigated pyramidal and wedge 
absorber structures that are used mainly to cover the walls, floor, ceiling and 
other static structures in an anechoic chamber. He concentrated on the backscatter 
characteristics of the absorber, but reported that many of the published solutions 
for scattering from dielectric wedges are either too limited or too complicated for 
practical use. In this chapter, a technique combining ray tracing and aperture 
integration is proposed; whereby, the forward scattering from a homogeneous lossy 
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Figure 141: Forward scattering from a semi-infinite homogeneous lossy dielectric 

slab with a wedge termination. 

dielectric slab with a wedge termination is determined. Consider the slab with a 
wedge termination and wedge angle (i/>) as shown in Figure 141. It is required to 
find the field in the region * > 0, with a source located in the region (x < x\ ,y < 
2/1 )• Let a planar aperture be defined along the entire t/-axis. The fields in the 
aperture are calculated through ray tracing, from which the equivalent currents 
in the aperture are found. An aperture integration is then performed over the 
aperture (— oo < y < oo) to find the total field at the field point. 

Consider now that part of the aperture where the slab is present (0 < y < 
oo). The aperture field in this region is found by tracing rays through the slab. 
The equivalent current in the aperture is then expressed in terms of the fields 
transmitted through the slab. Once the aperture current has been determined, the 
contribution to the total field at the field point from this part of the aperture is 
found by numerical integration. Since the integration would have to be performed 
numerically from y = 0 to y = oo, this is not very feasible from a numerical 
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viewpoint. However, since the slab is assumed to be very lossy (it is after all an 
electromagnetic absorber), it is reasonable to assume that the field transmitted 
through it will only be significant in that area of the aperture where the wedge 
termination is present; ie., in the area around (0 < y < y\). The aperture field 
in that part of the aperture where the parallel section of the slab is present (yj < 
y < o o) will reduce as y increases, since rays propagating through the entire 
width of the slab will be strongly attenuated. As an approximation, the numerical 
integration need thus only be performed from y = 0 up to some point (y = y sma u) 
where the transmitted field has been attenuated sufficiently to have an insignificant 
contribution. This point ( y small ) corresponds the aperture location where the field 
becomes smaller than a predetermined level (E™ tn ) and depends on the accuracy 
of the solution required, the geometry and material characteristics of the slab and 
the cost of running the integration on a computer. Note that the aperture field 
will be smaller than or equal to E™ m for all y > y amo //. In the case of a finite slab 
with wedge terminations on both ends, the integration can be performed over the 
entire length of the slab if the length of the slab is reasonable, as will be shown in 
Section 6.2. 

Consider now that part of the aperture where the slab is not present ( — oo < 
y < 0). It is shown in Appendix B that the integral in this region can be eval- 
uated by the method of stationary phase. This leads to closed form expressions 
corresponding to a direct ray (stationary term) and a diffracted ray (end point 
term). 

It is assumed that the homogeneous lossy dielectric material can be charac- 
terized by a complex permittivity; ie., 

e = CO (4 - j € r ) (6.1) 
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Table 19: Relative permittivities of typical absorber material. 


Frequency (GHz) 

n 

e" 
c r 

1 

3.00 

0.59 

3 

1.50 

0.69 

10 

1.45 

0.58 

100 

1.30 

0.03 


where e 0 is the permittivity of free space. The permittivity of absorber material 
is highly frequency dependent as shown in Table 19, where e' r and e" values for 
a typical absorber material[65] are listed for various frequencies. The permeabil- 
ity is assumed to be that of free space (yo)- The two-dimensional TE case is 
considered throughout, but the technique can also be applied to the TM and/or 
three-dimensional cases. 

0.2 Finite homogeneous lossy dielectric slab with wedge terminations 
and plane wave illumination 

Consider now a finite homogeneous slab of lossy dielectric material with wedge 
terminations and plane wave illumination as shown in Figure 142. The plane wave 
is assumed to propagate in the x-direction. In order to obtain the transmitted field 
in the aperture, incident rays must be traced as they propagate through the slab. 
Consider first a wedge termination. The wedge is defined in a xy-coordinate system 
as shown in Figure 143, with the wedge corners designated as (xi,yi) and (®2,y2)> 
For the purpose of this discussion (x 2 ,y 2 ) is located at the origin. The aperture 
is considered to be along the entire y-axis. Analysis of the ray path indicates 
that there is a transmission coefficient (Ti) at the free space/dielectric interface 
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Figure 142: Finite homogeneous lossy dielectric slab with wedge terminations. 
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Figure 143: Refraction through a homogeneous lossy dielectric wedge. 
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as well as a propagation constant (/?i) and attenuation constant (ai) as the ray 
propagates through the wedge. These parameters are determined in a x^-system 
that is compatible with the side of the wedge. Once the ray has passed through 
the first interface, it must be transformed into a coordinate system (xy-system) 
that is compatible with the aperture plane, since the two sides of the wedge are 
not parallel. A transmission coefficient (T 2 ) at the dielectric/free space interface is 
then found, after which the equivalent current in the aperture is determined, and 
the aperture integration performed. The GO spread factor (see Equation (2.2)) is 
equal to one in the case of a plane wave incident upon a planar interface. 

0.2.1 Free space/dielectric interface 

Consider now a ray propagating in free space that is obliquely incident at an 
angle (0*) on a lossy dielectric slab as shown in Figure 144. Let the incident field 
be a homogeneous plane wave which is given by 

E* = c e -jk(x' sin 6 l +i/ cos 9 l ) ( 5 .2) 


where k is the free space propagation constant and C is an arbitrary complex 
constant. The incident magnetic field is given by Maxwell’s equations [22] as 

H* - V x E* (6.3) 


or 


H* 


'* 9E ‘ 


;,9E‘ 


w/xo \ dy' dx' ) 
so that the tangential magnetic field at the interface (y 1 = 0) is given by 



kcos ° l Cc -ik(x' sin 0*+y' cos P) 
WHO 


(6.4) 


(6-5) 
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Figure 144: Refraction through a free space/dielectric interface. 


where u> is the angular frequency of the field. Since the reflection angle is known 
to be equal to the incident angle[22], the reflected fields are given by 

E r z = YiCe-j k ( x ' sin0l -y' cos0t ) , and 


HI, = 


:cos 


9 l 


u;/i 0 


r,<7e->*( x ' sin 01 cos 01 ) 


( 6 . 6 ) 

(6.7) 


where Ti is a complex reflection constant. Stratton[68] has shown that the field 
transmitted into the dielectric can be expressed as 


E l z — TiCe~ifo( x ' s * n0t+ d Cos 9 t ) e -ociy l 
Applying Equation (6.4) to Equation (6.8), it follows that 

H 1 , - — ^ 1 Ti C' e - J/31 ( ,r * sin ^ cos 9* ) e -a j ;/ 

* «* 7*0 


(6.8) 


(6.9) 


where 


T\ = complex transmission coefficient 
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Pi = propagation constant in the dielectric 
<xi = attenuation constant in the dielectric , and 
9 t = real transmission angle . 

Examination of Equation (6.8) shows that the transmitted field in the lossy dielec- 
tric is an inhomogeneous plane wave. The phase planes propagate in the direction 

Pl = x' sin 0 l + y' cos 6 t (6.10) 

while the amplitude attenuates in the y'- direction. The transmitted field has to 
satisfy the wave equation[22]; ie., 

V 2 E\ + k\E\ = 0 (6.11) 

where ki is the complex propagation constant in the dielectric, which is given by 

ki - wyfte = ksje' r - je'Ji . (6.12) 

Applying Equation (6.11) to Equation (6.8), one finds that 

a 2 — Pi = —k 2 e' r , and (6.13) 

2 aiPi cos $ i = k 2 e " . (6-14) 

The variables Tj, Tj, 0*, Pi and c*i can be determined by using Equations (6.13) 
and (6.14) and the boundary conditions that the tangential electric and magnetic 
fields should be continuous across the interface at y 1 = 0; ie., 

Ei+El = E*, and (6.15) 

HI, + HI, = Hi, . ( 6 . 16 ) 

Enforcement of the boundary conditions yield the relations 

i + r, = r, (6.17) 
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&COS0*(1 - Ti) = (Pi cos 6^ - jai)T\ , and 
k sin 9 l = Pi sin B t . 


(6.18) 

(6.19) 


Solving for Equations (6.13), (6.14) and (6.17)-(6.19), one finds that 


fll 

ai 

0 1 


Ti 

Ti 



2k cos 6 * 

k cos -f Pi cos 6^ — jc*i ’ &n< ^ 

Ti-l. 


6.2.2 Transformation of coordinates 


( 6 . 20 ) 

( 6 . 21 ) 

( 6 . 22 ) 

(6.23) 

(6.24) 


In order to transform the ray parameters to a coordinate system compatible 
with the aperture (xy-system), consider the transformation shown in Figure 145 
which is given by 


x 1 — x cos 6 — y sin 6 , and (6.25) 

y' = x sin£ + y cos 8 (6.26) 

where 6 is the angle with which the x'y '- system is rotated with respect to the xy- 
system. Note that 6 is related to the wedge angle (ip) by the following equation: 


6 + ip = 90° . 


(6.27) 


Let the electiuc field that is transmitted into the dielectric (as given in Equa- 
tion (6.8)) be expressed as 


E\ = TiCe^' 7 ' x+ ^ y y^ x '^ +y 'y^ 


(6.28) 
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Figure 145: Wedge coordinate transformation. 

where 


7x = 

—j(3\ sin 6 l , and 

(6.29) 

l'y = 

-(«! + jfix cos 6 l ) . 

(6.30) 


The transmitted electric field can, however, also be written as 

j 

El = T\Ce 


(6.31) 


Equation (6.31) represents the field that is transmitted into the dielectric in the 
*J/-system. Substituting Equations (6.25) and (6.26) into Equation (6.31), one 
finds that 


E\ — TiCe~ x ^°‘ x+ ^ x '> e~ y ^ a y + ^ y ^ 


(6.32) 


where 


sin 8 

(6.33) 

cos 8 

(6.34) 
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Figure 146: Refraction through a dielectric/free space interface. 

0x = sin(0* + 8) , and (6.35) 

(3 y = (3\ cos(e l + 8) . (6.36) 

6.2.3 Dielectric/free space interface 

Consider now a ray propagating through a lossy dielectric such that it is 
obliquely incident on a dielectric/free space interface as shown in Figure 146. Let 
the incident field have the following form: 

E l z = Ce~ x ( ax+jf3x ')e~ y ( a v+ jf3 v') (6.37) 

where C is an arbitrary complex constant, and a x , f3 x , a y and (3 y are given in 
Equations (6.33)-(6.36). Applying Equation (6.4) to Equation (6.37), one finds 
that the incident magnetic field is given by 

Hi = + j^ x )e~ x ^ x+ ^r.) e ~y(»y+)P y ) . (6.38) 

The reflected fields are then given by: 

El = T 2 Ce x ^ x+ j> 3x '>e-y^y + ^y) , and (6.39) 
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(6.40) 


Hi = -r 2 — (a* +jf3 x )e x ( a * +jl3x )e- y ( a v +j Py'> 

y iOfJ-Q 

where r 2 is a complex reflection coefficient. The field transmitted from the dielec- 
tric into free space is given by 

El = 'p^ e -jMxcos0 a +ysin6 a ) e -a o (-xsin0 a +ycos0 a ) ^ &n d (6.41) 

Hi = - T- -(p 0 cos 0 a + ja o sin 0 a ) 

y O7/U 0 

e -i/3 0 (* cos +y sin ^° ) e -ao ( - * sin +y cos ) ( 6 .42 ) 

where T 2 is a complex transmission coefficient and 0 a is a real angle. Note that 
the field in free space is attenuated in a direction perpendicular to the direction of 
propagation. This is allowed by Maxwell’s equations, as shown below. It is neces- 
sary to introduce a complex propagation constant of the form a + j/3 for the field 
transmitted into free space in order to match the complex propagation constant 
of the tangential field in the dielectric slab at the interface. The transmitted field 

has to satisfy the wave equation; ie., 

» 

V 2 E\ + k 2 E\ = 0 . (6.43) 

Applying Equation (6.43) to Equation (6.41) one finds that 

<4~P l = -k 2 . (6.44) 

The variables T 2 , T 2 , 0 a , (3 q and ao can be determined by using Equation (6.44) 
and the boundary condition that the tangential fields should be continuous across 
the interface at x = 0. Enforcement of the boundary condition yield the following 
relations: 

i + r 2 - T 2 

a y = olqcosB * 1 
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(6.45) 

(6.46) 


(3 y = (3 q sin 9 a , and 

(<*x+j@x){ 1 — T 2 ) = T 2 (— ao sin0° + j@Q cos# a ) . 
Solving for Equations (6.44)-(6.48) one finds that 


(6.47) 

(6.48) 


^ 

“0 = \!0q ~ k 2 
0 a = arcsin 


Po+ + k 2 + v^TiP^TTF 
2 


U/ 


where 


ji _ 2(a, + j fix) 

(a* - ao sin 0 a ) + j(/3 x + (3q cos 6 a ) 
^2 = T 2 -l 


, and 


00+ = <*y + 0y , and 

A>- = 


(6.49) 

(6.50) 

(6.51) 

(6.52) 

(6.53) 

(6.54) 

(6.55) 


0.2.4 Total aperture field in the region y 2 < y < y 2 

Let the incident ray intercept the wedge at (x p ,y p ) as shown in Figure 143. 
The refracted ray will be transmitted to a point (®2 iVq) in the aperture, where 

y q {x P ) = (*2 - Xp) tan(V> - 9 1 ) + y p . (6.56) 

The transmitted magnetic field in that part of the aperture which is illuminated 

by rays propagating through the wedge is given by 

rrapt \ -TiT 2 ( cos 6 a + jao sin 0 a \ -ait/ 

H y P ( y q ) = — ( J e Mf>abs e <*i Kb, (6.57) 

for y 2 <y q < y we , where 

ywe = y q (x i) 

Pabs = \/(*2 ~ x p ) 2 + (t/ g - t/ p ) 2 , and 
J/a6s — Pabs cos ^ 
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(6.58) 

(6.59) 

(6.60) 



if the plane y 1 = 0 is chosen as the phase reference. Note that Hy P as given in 
Equation (6.57) is applicable to rays entering the wedge part of the slab between 
y\ and t/2 as shown by Figure 142. A similar analysis yields the aperture field 
due to rays entering the wedge part of the slab between j/3 and 1/4. The aperture 
field due to rays entering the middle part of the slab (j/i < y < t/4) is obtained by 
setting B l = 0 when evaluating Equation (6.57). This is thus a special case where 
the plane wave is perpendicularly incident on the slab. One finds that the aperture 
field in the middle of the slab due to rays entering the slab between y\ and 3/4 is 
given by 

H 7 = for yi < y < t/ 4 (6.61) 

where k\ is given in Equation (6.12), 


T 3 


Ta 


2k 

k + k\ 
2 kj 
k\ -f- k 


and 


(6.62) 

(6.63) 


6.2.5 Total forward scattered field 


Consider now a homogeneous finite symmetric dielectric slab with wedge ter- 
minations as shown in Figure 142. Since a moment method solution will be used 
to compare the results obtained from the aperture integration technique, the sym- 
metry is utilized in the interest of computational efficiency. Symmetry is not, 
however, required for the technique to be applicable. A pattern cut is taken at a 
distance ( Xpat ) from the origin along a line parallel to the y-axis. The aperture is 
considered to be located along a line parallel to the y- axis, just outside the slab. 
Note that x\ = x<± — 0 in this case. The total field along the pattern cut is found 
by integrating the surface current in the aperture and is given by [11] 

E’Av) = - 1 r My, )Ho\kp) dy, (6.64) 

J —co 
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Figure 147: Composition of the aperture fields for a finite dielectric slab, 
where 

J z = electric surface current in the aperture 

( 2 ) 

Hq = Hankel function of the second kind and order zero 
y = field point position along the pattern cut , and 
P = \J{ X 2 ~ Xpat) 2 + ( Vq ~ V ) 2 • (6.65) 

To perform the numerical integration from — oo to +oo would not be feasible. 
This dilemma is solved, however, by subtracting and adding the incident field 
from the transmitted field in that part of the aperture where the slab is present. 
The total field in the aperture is then the superposition of the incident field, which 
is present from — oo to +oo, and the scattered field ( E\ — E\), which is present 
only in that part of the aperture where the slab is present (3/2 < y < 173). This 
situation is graphically illustrated in Figure 147. The aperture integration need 
thus only be performed over the length of the slab, since integration of El (which 
now extends over the entire aperture) will just give El (with the phase adjusted 
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accordingly) back at the field point. The total field at the field point will thus be 
the superposition of the incident field and the integration of the equivalent current 
produced by E\ — E\ over that part of the aperture where the slab is present. 
Since the aperture is planar, the equivalent surface current is given by[22] 


Jz{y q ) = 2[H% p (y q ) - H'y{x 2 )] for y 2 <y q < ys (6.66) 


where 

_ -jkx 

H' y (x) = (6.67) 

and Hy p (y q ) is given in Equations (6.57) and (6.61). Depending on the geometry 
and material characteristics, areas in the middle part of the slab may exist where a 
particular point in the aperture is illuminated by rays that are transmitted through 
the middle of the slab as well as through one or both of the wedge sides. This 
situation is shown in Figure 142. 

The total field at a field point along the pattern cut is thus given by 

EUv) = £}(*„») - ^ /VfW - ff*(* 2 )]/f‘ 2) (M iy, (6.68) 

2/2 

where 

E\{x) = e~ jkx . (6.69) 


Note that if kp is large, then the large argument form of the Hankel function can 
be used; ie., 

H^\ x ) — > \ —e~^ x ~ n i^ as i-»oo (6.70) 

V 7TX 

so that 

• /a /IT ry 3 e~i kf3 

El(v) = E , ;(x pat )-e”/ 4 Z 0] J—J n - H^)]— dy q (6.71) 

when kp is large. This equation is valid for all pattern points in this example. 
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Table 20: Parameters of finite absorber slabs. 



Ql 

<?2 

<h 

Qa 


14° 

26.6° 

14° 

26.6° 

*1,2/1 

0,0.1 

0,-0. 3 

0,-0.1 

0,-0. 3 

*2, 2/2 

0.1, -0.5 

0.1, -0.5 

0.1, -0.5 

0.1, -0.5 

*3,2/3 

0.1, 0.5 

0.1, 0.5 

0.1, 0.5 

0.1, 0.5 

*4,2/4 

0,0.1 

0,0.3 

0,0.1 

0,0.3 

e' r 

1.5 

1.5 

3.0 

3.0 

4 

1 

1 

0.59 

0.59 

n x 

8 

8 

6 

6 

n y 

60 

60 

55 

55 

n a 

600 

600 

300 

300 

x pat 

1 

1 

1 

1 

Dimensions in meters 


6.2.6 Results 

In order to verify the validity of the aperture integration technique, the scat- 
tered field obtained from Equation (6.71) is compared to that obtained by a mo- 
ment method solution as described in Section 5.5. Four cases are considered; viz., 
geometries Q\ and Q 2 (as discussed in Section 5.5) as well as geometries Q$ and 
Qa- The geometrical and material parameters are given in Table 20. Note that Q 3 
and Qa are the same size as Q\ and Q 2, respectively, but their material character- 
istics are different. Q 3 and Q 4 are not as lossy as Q\ and Q 2, but they are much 
denser. Simpson’s method is used to implement the integration in Equation (6.71) 
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with n a sample points in the aperture. Patterns are calculated at 3 GHz for Q \ and 
Q 2 , and at 1 GHz for Q 3 and Q 4 . Figures 148-152 compare the results obtained 
by applying the moment method and aperture integration techniques to Q\, Fig- 
ures 153-157 for Q 2 , Figures 158-160 for Q 3 , and Figures 161-163 for Q 4 . Due to 
the symmetry, the results are shown in the region y > 0 only. Figures 148, 153, 158 
and 161 show the magnitude of the total field along the pattern cut as well as the 
aperture field and incident field. Figures 149 and 154 show the detail amplitude of 
the ripple in the total field for Q\ and Q 2 at large distances. The scattered field 
can be obtained by subtracting the incident field from the total field in the case of 
the moment method, and is given by the integration term in Equation (6.71) in the 
case of the aperture technique. For the geometries, material characteristics and 
pattern cut considered, the scattered field will be equal to the diffracted field in 
the region y > y 3 . The scattered fields are shown in Figures 150, 155, 159 and 162. 
Figures 151, 156, 160 and 163 show the phase of the total field along the pattern 
cut as well as the phase of the incident field; whereas, Figures 152 and 157 show 
the ripple in the phase of the total field in detail for Q\ and Q 2 . The discrepancies 
between the results obtained from the moment method and aperture integration 
techniques indicate that the latter is indeed only an approximation. Neverthe- 
less, it does predict the scattered field, and thus also the fields diffracted by the 
wedges, to an acceptable accuracy. The main advantage of this technique is that 
electrically large structures can be analyzed in a computationally efficient manner. 
The inaccuracies in the results obtained by the aperture integration technique are 
probably due to the fact that not all mechanisms have been accounted for. The 
most obvious is the effect of reflected rays in the interior of the slab. These rays 
are multiply reflected in the interior of the wedge before being transmitted into 
free space. This effect would be most noticeable near the apex of a wedge, since 
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Figure 149: Detail of the ripple in the magnitude of the total fields for Q\ 

(3 GHz). 


260 



Figure 150: Scattered fields (relative to incident field) for Q\ (3 GHz). 
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Figure 151: Phases of the total fields for Q\ (3 GHz). 
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Figure 152: Detail of the ripple in the phases of the total fields for Qi (3 GHz). 
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Figure 154: Detail of the ripple in the magnitude of total fields for Q 2 (3 GHz). 
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Figure 155: Scattered fields (relative to incident field) for Q 2 (3 GHz). 
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Figure 156: Phases of the total fields for Q 2 (3 GHz). 
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Figure 159: Scattered fields (relative to incident field) for Q 3 (1 GHz). 
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Figure 160: Phases of the total fields for (1 GHz). 
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Figure 161: Normalized magnitude of the total fields for Q 4 (1 GHz). 
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Figure 162: Scattered fields (relative to incident field) for Q 4 (1 GHz). 
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Figure 163: Phases of the total fields for (1 GHz). 
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the path length of a ray in this region through the lossy absorber will be short. 
Since the attenuation is exponentially dependent on this path length, the attenu- 
ation will also be small. As the ray moves away from apex into the interior of the 
slab, the path length in the absorber and hence also the attenuation will increase 
so that multiply reflected rays will be too weak to have much effect. The effect 
due to this mechanism will be smaller in the case of lossier material. Since the 
refracted fields are discontinuous across the edges (*i,t/i) and (*252/2)) diffracted 
fields will emanate from these points to make the total field continuous. These 
diffracted fields will also contribute to the aperture field. Unfortunately, the form 
of this field is not known which leads to the motivation for the use of the aper- 
ture integration technique. The case of rays that are reflected from the side of 
the wedge back into free space have not been taken into account as they will not 
influence the pattern in the region of interest. 

0.3 Semi-infinite homogeneous lossy dielectric slab with a wedge ter- 
mination and line source illumination 

0.3.1 Aperture regions 

The results of the previous section indicate that the proposed aperture in- 
tegration technique can be used to analyze the forward scattering from a finite 
homogeneous lossy dielectric slab with wedge terminations that is illuminated by 
a plane wave. The technique will now be extended to the case of a semi-infinite 
slab with a wedge termination that is illuminated by a line source. Since the 
source is now assumed to be a line source, the incident fields on the slab have 
cylindrical wavefronts rather than the planar wavefronts as in Section 6 . 2 . The 
plane wave transmission coefficients can still be used, since the ray is assumed 
to exhibit local plane wave characteristics. Other than in the plane wave case, 


275 



however, the refracted rays in this case need to be multiplied by a spread factor 
(see Equation (2.2)) that will account for the conservation of energy in a ray tube 
refracted through the free space/dielectric interface. 

Let the wedge be defined in a a:t/-coordinate system as shown in Figure 164, 
with the source located at (x,,t/ s ), and let (x 2 ,y 2 ) be located at the origin for the 
purpose of this discussion. For the purpose of this analysis it is further assumed 
that the source is located such that y s < y\ and x t < x\. As in the previous 
section, the aperture is taken to be along the entire y-axis. Note that in the case 
of a semi-infinite slab, the aperture can be divided into two regions; viz., that part 
in which the slab is present (y > y 2 ), and that part where it is not (y < y 2 ). The 
total forward scattered field in the region x > 0 can thus be expressed as 

BIM) = - eW42 °/i; to < 6 - 72 > 

or 

E a z (p, <j>) = Ef + Ff (6.73) 

where 


<f> 

Hy P {y q ) 

E? 

Ef 


radial pattern distance 
pattern angle 

tangential magnetic field in the aperture 
contribution to El from the aperture where 
the slab is not present (y < y 2 ) 
contribution to El from the aperture where 
the slab is present (y > y 2 ) 
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FIELD POINT 



LINE SOURCE 


Figure 164: Refraction through a lossy dielectric wedge with line source 

illumination. 
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and 


p' = y(*2 ~ x ) 2 + (Vq ~ y) 2 ■ (6.74) 

Note that (x 2 ,y q ) is the point in the aperture and ( x,y ) is the field point; ie., 

x — p sin<^> , and (6.75) 

y — —p cos <f> . (6.76) 

The fields can thus be expressed as 

E? = -e W4 2o^/_“ K?(y,)£j=- iy, , (6.77) 

E f = (6.78) 

In the high frequency case (where k is large), the integral in Equation (6.77) can 

be evaluated by the method of stationary phase to yield 

E s z a ss E 3 z ai + E s z ad (6.79) 


where E z at is the direct incident ray resulting from the stationary point term; ie., 


o < <f> < 4 >q 

E sai _ y/pi -Y-yv 

0 otherwise 

(6.80) 

and E z ad is the ray diffracted from (z2>2/2) resulting from the end point contribu- 

tion, such that 

, e ~jkpo e -jkp 

Et ai =- — ^(*,0,1)— — 
VPo VP 

(6.81) 

where 


Ps = \/(* - x s ) 2 + (y - y,) 2 

(6.82) 

^0 = 7 ) 2 + {y-2 - Vs ) 2 

(6.83) 
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F = transition function (see Equation (2.21)) . 


Equations (6.80)-(6.87) are derived in Appendix B. The closed form expressions 
obtained from the stationary phase approximation method thus eliminates the 
need to evaluate the integral in Equation (6.77) numerically from -oo to y 2 . 

In order to evaluate E s z d , the integral in Equation (6.78) must, however, be 
performed numerically. Note that the integration need not be performed from y 2 
to oo, but as an approximation it need only be performed from y = y 2 to that 
point (y sma u) in the aperture where the aperture field falls below a specified level 
(E™ n ) as discussed in Section 6.1. The tangential magnetic field ( Hy P ) must be 
found at every point in the aperture by tracing rays through the dielectric slab as 
in the previous section. Let the point where a ray enters the slab be designated as 
( x p>2/p)> and the corresponding point where the same ray intercepts the aperture 
be designated as ( x 2 ,yq ). Two cases need to be considered; ie., rays entering the 
wedge (t/2 < Vp < 2 / 1 ) and rays entering the parallel section (y\ < y p ). 

0.3.2 Oblique incidence on a flat slab illuminated by a line source 


Consider now the case where the rays enter the parallel section of the slab at. 
an oblique angle ( 9 jc) as shown in Figure 165. The subscript (/) refers the flat part 
of the slab. In contrast to the case of plane wave incidence where 9^ is the same for 
all rays, the case of line source illumination discussed here results in a different 0 j. 
for each ray. These rays will illuminate the aperture in the region ( y fl < y q < oo), 
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INCIDENT RAY 


Figure 165: Refraction through a lossy dielectric slab with plane wave 

illumination. 


where the ray intercepting the aperture at (* 2 ,Vfl) enters the parallel section at 
(®i, 3 /l) as shown in Figure 164, such that 


Vfl = (*2 ~ X 1 ) tan + y\ 


( 6 . 88 ) 


where is the transmitted angle relative to the normal (x) at 1 , 2 / 1 ) (see Equa- 
tion (6.22)). The equations governing the free space/dielectric interface where the 
ray enters the slab have been derived in the previous section, and are given in 
Equations (6.20)-(6.24). Since the two faces of the slab are parallel, the coordi- 
nate transformation of Equations (6.25) and (6.26) does not apply here, and the 
transmitted field in the lossy dielectric slab thus attenuates in the ^-direction only. 
In order to derive the parameters applicable at the dielectric/free space interface 
where the ray intercepts the aperture, consider a field incident on the dielectric/free 
space interface from the dielectric side, such that 


gi __ £, e -x(<xi+jPl c °sO t j)^-y(j(3 1 sin 6^) 


(6.89) 


I 

I 

I 

I 
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where C is an arbitrary complex constant and ffj is a real angle. Applying Equa- 
tion (6.4) to Equation (6.89), one finds that the tangential incident magnetic field 
is given by 


fji _ g (/3l cos 9 l f J ' Q l) c -g(<>l+J/?i cosgp^-yQV?! sinfl^) 
y up o 


(6.90) 


The reflected fields are given by 


E; = r 5 Ce l(,,1+ift “ !9 / ) e“ ! ' (j *"" # / ) ,and (6.91) 

„ = F C Vl c ° 3 */ ~ J Q 1 ) e «(«n-jft ■!■»«}) (6 92) 

where Fs is a complex reflection coefficient. The fields transmitted into free space 
can be expressed as 


E\ = T 5 Ce-j k ( xcos6a+ y 8in0a ) , and (6.93) 

= _r 5 C'^^e-^(* co ^ 0 +ysine“) (6. 94 ) 

where Ts is a complex transmission coefficient and 6 a is a real angle. Note that 
in the case of the parallel slab, the field that is transmitted into free space does 
not attenuate, in contrast to the field that is transmitted into free space through a 
wedge. (Recall that such a ray attenuates in a direction perpendicular to the direc- 
tion of propagation). The variables Ts, T 5 and 0 a can be found by enforcing the 
boundary conditions that the tangential fields be continuous across the interface 
at x = 0 ; ie., 


0 a = 9\ 


n 

(6.95) 

k sin 6 a , and 

(6.96) 

T$k cos 0 a . 

(6.97) 

, one finds that 



(6.98) 
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so that the ray transmitted through the slab into free space propagates in the same 
direction as the incident ray. Solving for Equations (6.95) and (6.97) , one finds 
that 


r 5 

r 5 


2(/3i cos 0 l f - ja i) 
k COS 6 l f + /3\ COS — ja\ 
— T5 - 1 . 


and 


(6.99) 

( 6 . 100 ) 


The total tangential magnetic field in the aperture due to rays entering the parallel 
slab is thus given by 


H a /\y q ) = 


-Yq T\T§ cos d'f A e a ' y ab» 

0 


for yfi < y q 
elsewhere 

( 6 . 101 ) 


where 


Pair = y/ixi-x^+iVp-V.) 2 (6-102) 

Pab B = \/(*l - x 2) 2 + (yp-yq) 2 » and (6.103) 

Vabs = X 2 ~ Xi . (6.104) 


In the case of a cylindrical wave incident upon a planar interface the spread factor 
is given by[27, p. 170] 


A = 


Pair cos e l f 
\ m Pabs cos Of + Pair cos Of 


(6.105) 


where 

A- COS 

m = — tr • 

pi cos 6 f 

Note that the factor (1/ sjp a ir) ^ as been incorporated 
account for the spatial attenuation from the line source. 


(6.106) 

into Equation (6.101) to 
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The problem remains to solve for the entry point (®i,t/p), given any point 
(yq > Vfi) in the aperture and a fixed source location (» a ,y a ). The equation of the 
incident ray can be expressed as 

Vp = ( x \ ~ x a ) tan B l f + y a (6.107) 

and the equation of the refracted ray as 

yq = (*2 - ®i)tan^ + y p . (6.108) 

Equating y q in Equations (6.107) and (6.108), one finds that solving the following 
non-linear equation 

(®1 - * a ) tan + (x 2 - ®i)tan 0f + ( y s - y q ) = 0 (6.109) 

with the relation between 0 ^ and 0^ given in Equations (6.20) and (6.22), will 
yield 0j and 0j, This equation can be solved iteratively by using an appropriate 
computer routine, and the entry point (®i,y p ) can then be found by substituting 
0 l j into Equation (6.107). 

0.3.3 Oblique incidence on a wedge illuminated by a line source 

Consider now the case where a ray enters the wedge part of the slab at an 
oblique angle (0^) at ( x p ,y p ) as shown in Figure 164. The subscript (w) refers to 
the wedge part of the slab. Note that 0 X W is different for each ray. Such a ray will 
illuminate the aperture in the region (y 2 < y q < y we ), where the ray intercepting 
the aperture at (x 2 ,ywe) enters the wedge at (.Ti,yj) as shown in Figure 164. 
The aperture field resulting from a plane wave illuminating the wedge is given in 
Equation (6.57). In this case, however, the wedge is illuminated by a line source. 
The total magnetic field in the aperture due to rays entering the wedge is thus 
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given by 


Hy pw (y q ) = - y o T X T 2 


'Po cos 6 a + ja o sin d a ' 


A e Q i y abs 


; 3{^Pair+0\ Pabs ) 
y/ Pair 


( 6 . 110 ) 

for y 2 < y q < Vwe where p a f, a and y a i a are given in Equations (6.59) and (6.60), 
and 


A = 


m = 


Pair — 


Pair COS 6 l w 


\ mpabs cos 0 l w + p air cos 9 l w 
k cos 0^ 


Pi cos 


, and 


\]{ X P - x s ) 2 + (y p - y s p 


( 6 . 111 ) 

( 6 . 112 ) 

(6.113) 


Note that the factor ( 1 / ^p a i r ) has been incorporated nto Equation (6.110) to 
account for the spatial attenuation from the line source. 

The problem still remains to find the entry point ( x p ,y p ), on the wedge given 
a point (x 2 ,y q ) in the aperture and a fixed source location (x,,y s ). Consider the 
ray that intercepts the wedge perpendicularly. Let the entry point for this ray be 
{ x iiVi)i and the point where the ray intercepts the aperture be (x 2 ,y ws ) as shown 
in Figure 164. The entry and intercepting points are given by 


x i = 

(xi 

V a 

y t = 

a(x s 

y ws — 

( x 2 ' 


a = 

x 2 — 


(lt + X ’ a + y,- yi ) 


y i - 1/2 


+ y 3 , and 


2/2 - 2/1 

Having determined the point (a C{,yi), one finds that 


(6.114) 

(6.115) 

(6.116) 
(6.117) 


ywe = (*2 - *i)tan(V> T «. ) if x i > X 1 
hf = unit vector normal to parallel slab (=— x). 
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(6.118) 

(6.119) 


ip = angle between and n/ (= wedge angle) 
6w = angle of ray transmitted into the dielectric 
relative to — (see Equation (6.22)) 
h w = unit vector normal to wedge, and 


To solve for (xp,t/ p ), one uses the fact that (xp,y p ) lies on the wedge, 

_ (jg-flKjfa-gl) 


Xp — 


1/2 - 2/1 


+ Xl 


or 


Xp — ot/p + b 

where a is given in Equation (6.117) and 

6 = x\ - ay\ . 

The equation of the incident ray is given by 

Vp = (xp - x a ) tan 7 n j + y s 

where 

Tnt = V’ T @w ^ Vq ^ Vws • 
The equation of the refracted ray is given by 


y P = y q - (*2 - Xp)tan7„* 


where 


< 


Unt — 4' t ^ y q yws • 

Solving for Equations (6.121 )—( 6.126), one finds that 

Vg ~ x 2 tan 7 n t + b tan 7 n < _ y s + b tan 7,,,; - x s tan j ni 


y p 


1 — a tan 7 n t 


1 - a tan 7 n j 


so that 

( 6 . 120 ) 


( 6 . 121 ) 


( 6 . 122 ) 


(6.123) 


(6.124) 


(6.125) 


(6.126) 


(6.127) 
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Equation (6.127) can be solved iteratively by using an appropriate computer rou- 
tine, if one keeps in mind that 0 l w and 6 l w are related by Equation (6.22). Once 
0 l w and Oh, have been found, (x p ,y p ) can be determined from Equations (6.127) 
and (6.120). 

The total magnetic field in that part of the aperture where the slab is present 
(y > y 2 ) is thus given by 


H y P (yq) ~ Hy P3 {y q )U(y q - y f t ) + H“ pw (y q )[U (y q - y 2 ) - U(y q - yu, e )] (6.128) 
where 


U(x) 


(6.129) 


1 if x > 0 
0 if x < 0 . 

Note that Hy P3 is given in Equation (6.101) and Hy PW is given in Equation (6.110). 
Equation (6.128) can be substituted into Equation (6.78) to give the field scattered 
from the slab ( E* d ). The total field is then given by Equation (6.73). 

6.3.4 Results 


In this section the forward scattered fields from a semi-infinite absorber slab 
with a wedge termination is calculated for three cases. Consider the geometry of 
Figure 164. It is assumed that an electric line source with a uniform pattern is 
located at ( x s ,y 8 ) = (—4.576', 1.574') and that the absorber slab is 0.5' thick; ie., 
x 2 — xi = 0.5'. Note that the position of the source corresponds to the location 
of the bottom edge of the subreflector in Figure 123, where it is assumed that the 
coupling aperture is 2' wide. The origin in this case is located at the apex of the 
wedge, rather than at the focal point of the subreflector. 

The total field given in Equation (6.73) is calculated for three slabs, desig- 
nated as Q 5 , Qq and Q7, with wedge angles (V’) equal to 7.1°, 14.0° and 17.4°, 
respectively. The absorber parameters are given in Table 21. Patterns are cal- 
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Table 21: Parameters of semi-infinite absorber slabs. 



<?5 

<?6 

Qi 


7.1° 

14.0° 

17.4° 

*1,2/1 

l 

0 

b, 
— » 

-0.5', 2' 

-0.5', 1.6' 

*2,2/2 

o',o' 

0 

0 

o',o' 

4 

1.5 

1.5 

1.5 

4' 

1 

1 

1 


ciliated at 3 GHz along a line parallel to the y - axis at a constant ^-displacement 
of 5'. The aperture integration is performed by Simpson’s method, with aperture 
spacings of 0.1A. The aperture field is integrated from y = y 2 to y = 11.87', at 
which point the aperture field ( E™ m ) is 45 dB below the level that the field would 
have been if there was no absorber present. 

The field patterns for the total field (see Equation (6.73)), field due to the 
end point contribution from (a7 2 ,t/ 2 ) (see Equation (6.81)) and the scattered field 
obtained from the aperture integration over the slab (see Equation (6.78)) are 
shown in Figures 166-168 for Q5, Qq and Qt, respectively. Note that the scattered 
field obtained from the aperture integration over the slab does not have a peak 
along the geometrical ISB ( y isb = -1.72' at x = 5'), but that the peak is shifted 
towards the interior of the slab. This effect was also observed in Section 5.5 where 
finite slabs were illuminated by a plane wave. It is due to the ray bending associated 
with propagation through the dielectric medium. 

The total fields for the three cases are compared in Figure 169. It is clear from 
this figure that the magnitude of the ripple in the total field decreases as the wedge 
angle (t/>) decreases, as expected. Note, however, that although the ripple decreases 
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as ip decreases, the total field in the region around the apex (yj a {, < y < B 1 ) 
increases as ip decreases. This is due to the fact that a smaller wedge angle results 
in a thinner slab in the region of the apex, so that the transmitted field in this 
region is stronger. In the region away from the apex ( y > 8*), one finds that the 
field values are the same for the three cases (except for the ripple). This is to be 
expected since the rays are travelling primarily through the parallel section of the 
slab, and all the slabs have an equal thickness. In the region y < y^,, one finds 
that the field values are also the same for the three cases (except for the ripple). 
This is to be expected, since the total field in this region is given by the direct and 
end point contributions of the integration over the free space half plane, with the 
diffracted field from the slab superimposed. Since the wedge with the largest wedge 
angle ( Qj;ip = 17.4°) has the strongest diffracted field, one would thus expect the 
total field from this slab to have the largest ripple in the region y > y; s j,. Similarly, 
one would expect the wedge with the smallest wedge angle {Q^\ip = 7.1°) to have 
the smallest ripple in this region. This is indeed the case, as shown in Figure 169. 

As indicated in Section 5.6, the wedge terminations of the coupling aperture 
of the absorber fence has to be designed such that diffractions from the absorber 
fence is minimized. It has been shown that a smaller wedge angle results in smaller 
diffracted fields. In addition, the transmitted fields through the absorber fence have 
to be attenuated sufficiently. A smaller wedge angle leads to a thinner slab near the 
apex, which result in stronger transmitted fields. The design of the absorber fence, 
and in particular the wedge angles of the coupling aperture terminations, will thus 
be based on a compromise between small diffractions and weak transmitted fields. 
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Figure 166: Scattered fields from Q 5 (3 GHz). 
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Figure 167: Scattered fields from Qq (3 GHz). 
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Figure 168: Scattered fields from Q'j (3 GHz). 
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Figure 169: Scattered fields from semi-infinite absorber slabs. 
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CHAPTER VII 


SUMMARY AND CONCLUSIONS 


This study has shown that the performance of a compact range can be enhanced 
by improving the design of the main reflector and incorporating a dual chamber 
configuration with a Gregorian subreflector feed system. The quality of the plane 
wave illuminating the target zone is improved with respect to a reduction in the 
amplitude taper, ripple and cross-polarization errors. 

It has been shown that the diffractions from the edges of the parabolic main 
reflector are reduced if rolled surface terminations are added to the parabola. First 
generation rolled edges were elliptic in shape, resulting in a continuous surface slope 
across the parabola/rolled edge junction, but it has a discontinuity in the surface 
radius of curvature across the junction. The edge diffracted fields can be reduced 
even more by smoothly blending the rolled surface terminations from a parabola 
into an ellipse. Blended rolled edges ensure that the radius of curvature as well 
as the slope of the surface are continuous across the junction. By choosing an 
appropriate blending function, one can design the rolled edges such that the first 
and higher order derivatives of the radius of curvature are continuous across the 
junction as well. The impetus for obtaining a smoother junction is that a smoother 
junction results in a more continuous reflected field \yhich leads to a weaker edge 
diffracted field and ultimately less ripple in the total field. A relation between 
the n-th order blending function and the corresponding error term that results in 
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a discontinuity of the n-th order derivative of the radius of curvature across the 
junction was derived, and it was shown how the blending parameters should be 
chosen to minimize the error term. An algorithm for designing optimized blended 
rolled edges was consequently developed, based on an analysis of the error term. 
Several two-dimensional design examples were discussed, and the reduction in edge 
diffracted fields resulting from blended rolled edges were clearly illustrated. 

The edge diffracted fields from three-dimensional reflectors can further be reduced 
by shaping the contour of the junction between the parabolic section and rolled 
edge contours. It was shown that a concave edge contour has superior performance 
over conventional convex and rectangular edge contours. Design equations for a 
complete three-dimensional main reflector incorporating blended rolled edges and 
a concave edge contour were developed. Several three-dimensional main reflectors 
were analyzed using a new corrected physical optics method; whereby, the end 
point contributions from false scattering centers on the incident shadow boundary 
of the main reflector are removed. The results confirmed that the edge diffractions 
from a three-dimensional main reflector can be significantly reduced if the reflector 
has a concave edge contour and blended rolled edges. 

A dual chamber compact range configuration was proposed; whereby, the main 
reflector and target zone are contained in one chamber and an oversized Gregorian 
subreflector with feed assemblies in the second, smaller, chamber. The chambers 
are isolated from one another by an absorber fence, with a small coupling aperture 
to transmit signals between them. It was shown that the use of a subreflector sys- 
tem results in a significant reduction in the amplitude taper and cross-polarization 
errors in the illumination of the target zone. In addition, the absorber fence vir- 
tually eliminates the effects of diffraction from the subreflector and spillover from 
the feed on the desired illumination. 
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In order to minimize the secondary diffractions from the absorber fence, the cou- 
pling aperture should be shaped. It was shown that a wedge termination is most 
appropriate, and that the selection of a wedge angle depends on a compromise 
between weak diffracted and transmitted fields. In order to analyze the scattering 
from the coupling aperture, a new technique to calculate the forward scattering 
from a homogeneous, lossy absorber slab with a wedge termination was developed. 
The technique combines the tracing of rays through the absorber with aperture in- 
tegration. Several finite slabs with wedge terminations and plane wave illumination 
were analyzed, and the results compared with those obtained from the moment 
method. It was found that the new aperture integration technique can predict the 
scattering from the absorber slabs within an acceptable accuracy. The technique 
was then expanded to enable one to determine the forward scattering from a semi- 
infinite absorber slab with a wedge termination and line source illumination. The 
results confirmed that a smaller wedge angle results in less ripple in the scattered 
field but also in a stronger transmitted field. 

Several topics for further research have become apparent during the course of this 
study; viz., 

1 . Investigation of ways whereby the physical optics integration of three-dimensional 
main reflectors can be performed faster and cheaper. 

2. The extension of the equations for the derivatives of the blended rolled edges 
to three dimensions, and the optimization algorithm to three-dimensional 
reflectors. The present algorithm analyzes one plane, which enables one 
to design optimized circularly symmetric reflectors. A three-dimensional 
extension would enable one to optimize in three-dimensions, which would 
be particularly useful for designing optimized reflectors with concave edge 
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contours. One would probably find that the blending parameters would have 
to change continuously around the perimeter of the edge contour in order 
to obtain an optimized three-dimensionsal design; ie., a “blending of the 
blending parameters”. 

3. The optimization algorithm can be automated, so that a computer can deter- 
mine the optimized blending parameters with minimal human intervention. 

4. The development of an optimized shaping of the concave edge contour of the 
main reflector. 

5. Improvement of the aperture integration technique for determining the scat- 
tering from a lossy dielectric slab with a wedge termination. Higher order 
internal reflections can be taken into account, especially near the apex of the 
wedge. The technique can be extended to include multilayer cases and more 
complicated geometries. 

6. Development of an analysis of the dual chamber design which includes the 
effects of the coupling aperture diffractions in order to predict the field quality 
in the target zone. 
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APPENDIX A 


APPLICATIONS OF DIFFERENTIAL GEOMETRY 


In this appendix, differential geometry [69] is used to derive expressions for 
the unit vector normal to a 3-dimensional surface, the radius of curvature and 
tangential vector of a 2-dimensional curve as well as the radius of curvature and 
its derivatives for a blended rolled edge. 

A.l Unit vector normal to a surface 


Let a surface in three dimensions be given by the following parameterized 
equation: 

r(u,v) = [x(u,v), y(u,v), z(u,v )] . (A.l) 

The partial derivatives of r are given by 


r u = 


r„ = 


dx dy dz 
du' du' du 
dx dy dz 
dv' dv' dv 


, and 


The unit vector normal to the surface can then be expressed as 

n(u, v) 


r« x rt, 


(A.2) 

(A.3) 


(A. 4) 


I r„ x r„ | ' 

Depending on the application, one may be concerned with the inward or outward 
normal to the surface. They are related by 


n 


in 


' ^out 


(A.b) 
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Figure 170: Paraboloid. 


In the case of bodies of revolution, the inward normal points towards the axis 
of symmetry, and the outward normal away from it. Two special cases are now 
considered; viz., a paraboloid and a spheroid. 

A. 1.1 Paraboloid 


Consider a paraboloid with the z-axis as its axis of symmetry as shown in 
Figure 170. This surface can be expressed as 


r(tt,u) = 


u, v, 


u 2 + v 2 


(■ AA ) 


where 


u — x displacement 
v = y displacement , and 
f c — focal length . 
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In this case one then finds that 


r u 

*v 


1, o, 


u 

Wc 


, and 


0 , 1 , — 

2 f c 


so that 


(A.7) 

(A.8) 


**in{ u i v ) — 


—U 


—v 


2 f c 


1 \ju 2 + V 2 + 4/2 ’ sju 2 + V 2 + 4/2 ’ yju 2 + v 2 + 4/2 J 
Notice that n ln points towards the point (0,0, 2/ c ). 

A. 1.2 Spheroid 


(A. 9) 


A spheroid is an ellipsoid with two axes lengths the same, such that it is a body 
of revolution. Consider now a spheroid with the z-axis as its axis of symmetry. Let 
the axis length along the direction of the z-axis be a e , and along the direction of 
the x- and y-axes be b e . This configuration is shown in Figure 171. The spheroid 
can be described by the following parametric equation: 


r(u, v) = [ fe e cos v sin u, b e sin v sinu, a e cos u 


(A. 10) 


where 0° < v < 360° and 0° < u < 180°. Since the spheroid is a body of revolution 
around the z-axis, the parameter v is the spherical angle (<f>). The parameter u 
corresponds to the spherical angle (0), and they are related by 


tan 6 = — tanu 

tie 


In this case then, 


(A. 11) 


r u - [ b e cos v cos u, b e sin v cos u, -a e sin u ] , and 
r t> — [ — b e sin v sin u, b e cos v sin u, 0 ] 
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(A. 12) 
(A. 13) 




Figure 171: Spheroid. 


so that 




— cos v sin u — sin v sin u —cos v 


be b e cte 

The unit vector normal to the surface is then given by 


n 


m 


n; 


in 


(A14) 


M- 15) 


Notice that in the case where a e — b e , the spheroid becomes a sphere. In this case 
nout is equal to the spherical radial unit vector, and u = 0. 

A. 2 Radius of curvature and tangential vector of a 2-dimensional curve 


Consider a curve in the yz-plane, given by the parameterized equation 


r(7) = [0, y( m), =( 7 ) ] • 


(-4.16) 


The curvature (k) at 7 is given by 


I r '(7) x r" ( 7 ) | 
I r'(7) | 3 


(-4.17) 
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where 



r ( 7 ) = = [ 0 , y ( 7 ), z'( 7 ) ] , and 

(A.18) 


r (t) = d ^2 = i°» y (7)» 2 ( 7 )] • 

(A.19) 

The radius of curvature ( R c ) at 7 is then given by 



R c(7) = -y-r 
«(t) 

(A. 20 ) 

so that 

1 r' I 3 

JWT >“ |HxU| • 

(A. 21 ) 

However, 

| r' X r" |=| y'z" - y"z' \ 

(A. 22) 

so that 

Kv ') 2 + 

^ c(7) \y'z"-y"z'\ ■ 

(A. 23) 

The vector tangential to the curve is given by 



Hi) = r'(l) 

(A. 24) 

so that 

i= it] = 

(A. 25) 


and since n • t =0, one finds that 


n = (* 2 ,-<l). (-4.26) 


Two special cases will be considered; viz., an ellipse and a parabola. 
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Figure 172: Ellipse. 


A. 2.1 Ellipse 

Consider the ellipse described by 

r (7) = [0, b e sin 7, a e C0S7]. {A. 27) 

The semi-major and semi-minor axes lengths are a e and b e , respectively. Notice 
that the angle 7 is not a real angle, but a parametric one. In this case one, finds 
that 


so that 


y — b e cos 7 
y" = -b e sin 7 
z' = — o e sin 7, and 
z" = — o e cos 7 

t(7) — [ 0, cos 7, — a. € sin 7 ] , and 


(A. 28) 
(A. 29) 
(A. 30) 
(A. 31) 

M-32) 
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(>1.33) 


# c ( 7 ) = (aeSin 2 7 + 6^cos 2 7)3 
Qe^e 


At 7 = 0 and 7 = 90°, one finds that 


b 2 

Rc( 0) = ^ 


a: 


H c (90°) = ^ , and 
be 


^£u\o 

d'y 


(0°, 90°) = 0. 


(A.34) 
(A. 35) 
(A. 36) 


Note that it is assumed that a e > b e so that i2 c (0) < # c (90°). This is also evident 
from Figure 172. 

A. 2. 2 Parabola 


Consider the parabola described by 


r( 7 ) = 


0, 7> 


7 

4/c 


(A.37) 


In this case 7 is the height ( y ) shown in Figure 173, and f c is the focal length. 
Thus, one finds that 


so that 


2/'(7) = 1 


y "{ 7) = 0 


Al) 

z< \l) 


7 


2 f c 

1 

Wc 


, and 


t (7) = [°» Jr] » ancl 

^ JC 


Rc{y) — ‘Lfc 




(A. 38) 
(A. 39) 
(A. 40) 

(A.41) 

(A. 42) 
(A- 43) 
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Figure 173: Parabola. 
The derivatives of R c with respect to y are given by 


R'Ay) = 

+ (y/2/c) 2 


(A.44) 

R"(y) = 

3 y 2 

+ + (y/ 2 /c) 2 

6jC 

(A. 45) 

s/cV 1 + (y/ 2 M 2 

R'c(y) = 

9 y 


(A.46) 

8/cV 1 + (J//2 /c) 2 

32/|yi +W 2/c) 2 

Rc'(y) = 

9 

9y 2 

(A.47) 

8/c\/l + (y/2/c) 2 

16/1^1 + (y/2f c ) 2 


128 f c fl + (t,/2/e) 2 ' 

A. 2. 3 Blended rolled edge 

The blended rolled edge is described in Section 3.4. In this section expressions 


will be derived for the radius of curvature ( R . c ) and the first four derivatives of R c 
in the yc-plane. 




It is useful to evaluate the derivatives of R c of the blended surface with respect 
to y rather than 7. This serves as a standardization for comparison of different 
types of blending, and makes it convenient to compare the derivatives of R c of the 
blended surface at the junction with that of the parabola (see Equations (A.44)- 
(A.47)). In addition to deriving the general expressions for the radius of curvature 
of the rolled edge (R c ) and its derivatives for all 7, the derivatives of R c at the top 
junction (7 = 0) will also be evaluated as a special case. This is necessary in order 
to examine the continuity of the derivatives of R c across the junction. A similar 
analysis can be performed at the bottom junction. In the rest of this section a 
prime will denote differentiation with respect to 7, and (yj,Zj) will indicate the 
coordinates of the junctions. 

The general two-dimensional blended rolled surface is expressed in the yz- 
plane by Equations (3.73)-(3.75) such that 



r biend(i) = y(i)y + ^(t)* 

(A.48) 

where 



y{i) = 

(By + 2/j ) [ 1 - 6(7)] + [C sin 7 + Z?(l - cos 7) + yj]b(y) 

(A. 49) 

2(7) = 

(By + yA 2 

4fc * H7)] + [^sin 7 + F(l-cos7) + z#(7) 

(A. 50) 

where 

p _ i(*m/7in) 

V 1 + (Vj/2fc) 2 

(A- 51) 


n _ ±a e 

\/l + iVj/ifc) 2 

(A. 52) 


D - , b Jh 

2/c\/l + (*r>/ 2/c) 2 

(A.53) 


E - , and 

2/c\/l + (Vj! 2/c) 2 

(A. 54) 
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F = , e ■— • ( A - 55 ) 

\/l + (j/j/2/c) 2 

The (±) sign refers to the top and bottom junctions respectively. The following 
expressions are useful in evaluation of the derivatives of R c at the junctions: 


C yj 

= 0 

(A. 56) 

2 f c 



Dyj 
2 f c 

= -b eS /l + (yj/2f c ) 2 , and 

(A. 57) 

Eyj 
2 fc 

= ±u e \/i + (yj/2/c) 2 

(A. 58) 


where the (±) sign refers to the top and bottom junctions respectively. 


1 ) Radius of curvature 

The radius of curvature is given by Equation (A. 23) such that 


where 

R»(t) 

R 1(7) 

In this case, one finds that 

1/(7) = b(C cos j + Z)sin7 - B) 

+b'(C sinj + D( 1 - cos 7) - Bj) + B 
y"( 7) = b(D cosj — Csin7) + 2b'(Ccosj + I? sin 7 — B) 

+b"(C sin 7 4- D( 1 - cos 7) - £7) 

/ 5(57 + yy)\ 

z'(7) = 6 ( Ecosj + F sin~f ) 


He(7) = 


Ruh) 
Ri( 7) 


= ^') 2 + (*') 2 > and 


= I y'z"-y"z' 


(4.59) 

(A.60) 

(A.61) 

(A. 62) 
(A. 63) 
(A. 64) 
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+b' [ jEJsin 7 + F (1 — COS 7 ) + zj 


( By + yj) 2 \ B[By + yj) 


z"(y) — b ^F cosy — E siny — 

'( 


4/c 


*) 


+ 


2 f c 


B(By + yjY 

2 fc , 


+26 cos 7 + F sin 7 

+ 6 " sin 7 + F(1 - cos 7 ) + zj - — ~-^ Vj — j 

At the junctions (7 = 0) one finds that 

2 /( 0 ) = 5 
2 /"( 0 ) = 26'(<7-fl) 

z '(°) = Yf~ » and 


2 » = *(*-f*) + g 


, and 
(A. 65) 


(A. 66 ) 
(A.67) 
(A. 68 ) 

(A. 69) 


so that 


fl«(0) = |5| 3 ^/l + ( yi /2/ c )2 ,and (A. 70) 

^0) - I £ + **(*- f£) I- < A ' n > 

Substituting Equation (A. 56) into the equations above, one finds that 

R b J end = 2 f e ^1 + (y J /2/ c )2 = i2P“™ fcoia . (A. 72) 

The use of a blending function will thus ensure that the radius of curvature is 
always continuous across the junction. Note that is this true for the top and 
bottom junctions. 


2) First derivative 

The first derivative of R c with respect to y is given by 

dR c dR c dy 1 dR c 

dy dy dy y' dy 


(A. 73) 
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where 


dRc R{R ^ — R![Ru 

+" = Wi? 

In this case, one finds that 


<(7) 

= 3^1 (y 1 ) 2 + (z') 2 {y'y" + z z") 

(A. 75) 

-Rf(T) 

= ±(yV" - z'y'") if y'z" < zy" 

(A.76) 

y"'(7) 

= —b(C C0S7 + D sin7) + 36 , (D COS7 — C sin7) 

(A.77) 


-f36"(<7 cos 7 + D sin 7 — B ) 

+b"'(C sin 7 + D(1 — cos 7) — B7) , and 



. / 5 2 \ 

(A.78) 

^"(7) 

— — 6(Fsin7 + 2? cos 7) + 36 1 F cos 7 — E sin 7 2/ ) 


+36" cos 7 + F sin 7 - 

+6 W sin 7 + F(1 - cos 7) + Zj - ) • 


In order to examine the continuity of the first derivative at the top junction, it is 
neccesary to evaluate dR c /dy at 7 = 0. Assume that 6(7) is a first order blending 
function, so that 6(0) = 0 and b'(0) 7^ 0. One then finds that the following holds 


at the top junction: 

y'(o) = B (A. 79) 

y"(0) = 2b'(C-B) (A. 80) 

y'"(0) = 3b'D + 3b"(C - B) (A. 81) 

(A.82) 

z"( 0) = 26' ( £ - ff ) + ^ ■ antl (*.83) 

= 36 '( f -fi) +3t "( E - 3 • <A - 84) 
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Since B > 0 at the top junction, it follows that: 


One thus finds that 

Ri( o) = 


2 f c 


Rlivtop) = y'z" - y"z' . (4.85) 

(A. 86) 

*i(0) = -3b' B + ( Vi /2f c )l + Si) (A.87) 

i2 M (0) = B 3 + ( Vj/2f c ) 2 , and (A.88) 

^(°) = 6b' B 2 (M 1 + (yj/2f c ) 2 } - B tyl + ( yj/2f c ) 2 ^j (A. 89) 

, 3571 + ^/2/0)2 

4/1 


so that 


(dRc\ 

W y) 


blend 


( 7 = 0 ) = 


3 Vj\J 1 4- (yj/2f c ) 2 


2/c 


(A. 90) 


+ (»e[l + (»j/2/ c ) 2 ] 

+^(1 + (%/ 2/ c ) 2 ] 2 - f ^1 + ( y ,/ 2 / c ) 2 ^ . 

Using Equation (A.44), it follows that the following is true at the top junction: 


(dR c \ blend ( dR c \ parabola , 

W “Uv +b ^ 


(A.91) 


where 


12/c[l + (!/,/2/c0 2 I 2 / a e /ci+l+(»,/2/ c ) 2 A 

(W7m) ((*m/7m) (xm/jm) 2 2/ ' f' 1 ' 92 ) 


$ / Second derivative 

The second derivative of R c with respect to y is given by 

d 2 R c __ l d_ ( dR c \ 
dy 2 y' dry \ dy I 


(4.93) 
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where 


so that 


d_ (dRA = d_ ( dRc A = y'{d 2 R c /d-y 2 ) - y"{dR c /dy) 
dj \ dy ) d-y \ d-y y' ) (y 1 ) 2 

d 2 R c = y'(d 2 R c /dy 2 ) - y"{dR c /d-y) 
dy 2 ( y') 3 


It is found that 


where 


d 2 fl c = - R'/Ru) - 2 R\{R l R! u - fl'^) 

d-y 2 (R { ) 3 


KM = 3 v /(w ')2 + (z')W' + (/) 2 + (^") 2 + z'z'"] (A.97) 

+ %V+ i VV 

v ^( y ') 2 + ( z ') 2 

= ±(»'^ + y"z'" - /'z" -»‘V) if j,V' < j,V (A. 98 ) 

y tv ( 7) = (6 — 66^ )( (7 sin 7 — Dcos-y) — 4 b' (Ceos 1 + Dsin7) (A. 99 ) 

4 - 46 w (C cos 7 + D sin 7 — B) 


+b lv (C sin 7 + D{\ — cos 7) — B-y) , and 
2*^(7) = 6(£^sin7 — FCOS7) — 46 , (i r sin 7 + £Jcos7) 
+66” ^Fcos 7 — is 1 sin 7 — 

+ 4r(£cos 7 + Fsin 7 - g^Mj 

-LA” 1 ^ T? ei n ~ -L EV1 _ = .7 (^7 Vj) 


(A. 100 ) 


+6 ,l; I E sin 7 + F (1 — cos 7) + zy — 




Consider now 6(7) to be a second order blending function, so that 6 ( 0 ) = 0, 
6 , ( 0 ) = 0 and 6 ”( 0 ) / 0 . One then finds that the following holds at the top 
junction: 


y'(0) = b 


(A. 101) 
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y» = o 

y'"( 0) = 3b"(C - B) 
y™(0) = 6b" D + 4b'"(C - B) 

2 ' (0) = t 

B 2 


z"(0) = 


2 f c 


z"'(0) = 3 b" j , and 

**■(«> - 6h " ( F - £) + 46 "' i E - Tf c 


(A. 102) 
(A. 103) 
(A. 104) 
(A. 105) 

(A. 106) 

(A. 107) 

(A. 108) 


so that 


B 3 

w) - Wc 




0 , and 


(A. 109) 
(A. 110) 


*?(0) = -6 b"B (b'Jl + (vj/2/ c )* + |r) - W'irJ 6 ’ - B ) ■ ( A - m > 

Since y M (0) = 0 and i?j(0) = 0, one finds that 

d?R c - R','R U 

^ h=0) - (W (AU2) 

where 


fitt(O) = B 3 yi + (</j/2/ c ) 2 , and 

3BVl + (»i/Vc) 2 , 3i»V 


4/1 


16/cV 1 + (w/2/c) 2 

+9B 2 6" («.[! + (» > /2/«) 2 ] - B^l + (»y/2/ c )») . 


(A. 113) 
(A. 114) 


Combining the equations above with Equation (A. 45), one finds that the following 
is true at the top junction: 


d 2 R c \ blend _ (d 2 R c y arabola ( Lll 

dy 2 j 


= (^) 
\ / 


+ b"(0)e 2 


(A. 115) 
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where 


24/c ^/l + (t/j/2/ c ) 2 / Qe / c 6 e yi + (yj/2/ c ) 2 l\ 

(*m/7m) 2 \(®m/ 7 m) + (®m/7m) 2 2y ‘ 


Third derivative 

The third derivative of iZ c with respect to y is given by 



d 3 .R c 1 d 1 d 2 R c \ 

dy 3 y' d 7 \ dy 2 / 

(A. 117) 

where 

d /d 2 fl c \ d / y , (d 2 H c /d 7 2 ) — y"(dR c /d'y)\ 

d7 V dy 2 / d 7 \ (y') 3 / 

(A. 118) 

so that 



d 3 R c 

(y') 2 (d 3 J R c /d 7 3 ) - 3y'y n (d?R c / d 7 2 ) + [3(y") 2 - y'y"' 

](di? c /d 7 ) 

dy 3 

(y') 5 

(A. 119) 

It is found that 

d 3 i? c i 2 ? 3 

(A. 120) 


where 


(A. 121) 


(A. 122) 


= (RiflRiR" + R\Rl - R"R' U - R'i'R,,) 

-4RiR'[(r,Ru - n','R u ) 

+[6(Bi) 2 - - R',R U ) 

K(l) = 3\As/') 2 + (z') 2 (</V' + 3 < j/'V" + --V") + 

9(j/'/ + zVW' + (y") 2 + (--") 2 + .-'.-"'I 

\/W) 2 + (~) 2 

3 (!jV + zV ') 3 

y(y ) 2 +(=') 2 
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R'l'\ 7) 

y v ( t ) 


^( 7 ) 


±[yV + 2(y" z* v - z"y iv ) - z'y v ] if y'z" < y" z' 
(b — 10b" )(C cos 7 + D sin 7 ) 

+(10fe m — 5b')(D COS 7 — C sin 7 ) 

+5b iv (C cos 7 + D sin 7 — B) 

+ 6 v (Csin 7 + D( 1 - cos 7 ) — B 7 ) , and 
(b — 10 b")(E sin 7 + F cos 7 ) 

+56 , ( E sin 7 — FC 0 S 7 ) 

(■ 

+56™ ( E cos 7 -f - F sin 7 


B 2 \ 

+106 m ( F cos 7 — E sin 7 — — - ) 

2/c/ 

gCgr + yj) 


b v (esi 


2/c 


+ 6 W I E sin 7 + F(1 - C0S7) + Zj — 


{Bj_ + 

4/c 


5) Fourth derivative 

The fourth derivative of R c with respect to y is given by 

d A R c _ £ (<PRc \ = I £ (d*Rc\ 
dy 4 dy \ dy 3 / y 1 dy \ dy 3 / 

where 

A / \ _ jd_ / (y') 2 (d 3 R c /dy 3 ) - 3y'y"(d 2 R c /dy 2 ) + [3(y ") 2 - y' 
d 7 V dy 3 ) dy \ (y 1 ) 5 

so that 

d^ f2 c i? c 4 

~S+ = (7F 

where 

R ci = - Hy’f^Rc/dj 3 ) 


(A. 123) 
(A. 124) 

(A. 125) 


(A. 126) 

y ,n ] (dR c /dy) 
(A. 127) 

(A. 128) 
(A. 129) 
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+ [15y'(/) 2 - 4(y') 2 y"'](d 2 i? c /d 7 2 ) 
+[iOyW" - (y') 2 y iv ~ \h{y")\dR c /d 1 ) . 


Furthermore, one finds that 

d*R c _ Rgi 
d 7 4 ( Ri ) 5 

where 

R gi = (RifiRiR™ + 2(i?}< - r!('r! u ) - *{«£„] 
-GiRtfR'iiRM 1 + R'X - R'{R' U - R'I'Ru) 
+[lSRi(R' l ) 2 - 6(J^) 2 J2j')(/fcj£ - R'/Ru] 
MlSR^R'l - 2(R l ) 2 R' l " - 24(R l f](R l R' u - R' t R u ) 

K(i) = Sy/iv 1 ) 2 + { z ') 2 [y'y v + My"y iv + z"z iv ) 

+3([y'"] 2 + [z" f) + zV] 

< ^[y'y™ + 3(yV» + z"z"') + z'z iv }(y'y " + z'z") 

\f(y') 2 + (z') 2 

, 9[y'y , " + (y") 2 + (z") 2 + z'z"'} 2 

y/(y') 2 + {z') 2 

( y'y " + z'z") 2 [y'y'" + ( y ") 2 + (z") 2 + z'z'"] 
l/(y') 2 + (z') 2 

| 9 (y'y" + z'z")* 

fa n 2 + (*') 2 

*T(t) = ±[y'z vi + 3(yV-zV ) 

+2(y"V t -/'y i VzV i ] if y'z"<z'y" 
y m ( 7 ) = (6 — 156” + 15fc*")(£> cos 7 — C'sin 7 ) 

+(6i/ — 20b'" )(C cos 7 + Z?sin 7 ) 

+66 V (C cos 7 + D sin 7 — B ) 

+6 V *(C' sin 7 + £)(!— cos 7 ) — Bj) , and 


(A.130) 
(A. 131) 

(A. 132) 


(A. 133) 
(A. 134) 
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z vi ( , y) = (6 - 15b")(F cosy - Esin 7) 

+(66* — 2Qb n, )(Ecosy + Fsin7) 
+15 b iv ^Fcos 7 — E sin 7 — 


(A. 135) 


k ’( 


+66 v E cos 7 + Fsin7 — 


B(B 7 + j/j)' 


2/c 


+6 vt I i?sin7 + F(1 — cos 7) + Zj 


(Bl + yjf 

4/c 


Consider now 6 ( 7 ) to be a fourth order blending function, so that 6 ( 0 ) = 0 , 6 / ( 0 ) — 
0, 6 H (0) = 0, 6"'(0) = 0 and 6 * v ( 0 ) ^ 0. One then finds that the following holds at 
the top junction: 


y'( 0 ) 
y>) 
y"'(0) 

y tv ( 0 ) 
y v ( 0 ) 
y vt (0) 

z'(0) 

z"(0) 

z "'( o) 

z* v (0) 

z v (0) 

z vi { 0 ) 


= B 
= 0 
= 0 
= 0 

= 5b iv (C — B) 

= 156* V D + 66”(C' — B) 
Byj 
2 f c 

B 2 

2/c 
- 0 


= 0 



(A. 136) 
(A. 137) 
(A. 138) 
(A. 139) 
(A. 140) 
(A. 141) 
(A. 142) 

(A. 143) 
(A. 144) 
(A. 145) 
(A. 146) 

(A. 147) 
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so that 


*'<°> - v c 

R'ii 0 ) = 0 
R'l( 0) = 0 
R'i'( 0) = 0 , and 


Rj v (0) = -15Bb tv lb e] /l + (yj/2fc) 2 +2f~ 


(A.148) 


(A. 150) 
(A. 151) 
(A. 152) 


It is then found that 


-lSnC-B) — 


d*R C/ n . B|Bj l^R£Ru 

iy 4 (7 “ ' (*|)V) 4 


(A. 153) 


where 


R u { 0) = B 3 ^/lV(y j /2fc) 2 ,and (A.154) 

QR 7 

<«(0) = , . - (A. 155) 

Wts/l + tij/Vc) 2 

9B 7 y 2 

32/6 3/1 + ( yj /2/ c )2 
, 9-B 7 y 4 

+ 256/8^1 + ( yj /2/ c ) 2 

+15B 2 6* v (a e [l + (yj/2/c) 2 ] - 5 ^1 + (yj/2/c) 2 ) • 
Combining the equations above with Equation (A. 47), one finds that the fol- 


lowing is true at the top junction: 
/ >4 \ Wend 


(d*R c y len _ (d A R c \ 
\ dy 4 j \ dy 4 / 


parabola 


+ 5 l, ’(0)e 4 


(A. 156) 


where 


60/cv^l + (yj/2/c) 2 / a e / c 6 e ^1 + (yj/2/c) 2 1 


{Em/lmY 


( %mllm ) 


{xmllmY 
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APPENDIX B 


CLOSED FORM EXPRESSIONS FOR THE APERTURE 
INTEGRATION OVER A HALF PLANE IN FREE SPACE 


In order to solve the problem described in Section 6.3, it is required to de- 
termine the scattered field resulting from the aperture integration performed over 
a planar aperture extending over the half plane (-00 < y < 0 ) in free space as 
shown in Figure 174. The two-dimensional transverse magnetic (TM) scattered 
field at a field point (x,y) is then given by Equation (6.77) as 

E?[z,y) = iy' (BA) 

where y 1 is a point in the aperture, and 

P = \/ z 2 + (y 1 - y) 2 • (5.2) 


Note that the angle ( <f > ) which p makes with the — y-axis is given by 

<t> = arccos ( 5 . 3 ) 

as shown in Figure 174. Consider an electric line source to be located at (x s ,y s ). 
The electric field incident upon the aperture from this source is given by 


where 




(5.4) 


P = \J* 2 s +(.</'- Vs) 2 


(5.5) 
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Figure 174: Aperture integration over a half plane in free space. 
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so that the tangential magnetic field in the aperture is given by 

, — sin 


where 


K p (y') = 


<j>' = arccos 


Z 


Vs - y' \ 

p' ) 


The scattered field can thus be expressed as 






P 1 P 


(5.6) 


(5.7) 


(5.8) 


The integral in Equation (B.8) can be evaluated by the method of stationary phase 
as explained in Chapter 2, so that 


5f « E% p + E% p (5.9) 

where 5®“ p is the contribution from the stationary point, and 5*“ p is the contri- 
bution from the end point (at y — 0). The end point at y = — oo is too far away 
to have any effect. 

Following the notation of Equations (2.79)-(2.83), one finds that 


or 


It follows that 


or 


$(</) = -(p +p) 

(5.10) 

•Hy 1 ) = + W - y> ) 2 + y* 2 + (»' - y) 2 ) ■ 

(5.11) 

*v> = -(*'-,* + »' 7 ? ) 

(5.12) 

V p p / 


$'{y') — cos <t> — cos <j> . 

(5.13) 
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Figure 175: Stationary point in the aperture integration over a half plane 

space. 

The stationary point ( y ' sp ) occurs where 

*(y'sp) = o 

so that y' sp is that point in the aperture where 

COS $ = COS <f>sp 

as shown in Figure 175. 

Furthermore, one finds that 



so that 


\$"{y' sp )\ = sin V (^7^) ’ and 

sgn[$"(y' p )j = -1 


in free 

(5.14) 

(5.15) 

(5.16) 


(B.17) 

(B.18) 


where p and p' are interpreted as p — p(y' = y ' ap ) and p' = p'(y' = y' ap ). 

Substituting Equations (B.17) and (B.18) into Equation (2.79), one finds that 
the stationary point contribution to the scattered field is given by 

e ~jk(p'+p) 


ti sa 

&Z3P = 


\Jp' +P 

Equation (B.19) can thus be expressed as 


(5.19) 


El%{x,y) 


0<t'<<t>'o 

otherwise 


(5.20) 


where 


Ps = \f(x - x 3 ) 2 4 (y - y $ ) 2 , and (B.21) 

<Ao = arccos ^ — ^ . (B.22) 

This term represents the incident ray from the source at ( x s ,y 3 ) to the field point 

at ( x,y ). Note that E* ap is equal to zero when the stationary point ( y ap ) is no 

longer in the aperture. 

Consider now the end point contribution (E*“ p ) due to the end point at y = 0. 
Substituting Equation (B.15) into Equations (B.13) and (2.83), one finds that the 

end point contribution can be expressed as 

' 


E» a 

^ zep 


e ~j k p' , e ~i k P 


y/f 


where 


D a (<j>{ h 4>o) = 


xTp 


e • ?,r / 4 sin <f>Q 


(5.23) 


\/27r A:( cos (j > (j — cos </>o ) 
\/ I Us , and 


<t> 0 


p yj x 1 4 y 2 , and 

7 v\ 


arccos 


Pfi 


(B.24) 

( 13 . 25 ) 

(B.2G) 

(B.27) 
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Note that Equation (B.23) is cast in the form of a ray diffracted from the end 
point. 

However, the diffraction coefficient ( D a ) becomes singular at the shadow 
boundary; ie., where 0' o = 0 o- This singularity is similar to the one that occurs in 
Keller’s GTD diffraction coefficients [16]. The discontinuity in the diffraction coef- 
ficients cause the total field (ie. incident plus diffracted fields) to be discontinuous 
at the shadow boundary. Koumyoumjian and Pathak developed the UTD[17] to 
correct this problem, and found that the total field becomes continuous when the 
diffraction coefficients are multiplied by a transition function ( F ). The coefficient 
in Equation (B.24) can now be heuristically multiplied by the same transition 
function, so that modified diffraction coefficient ( D ep ) is given by 


i e sin <j>' Q 

D p {4>qi <t>o,L) = -■ °- 7T- 

V27rfc(cos0Q — cos 0o; 


kL 

( COS 4> Q — COS 00 \ ^ 

2 1 

\ sin< £o / 


(B.28) 


where the transition function is described in Section 2.2 and 


PoPO 
p'o + PO 


{B. 29) 


The end point contribition to the total scattered field can thus be expressed 


as 


e -jV , e ~jkp 

E s Z ep — f=f~D ep (<f)Q,<f>o, L)- 


(B. 30) 


/p' y/P 

The inclusion of the transition function causes the total field (E* a ) given in Equa- 
tion (B.9) to be continuous. 
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